UNCLASSIFIED 


_ AD  NUMBER _ 

AD873169 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited.  Document  partially  illegible. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
and  their  contractors;  Critical  Technology;  JUN 
1970.  Other  requests  shall  be  referred  to  Air 
Force  Weapons  Laboratory,  WLCT,  Kirtland  AFB, 

NM  87117.  Document  partially  illegible.  This 
document  contains  export-controlled  technical 
data . 


_ AUTHORITY 

afwl,  ltr,  30  nov  1971 


THIS  PAGE  IS  UNCLASSIFIED 


IEUillcopy  AD873169 


AFWL-TR-70-20 

a(\ 


t 


f  / 


AfWl-T*- 

70-20 


ANALYTICAL  MODEL  FOR  HIGH  EXPLOSIVE 
MUNITIONS  STORAGE 

H.  L.  Schrayar 
L.  E.  Romasbarg 


Mtchonici  Rtitarch,  Inc. 
Albuquarqu  •,  Naw  M  •xico 


TECHNICAL  REPORT  NO.  AFWL-TR-70-20 


Juna  1970 

AIR  FORCE  WEAPONS  LABORATORY 
Air  Forca  Systems  Command 
Kirtland  Air  Fore#  Basa 
Now  Moxico 


P  D  c 


t  .  1 


23 


1970 


S9/1/IS 


% 


Thl.  doctaan t  1.  .ubj.ct  to  .paclal  axport  control.  and 
to  foralgn  gov.mn.nt.  or  for.lgn  natlon.1.  nay  ba  nada 
approval  of  AFVL  (WLCT)  .  Klrtland  AfB.  m,  87117, 


•a  eh 

only 


trarralttal 
with  prior 


DISCLAIMER 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


AFWL-TR-70-2 


v 


* 


ANALYTICAL  MODEL  FOR  HIGH  EXPLOSIVE  MUNITIONS  STORAGE 

H.  L.  Schreyer 
L.  E.  Romesberg 


Mechanics  Research,  Inc. 
Albuquerque,  New  Mexico 


TECHNICAL  REPORT  NO.  AFWL-TR-70-20 


This  document  is  subject  to  special  export  controls 
and  each  transmittal  to  foreign  governments  or 
foreign  nationals  may  be  made  only  with  prior 
approval  of  AFWL  (WLCT) ,  Kirtland  AFB,  NM  87117. 
Distribution  is  limited  because  of  the  technology 
discussed  in  the  report. 


AFWL-TR-70-20 


* 


I  ♦ 


> 


FOREWORD 


This  report  was  prepared  by  Mechanics  Research,  Inc.,  Albuquerque,  New 
Mexico,  under  Contract  F29601-69-C-0034.  The  research  was  performed  under 
Project  1597,  Task  12. 

Inclusive  dates  of  research  were  February  1969  through  March  1970.  The 
report  was  submitted  24  March  1970,  by  the  Air  Force  Weapons  Laboratory  Project 
Officer,  Captain  Jacob  C.  Armstrong  (WLCT). 

Information  in  this  report  is  embargoed  under  the  U.S.  Export  Control  Act 
of  1949,  admlnstered  by  the  Department  of  Commerce.  This  report  may  be 
released  by  departments  or  agencies  of  the  U.S.  Government  to  departsients  or 
agencies  of  foreign  governments  with  which  the  United  States  has  defense  treatv 
commitments,  subject  to  approval  of  AFWL  (WLCT). 

This  report  has  b -;en  reviewed  and  is  approved. 


CLARENCE  E.  TESKE 
Major,  USAF 

Chief,  Aerospace  Facilities  Branch 


Chief,  Civil  Engineering  Division 


ii 


ABSTRACT 


Analytical  models  and  subsequent  computer  codes  have  been 
developed  for  predicting  peak  overpressure,  positive  unit  im¬ 
pulse,  the  distribution  and  impact  velocity  of  bomb  fragments, 
crater  dimensions  and  ejecta  thickness  from  the  detonations  of 
typical  bomb  stacks  used  by  the  Air  Force.  These  models  consider 
aboveground  barricaded  stacks  with  an  equivalent  net  weight  high- 
explosive  range  of  10  to  500  tons  of  TNT.  The  peak  overpressure 
and  impulse  from  a  detonation  are  obtained  by  modifying  the 
known  results  of  a  bare  hemispherical  charge  to  take  into  account 
the  stack  and  barricade  geometries  and  the  interaction  effect 
of  bombs.  Fragment  dispersion  patterns  are  predicted  by  com¬ 
bining  experimental  results  for  single  bombs  and  using  the 
trajectory  equations  for  the  motion  of  a  steel  fragment  in  air. 

By  using  basic  principles  and  experimental  data,  crater  and 
ejecta  dimensions  are  predicted.  Based  on  output  from  the  com¬ 
puter  codes,  illustrative  examples  are  given  together  with 
recommendations  for  future  tests  to  obtain  needed  data.  Programs 
for  optimizing  munition  storage  areas  are  also  suggested. 

(Distribution  Limitation  Statement  No.  2) 
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SECTION  I 


INTRODUCTION 

One  of  the  major  problemo  at  many  Air  Force  installations 
is  the  storage  of  large  quantities  of  munitions.  Safety  consid¬ 
erations  for  protection  of  personnel  and  material  in  the  event 
of  accidental  detonation,  although  overriding,  are  in  direct  con¬ 
flict  with  economics,  i.e.,  large  clear  zones  require  considerable 
real  estate  with  resulting  long  roads  and  utility  lines.  If 
munitions  storage  clearance  requirements  could  be  reduced  with¬ 
out  endangering  safety  requirements,  significant  economic  gains 
could  be  realized. 

Full-scale  tests  of  munitions  storage  concepts  have  been 
conducted  and  have  yielded  valuable  information  leading  toward 
more  rational  munitions  storage  criteria.  Such  tests  are,  how¬ 
ever,  expensive  and  exceedingly  time  consuming.  Availability  of 
analytical  procedures  which  could  be  used  with  some  confidence 
to  predict  the  effects  of  detonation  of  a  stack  of  munitions  would 
be  invaluable  in  analyzing  new  storage  concepts,  or  in  rational, 
effective  planning  of  such  future  full  scale  tests  as  may  be  re¬ 
quired.  This  study  is  a  step  toward  development  of  such  proce¬ 
dures  . 

The  most  significant  parameters  in  determining  munitions 
detonation  hazards  include  peak  overpressure,  unit  impulse,  mass 
and  velocity  of  projectiles  formed  from  bomb  fragments  and  their 
distribution,  crater  dimensions,  and  the  probable  ejecta  distrib¬ 
ution.  This  report  outlines  an  analytical  model  that  adequately 


predicts  these  parameters  for  the  range  of  25,000  to  500,000 
pounds  net  weight  of  high  explosives.  Associated  with  the  re¬ 
port  are  computer  programs  that  perform  the  numerical  analysis 
required  for  the  models.  Even  though  the  most  recent  sources 
have  been  consulted,  it  was  considered  advisable  to  construct  the 
model  and  the  computer  programs  in  a  manner  that  would  easily 
allow  alterations  as  new  experimental  and  theoretical  work  became 
available. 

Section  II  considers  the  peak  overpressure  and  unit  impulse 
emanating  from  bomb  stacks,  both  barricaded  and  unbarricaded. 
Section  III  handles  fragmentation  while  Section  IV  discusses 
cratering.  The  corresponding  computer  programs  are  listed  and 
described  in  Appendices  I,  II  and  III  respectively.  Examples  of 
results  obtained  from  the  computer  programs  are  presented  in 
Section  V.  As  a  result  of  the  extensive  literature  survey  that 
was  conducted  and  from  the  formation  of  the  analytical  model*, 
it  became  apparent  that  further  investigations,  both  experi¬ 
mental  and  analy tical ,  in  certain  critical  areas  would  be 
highly  beneficial.  Such  a  program  is  outlined  in  Section  VI. 

Throughout  the  report,  an  attempt  has  been  made  to  use  a 
judicious  combination  of  basic  principles  and  results  from  small 
and  full-scale  tests.  Such  an  approach  is  considered  necessary 
if  the  results  are  to  be  used  for  typical  situations  that  current¬ 
ly  confront  the  Air  Force.  Furthermore,  the  use  of  fundamental 
concepts  implies  that  new  situations  can  be  handled  with  some 
degree  of  confidence.  However,  it  should  always  be  kept  in  mind 
that  soil  conditions,  for  exc^.  le,  can  change  with  time  and 


2 


accordingly,  even  well  designed  experiments  produce  data  with 
a  considerable  amount  of  scatter.  Accordingly,  a  certain  amount 

of  engineering  judgment  is  required  in  connection  with  the  results 
of  this  study. 
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SECTION  II 
BLAST  EFFECTS 

1.  INTRODUCTION 

The  objective  of  this  section  is  to  develop  en  analytical 
nodel  which  will  predict  the  environment  produced  by  the  air 
blast  from  a  high  explosive  detonation.  The  primary  parameter, 
to  be  investigated  are  the  pea*  overpressure  and  positive  im¬ 
pulse  experienced  at  all  points  on  the  surface  surrounding  a 
high-order  surface  detonation.  The  effect,  to  be  investigated 
include:  a)  the  effect  of  substituting  conventional  Air  Force 
bombs  for  TNT  in  the  explosive  stac*,  b,  the  effect  on  the  pea* 
overpressure-scaled  distance  P-Z  and  scaled  positive  impulse- ' 
scaled  distance  I-z  relationship,  produced  by  the  explosive 
stack  geometry,  and  c)  the  effect  on  the  P-z  and  I-z  relation¬ 
ships  produced  by  a  barricade  surrounding  the  explosive  stack  on 
three  sides  (standard  open-end  b-ricade) .  The  scaled  distance 

2  13  defi"ed  to  be  the  di3tan-  *  from  the  point  of  detonation 
divided  by  the  equivalent  charge  weight  in  pounds  of  TNT  to  the 
one -third  power  W1^3 

Z  =  R/W1//3 

(1) 

The  scaled  positive  impulse  I  is  defined  to  be  the  positive  unit 

impulse  is  divided  by  the  equivalent  charge  weight  in  pounds  of 

^NT  uO  the  one-third  power  W3-^3 

I  =  i  /W1/3 

s  m 
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The  general  approach  followed  ia  to  develop  a  model  to  pre¬ 
dict  the  P-Z  and  I-Z  relationships  for  a  surface  detonation  of 
a  hemlepli  srlcal  stack  of  high  explosive  TNT,  and  then  to  modify 
these  relationships  to  account  for  the  Individual  effects  listed 
above.  The  effect  produced  by  changing  the  point  of  detonation 
in  the  stack  is  assumed  to  be  negligible.  (See  Reference  1). 


2.  BARE  CHARGE  PARAMETERS 

The  initial  task  in  the  development  is  to  model  the  environ¬ 
ment  produced  by  a  surface  detonation  of  a  bare,  i.e.,  unbarri¬ 
caded,  hemispherical  stack  of  '-’NT,  with  respect  to  peak  over¬ 
pressure  and  scaled  positive  impulse.  Curves  describing  the 

« 

pe*k  overpressure-scaled  distance  relationship  for  high  explo¬ 
sive  surface  detonations  are  available  throughout  the  litera¬ 
ture.  These  relationships  have  been  developed  through  many 
years  of  full  scale  testing  and  are  widely  accepted.  Discrep¬ 
ancies  do  appear  in  the  literature  when  comparing  the  relation¬ 
ships  published  by  one  testing  agency  with  those  published  by 
another;  however,  these  discrepancies  are  of  a  relatively  small 
order.  The  relationships  selected  for  the  model  development 
(See  Figure  1)  are  published  in  Reference  2. 

To  carry  out  the  objectives  of  this  section,  it  is  necessary 
to  have  these  results  available  in  a  numerical  form.  The  pro¬ 
cedure  used  in  modeling  the  P-Z  and  I-Z  relationships  is  as 
follows:  a)  Points  on  the  P-Z  and  I-Z  curves  are  selected;  b) 
the  coordinates  (P,  Z)  and  (I,  Z)  of  these  points  a;  e  t  ransformed 
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by  computing  the  natural  logarithm  of  each  coordinate  so  that 
the  coordinates  (In  P,  In  Z)  and  (In  I,  In  Z)  are  obtained; 

c)  the  coordinates  (In  P,  In  Z)  and  (In  I ,  In  Z)  are  used  in  a 

least  squares  polynomial  curve  fit  program  to  obtain  relation¬ 
ships  of  the  form 

In  P  «  aQ  +  a^  (In  Z)+  (In  Z)^  +  ....  and  (3) 

In  I  *  b  +  b.  (In  Z)  +  b~(ln  Z)2  +  _  (4) 

O  1  2 

d)  the  polynomial  coefficients  obtained  from  the  curve  fit  pro¬ 
gram  are  then  used  to  evaluate  the  value  of  peak  overpressure  P 
and  scaled  positive  impulse  I  at  the  desired  values  of  scaled 
distance  Z. 

This  procedure  yields  results  which  are  in  very  close 
agreement  with  the  original  relationships  (curves) .  The  maxi¬ 
mum  error  in  the  predicted  peak  overpressure-scaled  distance 
relationship  is  less  than  7%  for  .5<Z<10  ft/lb^3  (P>10  psi)  , 
less  than  3%  for  10<Z<45  ft/lb^3  (1.0>P>10psi)  and  less  than 
5.0%  for  45<Z<500  ft/lb1^3 (P<1. 0  psi).  The  maximum  error  in 
the  predicted  scaled  positive  impulse-scaled  distance  relation¬ 
ship  is  less  than  6%  for  .5<Z<10  ft/lb^3 (I>10.0  psi-ms/lb^3) , 
less  than  5%  for  10<Z<75  ft/lb 1/3  (10. 0>I>1.0  psi-ms/lb1/3)  and 
less  than  15%  for  75<Z<500  ft/lb 1/3 (I<1. 0  psi-ms/lb1/3) . 

The  resulting  coefficients  for  determining  the  overpressure 
and  scaled  positive  impulse  according  to  Equations  (3)  and  (4) 
are  given  in  Table  I. 
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TABLE  I 


POLYNOMIAL  COEFFICIENTS  FOR  DETERMINING  PEAK  OVERPRESSURE 
AND  SCALED  POSITIVE  IMPULSE  FOR  A  BARE  HEMISPHERICAL  CHARGE 

aQ  -  10.7036810  x  101 
ax  -  -0.1663724  x  101 
a2  =  -0.2516481  x  10° 
a3  =  -0.1137714  x  10° 
a4  =  +0.3818405  x  10_1 
a5  =  +0.5035198  x  10-1 
afi  =  -0.2756970  x  10"1 
a?  -  +0.5557968  x  10"2 
a8  =  -0.5108014  x  10“3 
a9  -  +3.1795565  x  10"4 


bQ  -  +0.3129288  x  101 
b1  -  -0.1295979  x  10° 
b2  »  +0.4112452  x  10° 
b3  =  -0.7687394  x  10° 
b4  =  +0.4969224  x  10° 
b5  =  -0.1684197  x  10° 
bg  =  +0.2805656  x  10" 
b?  =  -0.1791292  x  10" 
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3 .  BOMB  EFFECT 


The  effect  produced  by  substituting  conventional  Air  Force 
bombs  for  TNT  in  the  explosive  stack  has  been  studied  by  several 
authors  for  many  years.  This  effect  is  accounted  for  in  the 
model  development  through  the  use  of  a  bomb  factor.  This  factor 
is  multiplied  by  the  total  explosive  weight  in  the  bomb  stack 
to  yield  an  equivalent  weight  of  TNT.  The  bomb  factor  includes 
the  confined  explosion  effect,  the  surface  reflectivity  effect, 
and  the  individual  bomb  geometry  effect.  Typical  bomb  factors 
can  be  found  in  Reference  3.  No  attempt  has  been  made  to  account 
for  the  explosion  confinement  effect  caused  by  stacking  bombs. 
There  will  be  some  confinement  effect  caused  by  surrounding  a 
bomb  by  other  bombs.  However,  since  there  apparently  is  no 
empirical  or  theoretical  data  available,  this  effect  cannot  be 
accounted  for  in  the  present  model. 


4.  STACK  GEOMETRY  EFFECT 

It  is  a  well  accepted  fact  that  the  geometry  of  the  explo¬ 
sive  stack  has  a  great  effect  on  the  peak  overpressure  ±md 
positive  impulse  at  positions  "close"  to  the  stack.  This 
effect  diminishes  with  distance  from  the  stack.  Even  though 
this  fact  is  well  known  and  accepted ,  there  is  very  little 
data  available  in  the  current  literature  which  quantitatively 
describes  such  a  variation. 

Apparently  Reference  4  describes  one  of  the  few  attempts 
to  measure  the  effect  of  the  geometry  of  charges  on  peak 
overpressure  and  positive  impulse.  The  report  is  composed 
basically  of  peak  overpressure  and  positive  impulse  measure¬ 
ments  using  eight  charge  shapes  composed  of  50  pounds  of  RDX 
composition  C-3  explosive  which  is  equivalent  to  54.5  pounds 
of  TNT. 

It  will  be  assumed  that  a"standard"  high  explosive  bomb 
stack  can  be  approximated  by  a  solid  stack  of  explosives  with 
a  rectangular  solid  shape.  Therefore,  the  only  shapes  consid¬ 
ered  by  the  above  mentioned  report  which  are  pertinent  to  the 
model  development  are  the  cubical  shaped  charge  and  the  plate 
shaped  charge.  The  length  of  a  side  of  the  cubical  shaped  charge 
used  in  the  test  was  9.6  in.  x  9.6  in.  while  the  plate  dimensions 
were  54.1  in.  x  9.0  in.  x  1.8  in.  Pressure  and  impulse  measure¬ 
ments  were  made  at  35  ft.,  45  ft.,  60  ft.,  70  ft.,  and  80  ft., 
from  the  center  of  the  charges,  along  lines  perpendicular 
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to  each  face  and  through  the  charge  center  for  the  cube  and  the 
plate,  and  along  lines  thrcugh  the  edge  and  the  charge  center 
in  a  horizontal  plane  for  the  cube.  These  distances  from  the 
charge  center  correspond  to  scaled  distances  of  9.22,  11.85, 
15.80,  18.43,  and  21.06  ft/lb1^,  respectively.  The  average 
peak  overpressures  and  positive  impulses  measured  at  the  above 
locations  for  the  sphere,  cube  and  plate  are  shown  in  Tables 
II  and  III,  respectively  (See  Tables  la,  lb,  2a,  2b,  7a,  and 
7b  of  Reference  4)  . 

The  author  of  Reference  4  commented  that  the  accuracy  of 
the  data  did  not  warrant  an  attempt  at  curve  fitting.  Since 
the  development  of  the  geometry  effect  portion  of  the  model  is 
based  solely  on  the  data  from  this  one  report,  steps  were  taken 
to  smooth  out  some  of  the  inaccuracies  of  the  test  data.  To 
reduce  the  effect  on  readings  by  individual  differences  in 
recording  instruments,  system  circuitry,  and  drift  from  zero 
calibration  point  before  testing,  the  ratio  of  measurements 
for  shaped  charges  divided  by  measurements  from  spherical 
charge  rather  than  actual  measurements  was  used.  In  other  words, 
the  model  development  deals  with  the  effect  of  going  from  a 
spherical  charge  to  a  rectangular  charge  (i.e.  ,  cube  and  plate) 
rather  than  dealing  with  the  rectangular  charges  at  face  value. 

A  parametric  study  performed  on  the  data  revealed  that  a 
reasonable  approach  to  the  model  development  would  be  to  analyze 
the  data  with  respect  to  an  "area  ratio"  scheme.  The  reasons 
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for  adopting  this  scheme  for  a  rectangular  charge  are  as  follows 

a.  The  total  energy  released  to  the  surrounding  atmos¬ 
phere  by  a  charge  resting  on  a  surface  should  be 
proportional  to  the  total  surface  area  exposed  to  the 
atmosphere  (this  excludes  the  face  of  the  charge  which 
is  in  contact  with  the  surface)  , 

b.  The  energy  experienced  by  a  point,  which  lies  on  a 
line  perpendicular  to  one  of  the  exposed  faces  of  the 
charge  and  through  the  charge  center,  should  be  pro¬ 
portional  to  the  total  energy  released  by  the  detona¬ 
tion  of  the  charge  times  the  ratio  of  the  charge  face 
area  nearest  to  the  point  to  the  total  charge  surface 
area  exposed  to  the  atmosphere.  This  ratio  is  here¬ 
after  referred  to  as  the  "face  area  ratio"  (FAR),  and 
is  computed  as  follows  (See  Figure  2)  : 


PARP1 


_ SH  X  SD _ 

2  (SH  x  SD)  +  2(SH  x  SL)  +  (SL  x  SD) 


(5a) 


_ SH  X  SL _ 

AKP2  2  (SH  x  SD)  +  J(SH  x  SL)  +  (SL  x  SD) 


(5b) 
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The  values  of  FAR, which  correspond  to  the  charge  orienta¬ 
tions  and  measurement  directions  for  the  cube  and  plate  data 
of  Tables  II  and  III,  are  FAR2  -  0.20,  FAR^  •  0.0147,  FARg  - 
0.440,  FARg  =  0.0227,  and  FAR^  ■  0.1365,  where  the  subscripts 
correspond  to  the  numbered  columns  of  both  Tables.  The 
magnitude  of  FAR  will  always  be  greater  than  0.00  and  less 
than  0.50.  The  data  in  columns  4  and  6  of  Tables  II  and  III 
are  ignored  in  the  analysis  because  it  is  felt  that  the  values 
of  FAR,  for  these  conditions,  are  much  lower  than  will  ever  be 
experienced  in  an  actual  bomb  stack.  Figure  3  shows  the  pro¬ 
duct  of  overpressure  ratios  (face  overpressure  from  rectangular 
stacks  divided  by  overpressure  from  spherical  stack)  and  FAR 
plotted  against  distance  from  stack  center.  Figure  4  shows 
similar  curves  for  impulse  ratios.  Note  that  for  large  dis¬ 
tances  from  charge  center  these  curves  approach  the  value  of 
FAR  in  each  case.  These  curves  not  only  vary  with  the  value 
of  FAR,  but  also  with  the  distance  from  the  charge  center. 

At  each  value  of  R,  or  the  corresponding  value  of  scaled 
distance  Z,  however,  a  relationship  between  pressure  ratio  or 
impulse  ratio  and  FAR  can  be  established  for  the  three  curves 
shown  in  each  figure.  These  relationships  can  be  expressed  in 
the  form 

_ l _  -  B  +  C  (FAR)  +  D  (FAR)2  (6a) 

Psphere 

j -  *  F  +  G  (FAR)  +  H  (FAR)2  (6b) 

sphere 
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Figure  3.  Overpressure  Ratio  x  FAR  vc.  Distance  from 
Charge  Center 
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FAR 


where  B ,  c,  D,  F,  G,  and  H  are  coefficient*  which  can  be  solved 
exactly  for  any  fixed  value  of  Z  since  in  each  case  there  are 
three  unknowns  and  three  equations. 

In  these  relations  data  from  charge  shapes  with  different 
equivalent  amounts  of  TNT  can  be  used  since  the  scaled  distance 
Z  can  be  used. 

The  solutions  to  these  equations  for  the  values  of  R  or  Z 
given  in  Tables  II  and  III  are  shown  in  Table  IV  which  indicates 
that  B ,  C ,  D,  F,  G ,  and  H  vary  with  the  value  of  Z.  It  was 
assumed  that  fifth  order  polynomials  would  adequately  describe 
this  variation  in  Z.  A  least  squares  polynomial  curve  fit 
program  was  used  to  fit  the  available  data  from  Reference  4. 

Also,  the  geometry  effect  must  vanish  for  large  Z.  Hence  there 
were  six  data  points  available  for  the  curve  fitting  routine 
(Table  IV)  . 

Denote  the  general  form  of  these  polynomials  by 

Br  C,  D,  F,  G,  H  ■  C__  +  C,Z  +  C-Z2+  C-Z3+  C.Z4+  CcZ5  (7) 

O  X  /  j  4  2 

The  computed  values  of  the  coefficients  are  given  in 
Table  V. 

Figure  5  gives  the  overpressure  and  impulse  ratios  based 
on  Equations  (6)  and  (7)  for  points  along  a  line  perpendicular 
to  the  center  of  one  of  the  vertical  face*  of  a  cubical  shaped 
explosive  stack,  that  is,  FAR  -  0.2.  Note  that  the  predicted 
curves  fit  the  data  of  Reference  4  quite  accurately 
as  well  as  approaching  the  value  1  for  large  Z.  This  is  to  be 
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POLYNOMIAL  COEFFICIENTS  FOR  B 
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5.  Comparison  of  Predicted  Overpressure  and  Impulse 
Ratios  to  BRL  Data  for  Cubical  Shaped  Explosive 
Stack. 
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expected  since  these  points  were  used  in  determining  the 
coefficients.  A  similar  conclusion  holds  when  FAR  *  0.44  and 
FAR  »  0.1365.  The  polynomial  expressions  are  necessary  for 
determining  the  overpressure  and  impulse  ratios  for  values  of 
Z  and  FAR  not  given  explicitly  by  Reference  4. 

Reasonable  results  for  extrapolations  to  values  of  Z  in 
the  region  Z>9  can  be  expected  because  of  the  experimental  data 
that  is  available.  However,  for  values  of  Z  lees  than  9  the 
predicted  values  may  not  be  very  accurate. 

To  obtain  approximate  values  for  the  actual  peak  over¬ 
pressure  and  positive  impulse  for  positions  along  a  line  per¬ 
pendicular  to  the  center  of  a  vertical  stack  fae*,  the  over¬ 
pressure  and  positive  impulse  ratios  are  multiplied  by  the 
peak  overpressure  and  positive  impulse  respectively.  These  are 
obtained  from  the  polynomial  fit  for  curves  associated  with 
the  surface  detonation  of  a  hemispherical  stack  at  the  same 
values  of  scaled  distance  Z.  Although  the  overpressure  and 
impulse  ratios  involve  spherical  charges,  the  results  from  a 
hemispherical  detonation  are  used  to  convert  these  ratios  to 
true  peak  overpressures  and  positive  impulses.  The  basic 
reason  for  using  data  associated  with  a  hemispherical  charge 
is  that  the  positive  impulse  is  larger  than  that  associated 
with  a  spherical  charge  for  the  same  value  of  Z.  Presumably 
this  would  yield  conservative  values  of  positive  impulse,  which 
is  especially  needed  in  the  region  Z<9  where  no  experimental 
data  is  available. 
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It  is  assumed  that  there  is  very  little  difference  in 
peak  overpressure  and'  positive  impulse  between  the  center  line 
mentioned  above  and  the  corresponding  line  on  the  surface  ver¬ 
tically  below  the  center  line.  Such  an  assumption  will  be 
implicitly  assumed  from  now  on  for  all  other  horizontal  lines 
emanating  from  the  center  of  a  rectangular  stack. 

Now  that  the  overpressure  and  impulse  at  positions  out 
from  the  center  of  a  vertical  face  of  n  rectan" -iar  stack  have 
been  modeled,  the  next  step  is  to  predict  these  parameters  out 
from  the  vertical  edges.  The  only  data  available  on  which  this 
development  can  be  based  are  the  data  shown  in  colunn  3  in 
Table  II  and  Table  III. 

The  technique  employed  in  the  development  of  the  edge 
peak  overpressure  and  scaled  impulse  versus  scaled  distance 
was  to  establish  edge  to  face  peak  overpressure  and  edge  to 
face  positive  impulse  relationships  as  functions  of  scaled 
distance.  This  was  achieved  by  dividing  the  values  in  column  3 
of  Tables  II  and  III  by  the  values  in  column  2  of  the  corres¬ 
ponding  tables,  thereby  establishing  coordinates  for  five 
points  for  the  ratio  relationships.  These  coordinates  were  put 
into  a  least  squares  polynomial  curve  fit  program  to  obtain  the 
polynomial  coefficients  for  the  polynomials  in  Z  to  describe  the 
ratio  relationships.  Table  VI  shows  the  polynomial  coefficients 
for  the  relationships  that  yield  the  ratio  of  edge  peak  over¬ 
pressure  to  face  peak  overpressure  and  edge  impulse  to  face 
impulse.  Figure  6  shows  these  relationships  as  functions  of 
scaled  distance. 
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TABLE  VI 


POLYNOMIAL  COEFFICIENTS  FOR  THE  RATIOS 
EDGE  PRESSURE/FACE  PRESSURE  AND 
EDGE  IMPULSE/FACE  IMPULSE 

6 


EDGE  PRESSURE 


dQ  -  -0.5442047  x  10_1 
dx  -  -0.3279577  x  10~2 
d2  -  0.3172064  x  10_1 

d3  -  -0.2700095  x  10~2 

d4  -  0.8668014  x  10~4 

d5  -  -0.1294399  x  10“5 

dfi  -  0.7019624  x  lo"8 

4 

EDGE  IMPULSE  J  i 

FACE  IMPULSE  ",  2 

i-0 


eQ  -  -0.3879756  x  10° 
e2  -  0.2953776  x  10° 
e2  “  “0.1519496  x  10"1 
e3  -  0.2943603  x  10_3 
e4  -  -0.1946069  x  lo“5 
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Edge  Pressure/Face  Pressure  and  Edge  Impulse/Face  Impulse  vs.. Scaled 
Distance  from  Charge  Center. 


These  relationships  were  established  from  data  for  a 
cubical  shaped  charge  where  the  pressures  and  impulses  off  of 
the  faces  as  functions  of  Z  are  equal.  The  technique  employed 
in  the  model  for  predicting  edge  overpressures  and  impulses  for 
stacks  which  are  not  cubical  (i.e.,  face  pressure  and  inpulse  as 
functions  of  Z  are  not  equal  for  adjoining  faces)  is  to  multiply 
the  ratios,  at  a  given  value  of  Z,  by  the  average  face  value  of 
peak  overpressure  and  impulse  at  the  same  value  of  Z. 

Now  that  the  peak  overpressure  and  positive  impulse  rela¬ 
tionships  as  functions  of  Z  have  been  established  along  lines 
perpendicu&lar  to  the  stack  faces  through  the  stack  center,  and 
along  lines  extending  from  the  stack  center  through  the  stack 
edges  (all  lines  in  a  horizontal  plane)  the  parameter  values 
along  intermediate  lines  through  the  stack  center  can  be  esta¬ 
blished  by  linear  interpolation.  This  technique  is  illustrated 
in  Figure  7. 


5 .  BARRICADE  EFFECT 

It  is  a  well  known  and  accepted  fact  that  a  barricade  in 
close  proximity  to  an  explosive  detonation  will  significantly 
affect  the  peak  overpressure  and  positive  impulses  at  positions 
’’close"  to  the  barricade.  Although  there  has  been  considerable 
study  dealing  with  qualitative  (i.e.,  amount  of  destruction) 
effects  produced  by  barricaded  explosive  charges,  there  has  not 
been  much  study  concerning  quantitative  (i.e.,  actual  pressure 
and  inpulse  measurement)  effects  produced  by  barricaded  explo¬ 
sive  charges  in  the  current  literature. 
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Figure  7.  Geometry  Effect  and  Interpolation  Technique 
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The  only  reference  revealed  by  the  literature  search  which 
deals  directly  with  the  effects  produced  by  a  standard  (i.e., 
three  adjoining  walls  perpendicular  to  each  other)  barricade  is 
Reference  2.  This  report  qives  incident  pressures  as  a  function 
of  scaled  distance  (See  Figure  8)  for  various  directions  of 
propagation  from  a  three  sided  barricade.  The  barricade  length 
to  depth  ratio  is  approximately  one  (1)  and  the  weight  of  charge 

3 

to  voir..  ..e  of  structure  ratio  (W/V)  (pounds  of  TNT/ft.  )  is  in 
the  range  0.2  to  2.0.  For  very  large  or  small  W/V  values,  the 
incident  pressure  versus  scaled  distance  in  all  directions  of 
propagation  from  a  barricade  will  be  very  nearly  equal  to  the 
results  for  an  unbarricaded  charge.  Results  of  barricade  effects 
on  positive  impulse  versus  scaled  distance  for  various  direc¬ 
tions  of  propagation  from  the  barricade  considered  are  not 
reported  here  or  elsewhere  in  the  literature. 

The  technique  employed  in  the  development  of  this  portion 
of  the  model  was  to  establish  peak  overpressure  ratio  relation¬ 
ships  for  the  four  directions  from  the  barricade  center  shown  in 
Figure  8  as  functions  of  scaled  distance.  This  was  accomplished 
by  dividing  the  pressure  values  from  the  curve  for  an  unconfined 
surface  burst  shown  in  Figure  8  by  the  pressure  values  from 
the  curves  for  the  four  directions  from  the  barricade  center  at 
selected  values  of  scaled  distance.  These  point  coordinates 
were  then  put  into  a  least  squares  polynomial  curve  fit  program 
to  obtain  coefficients  for  fourth  order  polynomials  in  terms 
of  scaled  distance  Z  which  describe  the  overpressure  ratios  for 


29 


.1 


1000 


I 


Figure  8 


10  100 

SCALED  DISTANCE,  Z,  (FT/l.BS^) 

Exterior  Leakage  Pressure  vs  Scaled  Distance 


30 


each  of  the  four  directions  from  the  barricade  center  shown  in 
Figure  8.  The  polynomial  coefficients  are  shown  in  Table  VII. 

Pressure  ratio  relationships  for  directions  between  those 
shown  are  established  through  linear  interpolation  as  was  done 
in  the  geometry  effects  model  development. 

Since  the  barricade  effect  on  positive  impulse  cannot  be 
established  because  of  the  lack  of  reliable  data,  the  model 
assumes  that  positive  impulse  is  affected  in  the  same  way  as 
is  the  pressure.  Therefore,  the  above  developed  overpressure 
ratio  relationships  are  reused  as  positive  impulse  ratio 
relationships. 

The  technique  to  establish  the  peak  overpressure  and  posi¬ 
tive  impulse  produced  by  a  barricaded  explosive  detonation  is 
to  evaluate  the  above  ratio  relationships  at  the  desired  direc¬ 
tion  and  scaled  distance  from  the  barricade  and  multiply  the 
computed  ratio  value  by  the  overpressure  and  impulse  of  the 
bare  charge  at  the  same  value  of  scaled  distance.  If  geometry 
effects  are  included  in  the  problem,  the  computed  ratio  value 
is  multiplied  by  the  peak  overpressure  and  positive  impulse 
which  has  been  previously  modified  to  account  for  the  stack 
geometry  effect. 

It  should  be  emphasized  that  this  approach  assumes  that 

the  explosive  stack  and  barricade  are  rectangular  with  stack 

and  barricade  sides  parallel  to  each  other,  that  tb**  same  vertical 

line  passes  through  the  center  of  the  stack  and  the  barricade,  that 

the  ratio  of  barricade  length  to  barricade  depth  is  approximately 

equal  to  o.ie,  and  that  the  ratio  of  the  weight  of  charge  to  the 

3 

volume  of  barricade  (lbs/ft  )  is  in  the  range  of  0.2  to  2.0. 


6 .  SUMMARY 


This  section  has  outlined  the  theory,  experimental  data 
and  assumptions  that  have  been  utilized  in  developing  a 
workable  model  for  predicting  peak  overpressures  and  positive 
impulses  associated  with  a  barricaded  HE  detonation.  The 
basic  approach  was  to  use  well  known  results  from  detonations 
of  bare  hemispherical  TNT  charges  and  to  modify  these  results 
to  account  for  the  effects  of  bombs,  rectangular  stack  con¬ 
figurations  and  barricades.  For  those  areas  where  a  certain 
amount  of  uncertainty  prevailed,  a  conservative  approach  was 
adopted  so  that  the  predicted  values  would  be  larger  than 
what  might  be  actually  experienced  in  practice. 

The  computer  program  which  follows  the  theory  outlined  in 
this  section  is  described  in  Appendix  I. 


33 


SECTION  III 

fragmentation 

1  •  INTRODUCTION 

The  purpose  of  the  fragmentation  portion  of  the  program 

formulate  an  analytical  model  to  describe  the  fragment 

aversion  pattern  resulting  from  the  explosion  of  barricaded 

The  model  muat  consider  the  di,persions  an(J  patterns 

^  ter,”S  °£  fra?ment  vel-itiea,  weights,  trajectories  and 
ranges . 


2.  STATEMENT  OF  PROBLEM 

The  fragment  dispersion  pattern  of  a  barricaded  explosion 

18  aff8Cted  Primarily  *  the  P«»eters:  a,  initial  fra -men  t 
velocity,  b,  fragment  mass,  c,  spatial  position  on  the  bomb, 

d)  lnltlal  departure  an9le'  e>  fra9ment  trajectories,  f,  bomb 
staot  geometry,  g,  and  barricade  geometry.  The  analytical 

-del  developed  with  these  parameter,  must  yield  fragment  die- 

patterns  that  are  at  least  comparable  to  dispersion 

data  from  barricaded  explosive  fOBfa 

explosive  tests.  The  model  must  be  able 

to  yield  the  range  and  the  striking  velocity  of  the  fragments 

order  to  determine  the  danger  of  sympathetic  detonation  of 
adjacent  bomb  stacks. 
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3.  SCOPE  OF  INVESTIGATION 

The  basic  approach  for  predicting  the  above  parameters  for 
the  explosion  of  a  barricaded  bomb  stack  is  to  proceed  from 
known  fragment  behavior  for  a  single  bomb.  Representative 
single  bomb  fragment  data  that  are  available  (References  5-7) 
for  a  large  variety  of  bombs  is  shown  in  Figure  9.  These  data 
consist  of:  a)  the  number  of  fragments  per  steradian  (solid 
angle),  b)  the  average  fragment  mass  (grams)  ,  as  a  function  of 
the  polar  angle  measured  from  the  nose  of  the  bomb. 

Once  the  initial  fragment  velocity  and  fragment  weight  are 
known,  the  trajectory  range  and  striking  velocity  may  be  pre¬ 
dicted  for  various  angles  of  departure.  These  data  are  first 
generated  for  the  explosion  of  a  barricaded  single  bomb  and 
is  then  correlated  with  fragment  survey  data  from  experimental 
bomb  stack  tests  such  as  the  "BIG  PAPA"  tests  (Refers nee  1).  The  corre¬ 
lation  with  experimental  fragment  survey  data  yields  interaction 
coefficients  which  indicate  the  number  of  single  bombs  required  to  yield 
the  fragment  data  from  the  stack  tests. 


35 


Figure  9.  Representative  Fragmentation  Data 
For  The  Explosion  of  a  Single  Bomb. 
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4.  ANALYTICAL  FRAGMENT  MODEL  FOR  A  SINGLE  BOMB 
a.  Fragment  Parameters 

To  develop  a  model  for  a  bomb  stack,  it  is  necessary 
to  first  consider  the  dispersion  pattern  for  a  single  bomb. 
The  parameters  that  are  required  include  the  mass  and 
initial  velocity  of  the  fragments  and  the  number  of  frag¬ 
ments  per  steradian  that  are  emitted  from  the  bomb. 

Gurney's  Theory  and  Mott  and  Shapiro's  Theory  (Reference 
7)  are  available  for  predicting  the  initial  fragment  velo¬ 
city  and  mass  distribution  respectively.  Both  of  these 
theories  have  been  favorably  correlated  with  experimental 
data.  A  significant  disadvantage  of  these  theories  is  that 
there  is  no  method  for  predicting  how  the  mass  of  the  frag¬ 
ments  vary  along  the  length  of  the  bomb.  Hence,  a  distrib¬ 
ution  would  have  to  be  assumed  or  the  average  fragment  mass 
determined  from  Mott's  equation  could  be  used. 

Since  experimental  fragment  data  for  several  bombs  are 
available  in  the  literature,  as  shown  in  Figure  9,  it  was 
considered  appropriate  to  use  this  information  as  needed 
in  the  model  instead  of  using  an  exclusively  theoretical 
approach . 

By  using  the  data  from  fragment  tests  on  single  bombs, 
the  following  parameters  are  known  as  functions  of  the 
polar  angle  y  (Figure  10):  a)  fragment  mass,  b)  fragment 
initial  velocity,  c)  fragment  mass  distribution  along  the 
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for  a  Bomb  Fragment. 
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21 


bomb,  d)  and  number  of  fragments  per  steradian  or  solid 
angle.  It  is  assumed  that  each  fragment  initially  departs 
along  the  line  through  the  bomb  center  and  the  point  on 
the  bomb  casing  at  which  the  fragment  is  located.  The 
angle  that  this  line  makes  with  the  horizontal  plane  is 
denoted  as  the  departure  angle  0. 
b.  Fragment  Trajectories 

The  trajectory  equations  are  necessary  to:  a)  determine 
if  fragments  clear  the  barricade,  b)  predict  the  fragment 
range  and  impact  velocity,  c)  predict  the  number  of  frag¬ 
ments/unit  area  and  mass/unit  area  as  a  function  of  azimuth 
angle  and  range  from  blast.  It  will  be  assumed  in  the  model 
that  if  fragments  do  not  clear  the  barricade,  they  are 
stopped  and  no  longer  considered. 

The  range  and  impact  velocity  will  be  predicted  by  using 
basic  equations  of  mechanics  for  the  trajectory  in  a  finite 
difference  form. 

The  trajectory  of  a  fragment  is  shown  in  Figure  11. 

equations  below  are  taken  from  References  8  and  9  and 
modified  for  purposes  of  this;  study. 
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The  forces  that  act  on  a  body  in  flight  are  the  drag 
PD  and  the  force  of  gravity  F^.  Assume  that  the  mag¬ 
nitude  of  the  drag  force,  acting  in  the  direction 
opposite  to  the  velocity,  is  given  by 

Fd  -  1/2  paACdV2  (8) 

where 

3 

PA  ■  density  of  air,  slugs/ft. 

A  «  cross-sectional  area  perpendicular  to  the 

2 

direction  of  propagation,  ft. 

Cp  «  drag  coefficient  which  is  a  function  of  shape 
V  ■  speed  of  fragment,  ft. /sec. 


The  components  of  the  drag  force  in  the  X  and  Y  directions 
are  given  respectively  by 


FdX  -  -PD  cosa  (9a) 

FdY  -  -PD  sina  (9b) 


where  a,  as  shown  in  Figure  11,  is  the  angle  between  the 
horizontal  line  and  the  tangent  to  the  trajectory  and 
varies  such  that 


B  >  a  >  m~aJ 
—  —  impact 


(10) 
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The  magnitude  of  the  force  due  to  gravity  is  given  by 


Fg  -  mg  (11) 

where 

m  «  mass  of  fragment/  slugs 

2 

g  »  gravitational  constant,  ft. /sec. 

Newton's  law  yields  the  following  equations  of  motion: 


„  X 

mx  -  fd 

(12a) 

..  ‘  Y 
my  -  Fd  -  Fg 

(12b) 

where 


x  ■  acceleration  in  X-direction 
y  ■  acceleration  in  Y-direction 


If  the  ballistic  coefficient  c  is  defined  by 


»aa  cd 


then  Equations  (12a)  and  (12b)  become 

x  ■  -cV2  coso 
y  ■  -cV2  sina  -  g 


(13) 

(14a) 

(14b) 


(21) 


By  combining  Equations  (16a)  and  (20) ,  we  get 

3 

d(v  cosot)  -  da 
9 


Multiply  Equation  (20)  on  the  left  hand  side  by 

(jE1)  and  on  the  right  by  its  equivalent  (V  cosa)  from 
Equation  (15a) 


dx  - 


(22) 


Similarly  from  Equations  (15b)  and  (20) 

_  -V2  tana  * 

*  - g -  da  (23) 

In  summary,  the  governing  equations  for  the  time, 
velocity  and  coordinates  of  the  fragment,  with  a  chosen 
to  be  the  independent  variable,  are 


dt 


g  cosa 


da 


<*(V  cosa)  -  da 
9 


dx  -  21  da 

A  “V2 

dy  ■  tana  da 


(24a) 

(24b) 

(24c) 

(24d) 
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The  initial  conditions  are 


For  a  numerical  solution  of  V,  x  and  y  in  terms  of 
a,  Equations  (24b)  ,  (24c)  and  (24d)  are  written  in  a 
finite  difference  form  as  follows: 


4V  -  V(t«na  ♦  SC—)  Ac 

Ax  *  Aa 

9 

-V2 

Ay  -  -  tana  Aa 


(26a) 

(26b) 

(26c) 


After  each  increment  in  a ,  the  coordinates  of  the 
fragment  can  to  determined  and  the  computation  stopped 
cnce  the  particle  strikes  the  barricade  or  the  ground. 

c.  Fragment  Ballistic  Coefficients 

In  order  to  integrate  the  trajectory  equations  of 
the  previous  section,  the  ballistic  coefficient  c  must 


I 


be  determined  for  each  fragment.  It  is  more  convenient 
to  express  Equation  (13)  in  the  form 


w 


F 


(27) 


where 


Ya  «  weight  density  of  air,  lb. /ft. 3 
wF  ■  weight  of  fragment,  lb. 

The  weight  density  of  air  is  assumed  known.  CD  is 
expected  to  take  on  values  ranging  from  0.3  to  2.0  and 
will  be  assumed  constant  for  any  fragment.  Reference  10 
indicates  a  value  of  0.6  for  the  drag  coefficient,  CD, 
as  most  appropriate  for  random  steel  fragments  from  an 
exploded  bomb  case.  The  remaining  term  that  is  needed 
is  the  ratio  AAp. 

Since  there  is  some  disagreement  as  to  appropriate 
values  for  this  ratio,  consider,  for  purposes  of  illustration 
only,  a  steel  cube  where  the  length  of  any  edge  is  i. 

Then  the  cross-sectional  area  A  must  lie  in  the  region 
l2  <  A  <  3i2cos  (54#44 ' ) 

If  the  density  is  taken  to  be  500  lb. /ft.3,  then 

w_  -  500i3 
F 

and 


The  range  for  the  area  to  weight  ratio  ie  then  given  by 


.0159  .  A  „  .0275 

-rrn  *7  <- 

P  F 


For  irregular  fragments  such  as  those  formed  by  the 
rupture  of  a  bomb  case/  the  area  to  weight  ratio  would 
be  expected  to  be  larger.  An  average  value  of 


.0345 


(28a) 


or 


.232 


(28b) 


has  been  suggested  by  Reference  11  where  wpG  is  the 
weight  of  the  fragment  in  grams.  This  average  value  has 
been  used  in  connection  with  the  numerical  analysis  for 
the  trajectory  portion  of  the  program  in  which  case  the 
fragment  ballistic  coefficient  becomes 


c  -  .00532/w, 


1/3 


(29) 


d.  Coordinate  System 

Most  typical  bomb  stack  barricades  are  rectangular 
with  three  closed  sides  and  one  open  as  shown  in 

Figure  12.  For  convenience,  we  adopt  the  coordinate  system 
shown  in  Figure  13  where  0  ■  0  corresponds  to  the  line 
coming  from  the  center  of  the  bomb  out  the  open  side  of 
the  barricade.  As  before,  6  is  the  angle  that  a  given 
line  makes  with  the  horizontal  plane.  The  bomb  is 
assumed  to  be  horizontal  and  0Q  denotes  the  orientation 
of  the  polar  axis  of  the  bomb. 

The  number  of  fragments  per  steradian  that  are 
emitted  from  a  bomb  (Figure  9)  is  given  as  a  function 
of  the  polar  angle  y.  Since  the  analysis  will  be  per¬ 
formed  using  the  coordinates  9  and  B  it  is  necessary  to 
obtain  a  relation  which  expresses  y  in  terms  of  these 
two  coordinates.  This  is  easily  handled  using  vector 
algebra. 

In  connection  with  Figure  13,  let  e^,  e^,  and  e^ 
be  unit  base  vectors  in  the  directions  X,  Y,  and  Z  respectively. 
Denote  a  unit  vector  in  the  direction  of  Vg  by  eg  and 
a  unit  vector  out  the  nose  of  the  bomb  by  efi.  Then 

eQ  ■  cosB  cos0  5^  +  cosB  sin©  e2  +  sinB  e"3  (30a) 


®B  ‘  °°»SB  *1  +  ,in6B  *2 


(30b) 
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Figure  13.  Coordinate  System  for  a  Regular 

Barricade  Enclosing  a  Single  Bomb. 
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Using  the  definition  of  the  cross-product , 

l«o  x  *inY  (31) 

Since  the  magnitudes  of  Sq  and  eB  are  both  one, 
this  relation  yields 

y  «  arcsin  (sin20  +  cos20  sin2(0  -  Og))1^2  *32* 

Let  V(y)  denote  the  number  of  fragments  per  steradian 
ejected  by  the  bomb.  Then  the  total  number  of  fragments 
ejected  out  an  arbitrary  region  is  given  by 

N  -  /  »(y)  dw  (33) 

where  the  element  of  steradian  du  is  given  by 

dw  *  cosB  d6  dB  (34) 

For  the  region  bounded  by  the  coordinates  0  ■  0 1 , 

0  *  02,  B  ■  8i  and  0  ■  02,  the  total  number  of  fragments 
would  be 

01,02  02  02 

N "  /  /  'My)  cos0  d0  d0  (35) 

6"  62  s,  e, 

If  the  average  number  of  fragments  per  steradian  over 
this  region  is  denoted  by  ¥12,  then  the  total  number 
of  fragments  would  be 

0  1  0  2 

Nft  '  -  ¥12  (0*  -  0i)  (sinB 2  -  Sin0,)  (36) 

p  1 ,  P  2 
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Usually  the  coordinates  of  a  region  such  as  the  one 
defined  above  are  given  by 

e2  »  e  +  A|  (37a) 

0 1  ■  0  —  —  (37b) 

2 

02  -  B  +  —  (37c) 

2 

0i*6—  —  (37d) 

2 

If  A0  and  A0  are  small  enough,  it  would  be  reasonable  to 
choose  as  average  values  for  the  initial  velocity,  number 
of  fragments,  and  fragment  mass,  those  values  given  in 
Figure  9  for  the  polar  angle  y(0*B)*  These  parameters 
could  then  be  used  in  connection  with  the  trajectory  equations 
derived  previously  to  determine  probable  impact  velocities 
and  coordinates. 

e.  Barricade  Geometry  Considerations 

To  include  barricade  geometry  affects,  a  more  general 
barricade  composed  of  straight  wall  segments  was  considered. 
This  introduces  very  little  additional  complicating  features 
and  allows  some  flexibility  so  that  optimization  of  barricade 
design  could  be  considered  in  the  future. 

The  geometry  of  the  wall  segments  of  the  barricade 
is  described  in  terms  of  cylindrical  coordinates  R,  0, 
and  Z  with  the  origin  placed  at  the  bomb  center.  It  is 
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assumed  that  the  wall  remains  intact  as  far  as  the  fragments 
are  concerned  so  that  if  the  fragment  strikes  a  barricade 
wall,  it  stops.  Thus,  for  particular  azimuth  and  departure 
angles/  the  existence  of  a  wall  for  that  azimuth  must 
first  be  determined;  if  a  wall  is  present  its  height 
must  be  known  so  that  the  question  of  whether  or  not  the 
fragments  have  cleared  that  portion  of  the  barricade  can 
be  answered. 

The  problem  to  be  solved  is  illustrated  in  Figure  14. 
For  a  given  azimuth  angle  of  trajectory  6/  the  distance  to 
the  wall  must  be  determined.  The  known  quantities 
are  R1#  R2 ,  81#  62/  and  the  height  of  the  barricade  z^, 
z2#  and  z3  Since  two  sides  (R^,  2)  and  the  included 
angle  (A8)  are  known,  the  law  of  tangents  can  be  used  to 
find  $  and  ¥ : 


(180  - 


|  (180 


R  —  R 

48)  +  tan'1  - A  tan  i  (180 

R1  +  R2 
R  —  R 

A8 )  -  tan"1  {  [— - -]  tan  i  (180 

*1  +  *2 


-  A8) }  (38a) 


-  A8) }  (38b) 


The  included  angle  between  R^  and  the  barricade  wall  is 


then 


6  «  180  -  <♦  +  (82  -  8 ) ) 
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The  law  of  sines  for  a  plane  triangle  can  then  be  used 
to  obtain  the  required  distance  to  the  barricade 


R3  ‘  R2  SiSy  ,40> 

f.  Distribution  of  Fragments 

The  previous  sections  have  outlined  a  method  for 
computing  the  impact  point  and  velocity  for  a  group  of 
fragments  which  are  assumed  to  have  identical  mass  and 
velocity.  The  total  number  of  fragments  in  this  group 
depend  on  the  azimuth  angle  6,  departure  angle  $  and  the  size 
of  the  region  defined  by  0i,  0 2 >  61  and  02.  In  actual 
fact  all  of  these  fragments  will  not  land  at  one  spot 
but  in  general,  they  will  be  distributed  over  some  area. 

The  following  discussion  presents  a  method  that  should 
give  reasonable  results  that  can  be  compared  with  exper¬ 
imental  data. 

For  a  given  0,  a  series  of  impact  ranges  will  be 
determined  together  with  corresponding  values  of  mass, 
number  of  fragments  and  impact  velocities.  The  number  of 
impact  ranges  will  correspond  directly  to  the  number  of 
increments  used  to  cover  the  range  of  0.  Suppose  the  ranges 
are  ordered  in  an  increasing  sequence 

0  £  R!  1  r2  i  •••  £  Rj  i  •••  £ 

To  illustrate  the  procedure,  consider  the  fragment  parameters 
associated  with  the  range  R^.  A  reasonable  approximation 
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is  to  distribute  these  fragments  over  the  area  bounded 
by  the  coordinates  2>  R2  3#  8t  and  02(See  Figure  15) 
where 


0i-0- 

0  2  —  9  + 


A0 

“7 

A0 


R1  +  R2 


1,2 


R2  +  R3 


2,3 


(41a) 

(41b) 

(41c) 

(41d) 


The  area  A2  covered  by  such  a  segment  is 


)  A0 


(42) 


Similarly,  the  area  associated  with  the  range  R_ 

J 

would  be 


AJ  “  7  (RJ,J+1  "  RJ-1,J)  A0 

where 

RJ-lfJ  *  7  ^RJ-1  +  RJ^ 
RJ,J+1  “  7  (RJ  +  RJ+1) 


(43) 


(44a) 


(44b) 


For  the  furthest  impact  point  R^,  choose 


*N-1,N  “  J  (RN-1  +  *N) 

(45a) 

RN,N+1  2rn  *^-1^ 

(45b) 

With  the  impact  areas  defined  by  Equation  (43) , 
both  the  number  of  fragments  and  the  total  weight  per 
unit  area  can  be  determined  by  dividing  the  total  number 
of  fragments  and  the  total  weight  landing  at  a  particular 
range  by  the  corresponding  area,  that  is. 


fonTtTrea]  “  ^Numbor  °*  fragment*  landing  at  Rj)/Aj  (46a) 

(Weight  )  m  f Number  1  (Average  weight  of  fragment 
[unit  AreaJj  [unit  AreaJj  landing  at  Rj)  (46b) 
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5.  ANALYTICAL  FRAGMENT  MODEL  FOR  A  STACK  OF  BOMBS 

The  fragmentation  model  developed  in  the  preceding 
diacuaaion  will  yield  fragment  dispersiona  for  a  single  bomb. 

The  theory  ia  expected  to  hold  for  a  stack  of  bomba  where 
the  origin  of  the  coordinate  system  is  placed  at  the  center  of 
the  stack.  However,  interaction  effects  should  produce  a  number 
of  fragments  for  a  given  impact  area  somewhat  less  than  the 
product  of  the  number  of  bombs  and  the  number  of  fragments 
produced  by  a  single  bomb.  The  exact  effect  must  be  determined 
experimentally . 

An  illustration  of  the  fragment  survey  areas  from  Phase  II 
of  the  BIG  PAPA  Tests  is  shown  in  Figure  16.  For  each  one  of 
the  fragment  survey  areas,  fragments  were  counted  and  weighed. 

By  computing  the  number  of  fragments  and  their  weights  for  the 
same  area,  a  "correlation”  or  "effective  number  of  bombs”  factor 
may  be  computed  for  each  e**ea.  Hopefully,  the  variation  of  the 
correlation  factor  from  survey  area  to  survey  area  would  lie 
within  a  tolerable  limit.  With  the  "effective  number  of  bombs" 
EFNB,  determined,  the  fragment  distribution  for  a  single  bomb  can  be 
nultiplied  by  EFNB  to  obtain  the  distribution  pattern  for  the  stack  of 
bombs. 
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This  section  has  presented  a  model  that  predicts  the  impact 
velocity  and  distribution  of  fragments  for  a  barricaded  or  unbarri¬ 
caded  boob  stack. 

Data  in  the  form  of  initial  fragment  velocity,  size  and  dis¬ 
tribution  for  each  bomb,  which  is  available  in  References  5  and 
6,  is  used.  The  fragments  thr.t  are  ejected  within  a  small  area 
of  the  bonbcase  are  assumed  to  have  identical  initial  conditions 
as  far  as  the  trajectory  equations  are  concerned.  Average  values 
for  the  mass  and  initial  velocities  of  the  fragments  are  used 
together  with  a  ballistic  coefficient  that  was  determined  experi¬ 
mentally.  if  the  group  of  fragments  clear  the  barricade  then  the 
point  and  velocity  of  impact  can  be  determined.  For  comparison 
with  experimental  data,  this  group  of  fragments  is  assumed  to  be 
distributed  over  an  adjacent  area  determined  by  the  increments 
in  azimuth  angles  and  points  of  impact  of  other  groups  of  fragments 
along  the  same  azimuth  angle. 

The  same  procedure  is  used  for  a  stack  of  bombs  except  that 
the  total  number  of  bombs  must  be  modified  to  account  for  inter¬ 
action  effects.  The  degree  of  modification  must  be  determined 
emperically  from  existing  experimental  data. 

The  computer  program  that  follows  the  theory  of  this  section 


is  described  in  Appendix  II. 


SECTION  IV 


CRATERING 


1 .  INTRODUCTION 

A  considerable  amount  of  data  have  been  accumulated  that 
relate  the  size  of  craters  to  the  yield  of  explosives ,  primarily 
for  bare  spherical  and  hemispherical  charges.  Most  of  the  work 
accomplished  prior  to  1961  has  been  summarized  in  Reference  12. 
Crater  measurements  of  later  detonations  of  major  importance 
have  been  made  and  include  "Operation  Snow  Ball"  (Reference  13) , 
"Operation  Distant  Plain"  (References  14  and  15)  and  "Operation 
Sailor  Hat"  (Reference  16).  These  shots  cover  a  wide  range  of 
explosive  yield  and  soil  types,  a  fact  which  produced  a  con¬ 
siderable  amount  of  diversity  in  crater  sizes. 

Predicting  the  shape  of  explosion-produced  craters  and 
the  distribution  of  the  ejected  material  has  been  a  matter  of 
some  concern  for  several  years  (References  17-20) .  Generally, 
the  approach  has  been  to  assume  a  non-dimensional  relation 
between  characteristic  dimensions  of  the  crater  and  the 
explosive  yield  raised  to  some  power.  The  resulting  relation, 
called  a  scaling  law,  is  then  used  to  extrapolate  to  new  regions 
of  interest. 

These  scaling  laws,  even  though  set  up  to  form  a  "best 
possible"  fit  for  existing  data,  must  be  modified  for  changes 
in  earth  media.  Furthermore,  it  has  become  rather  apparent 
from  the  shots  at  the  Suffield  Experimental  Station  and  from 
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"Operation  Sailor  Hat"  that  any  one  scaling  law  will  produce 
reasonable  results  over  a  limited  yield  range. 

The  charge  shapes  in  the  tests  mentioned  above  are  exclu¬ 
sively  spherical  or  hemispherical.  Results  for  other  shapes 
such  as  rectangular  parallelepipeds  that  would  be  of  more  sig¬ 
nificance  for  this  report  are  almost  non-existent  except  for 
the  ,fBIG  PAPA"  test  (Reference  1) . 

In  light  of  the  almost  total  absence  of  analytical  work 
in  this  area  and  because  of  the  different  shaped  charge  of 
interest  in  this  project,  specifically  that  associated  with 
bomb  stacks,  it  seemed  appropriate  to  develop  an  elementary 
model  using  basic  principles  of  mechanics.  Although  the 
approach  adopted  in  this  section  freely  uses  past  empirical 
relations,  such  an  analysis  could  form  the  basis  of  a  more 
rigorous  development  in  the  future  if  experiments  could  be 
devised  to  adequately  determine  the  governing  parameters. 

As  for  the  past  sections,  all  charges  will  be  assumed  to 
be  resting  on  the  surface. 


THE  EFFECT  OF  CHARGE  SHAPE 


a.  Preliminary  Comments 

When  an  explosion  is  detonated,  the  total  available 
energy  is  divided  into  various  categories!  blast  wave, 
heat,  and  kinetic  energy  of  the  material  itself  to  mention 
the  most  obvious  ones.  Apparently  the  blast  wave  does  not 
contribute  to  cratering,  but  rather,  it  causes  a  shock 
wave  to  be  instigated  in  the  earth  (Reference  21).  The 
major  source  of  cratering  action  must  then  be  the  kinetic 
energy  of  the  explosive  material.  By  momentum  transfer, 
energy  is  transferred  to  the  earth  media  elements  located 
on  the  surface  and  adjacent  to  the  explosive,  and  by 
propagation  within  the  neighboring  region,  earth  particles 
are  ejected  and  a  crater  is  formed. 

Because  of  f-he  confining  effect  of  the  explosive 
material  on  itself,  it  seems  plausible  that  the  initial 
direction  of  propagation  of  explosive  elements  would  be 
i.owards  the  surface  of  the  explosive.  Furthermore,  because 
of  air  friction,  interaction  with  other  elements  and  so 
forth,  the  velocity  of  each  element  would  decrease  with 
time.  Hence,  those  elements  closest  to  the  earth  would  be 
the  most  effective  as  far  as  cratering  is  concerned. 

b.  Basic  Assumptions 

The  above  observations  suggest  that  it  is  appropriate 
to  make  several  simplifying  assumptions  to  make  an  investi- 
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gat  ion  of  charge  shape  amenable  to  analysis.  Because  of  the 
preliminary  nature  of  this  work,  the  following  assumptions  are  made 
with  the  full  knowledge  that  they  may  be  unjustified;  however,  the 
resulting  analysis  should  yield  a  reasonable  approximation: 
approximation : 

(1)  Assumption  1 

Just  the  bottom  half  of  the  charge  contributes  to 
the  cratering  phenomenon. 

(2)  Assumption  2 

The  velocities  of  all  elements  in  the  bottom  half 
of  the  charge  are  the  same  immediately  after  detonation 
and  are  directed  vertically  downward.  This  implies  that 
if  \Q  is  the  initial  speed  of  all  elements  in  the  charge, 
then  the  initial  kinetic  energy  per  unit  mass 

«o  -  7  vo  <47> 

is  independent  of  position. 

(3)  Assumption  3 

The  friction  force  is  constant  and  is  the  same  with” 
in  and  outside  the  original  outline  of  the  charge. 

Such  an  assumption  yields  a  couple  of  interesting 
results.  First  it  can  be  shown  that  the  kinetic  energy/unit 
mass  e  decays  linearly  with  distance  by  letting  z  denote  the 
distance  of  an  element  above  the  surface  at  any  time  t 
(See  Figure  17).  From  Newton's  Law,  the  acceleration  is 
constant 

5s  -  aQ  (48) 
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so  that  the  velocity  and  position  respectively  becone 

v  ■  t  ■  at  -  v  (49a) 

o  o 


V 


+  V  t  +  Z 

o  o 


By  a  series  of  substitutions 


M- 


lo*o 


ao7 


(49b) 


(50) 


and  hence ,  e  decreases  as  z  decreases. 

Secondly,  if  an  element  is  initially  at  the  critical 
distance  L  above  the  surface  where 

'o 

L  -  —  (51) 

2®o 

then  v  (and  hence,  e)  is  zero  when  that  particular  element 
reaches  the  surface.  In  other  words,  all  elements  above-* 
the  plane  z«l.  tnd  in  the  bottom  half  of  the  charc/e  will 
not  contribute  to  the  cratering  action. 

(4)  Assumption  4 

For  the  range  of  explosive  yield  considered  in  this 
report,  the  height  of  the  center  of  gravity  of  the  charge 
zCG  is  below  the  plane  z-L  so  that  all  parts  of  th«i  bottom 
half  of  the  charge  will  contribute  to  cratering. 
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(5)  Assumption  5 


At  the  interface  between  the  surface  of  the  earth 
and  the  explosive ,  the  loss  of  kinetic  energy  is  negligible. 


c.  Cratering  Factor 

The  energy  es  per  unit  mass  delivered  at  the  surface  of  the  earth 
by  elements  originally  at  a  distance  zQ  above  the  surface  is  obtained 
from  Equation  (50) : 


¥1 


(i 


(52a) 


or,  after  using  Equation  (51) 

es  *  eo  (1  *  VL> 

The  total  kinetic  energy  of  the  charge  is 


I  eo 

B 

^O 


(52b) 


(53) 


where  M  denotes  the  total  mass  of  the  charge  and  B  the 
region  occupied  by  the  total  charge.  On  the  other  hand, 
the  kinetic  energy  delivered  to  the  surface  of  the  earth 
is 


dm 


(54) 
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where  denotes  the  bottom  half  of  the  charge  and  zCQ 
is  the  distance  to  the  center  of  mass  of  the  total  charge. 

A  cratering  factor  Cp  can  be  defined  as  the  ratio  of 
the  kinetic  energy  reaching  the  earth's  surface  to  the  total 
kinetic  energy ,  that  is 

CP  1  V*T 


According  to  assumption  4,  the  cratering  factor  must  lie 
in  the  range 


(56) 


Equation  (55)  Implies  that  for  two  charges  with  the 
same  total  kinetic  energy,  the  charge  with  the  lowest  center 
of  mass  will  be  the  most  effective  as  far  as  cratering  is 
concerned. 

d.  Charge  Shapes  and  Non-dimensional  Variables 

If ,  W  denotes  the  yield  of  an  explosion  in  terms  of 

an  equivalent  weight  of  TNT,  then  immediately  after 

detonation,  assume  that  the  total  kinetic  energy  of  the 

explosive  i®  directly  proportional  to  W  and  is  independent 

of  the  charge  shape,  that  is 

Et  -  wW  (57) 
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Wh.re  V  i.  a  con.tant.  The  energy  delivered  to  th.  .urf.c. 
is  then 


Ea  -  Cp  VW 


(58) 


To  eliminate  the  unknown  factor  u,  reference  Chargee 

can  be  introduced.  For  each  charge  ahape,  let  w  denote 

o 

a  reference  yield  and  define 


Es  =  CFo  ""0 


(59) 


where 


'Fo 


(60) 


and  rce  is  the  height  of  the  center  of  meet  of  the  reference 
charge . 

If  the  following  non-dimensional  variables  are 
introduced : 


T  -  w/w 


(61a) 


K  -  Ea/E” 


CF/CFc 


(61b) 

(61c) 


then 


K  -  cj  T 
o 


(62) 
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and  tha  factor  u  is  not  present.  If  zCG  and  z°G  are  both 
much  smaller  than  L,  then  is  approximately  equal  to 
one. 

For  a  given  class  of  charge  shapes,  the  total  yields 
cam  be  used  rather  than  the  dimensions  L,  z__  and  z^L. 

Ui  Cv 

One  such  eet  is  illustrated  in  Figure  18  and  is  eharacter- 
ized  by  the  fact  that  the  volume  is  proportional  to  the 
cube  of  one  dimension  in  each  case.  This  restricts  the 
group  of  triangular  and  cylindrical  prisms  (which  includes 
rectangular  parallelepipeds)  to  that  for  which  the  surface 
contact  area  Afi  is  proportional  to  the  square  of  the  height 
h.  In  effect  this  implies  that,  for  exasiple,  in  the  case 
of  triangular  prisms,  we  can  let  the  sizes  change  but  the 
shapes  must  be  similar  to  the  reference  shape. 

Since  the  yield  is  directly  proportional  to  the  vol¬ 
ume,  it  can  also  be  said  that  for  this  class  of  shapes 
the  distance  to  the  center  of  mass  is  proportional  to  the 
cube  root  of  the  yield.  Thus,  if  we  let 


W. 


W 


CG 


(63) 


where  WL  is  the  weapon  yield  for  a  charge  shape  in  the  same 
class  as  those  in  Figure  18  and  whose  center  of  mass  is  at  a  height 
L  above  the  surface,  then 
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each  of  the  charge  shapes  shown  in  Figure  18. 


\ 
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3.  CRATER  AND  EJECTA  FORMATIONS 
a.  Basic  Shape  Parameters 

The  most  significant  parameters  associated  with  the 
description  of  the  crater  and  ejecta  shape  are  shown  in 
Figure  19.  If  the  origin  of  a  cylindrical  coordinate 
system  is  placed  on  the  original  surface  at  the  center  of 
the  crater ,  then  the  crater  depth  is  assumed  to  be  adequately 
described  by  the  parabola 

#  ■  a  +  bR  +  cR2  (65) 

where  a,  b  and  c  are  constants  and  R  is  the  distance  from 
the  origin.  If  denotes  the  apparent  depth  of  the  crater 
at  the  origin,  R  the  apparent  radiuc  at  the  original 

A 

surface  and  if  we  assume  that 


da 

dR 


R-0 


(66) 


then  the  crator  depth  is  described  by 

0  1  R  1  Ra 


Da  |1  - 


(67) 


The  maximum  slope  of  the  crater,  which  will  be  of  sig¬ 
nificance  later  in  connection  with  energy  dissipation, 
occurs  a*  the  intersection  of  the  crater  with  the  original 
surxw.ce  and  is  given  by 


da 

cT£ 


R-R. 


(68) 
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An  aid  in  visualizing  this  condition  is  to  note,  as  shown 
in  Figure  20,  that  a  cavity  with  this  slops  svsrywhsra 
would  be  a  cone  with  twice  the  depth  of  the  actual  crater. 

After  a  detonation,  particles  that  were  originally  within 
the  boundary  of  the  crater  are  located  immediately  beyond  the 
top  rim  of  the  crater  on  the  surface  of  the  earth  and  this 
material  that  has  been  thrown  out  is  called  the  ejecta.  The  depth 
of  the  ejecta  as  a  function  of  distance  from  the  center  of  the 
crater  was  assumed  to  be  the  same  as  that  given  in  Reference  22  for 
nuclear  detonations.  This  relation  is 


{>} 


R>  R. 


where  a  and  6  are  parameters  that  depend  on  the  earth  media. 

If  the  earth  media  is  assumed  to  be  incompressible,  then  the 
volume  of  the  crater  should  be  the  same  as  the  volume  of  the 
ejecta.  Such  a  relation  can  be  used  to  express  o  in  terms  of 


The  volume  of  the  crater  is 

R  2it  _  a 


/*./.  / 


R  dx  d6  dR 
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or 


TT  2 
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(70) 


Similarly,  if  8  >  2,  the  volume  of  the  ejecta  is  given 
by 


|  |  |H*  R  dz  d0  dR 


R  o 


2tt 


r 


a  Da  R8  R1"8  dR 


or 


— 0  D.  *1 


6-2 


(71) 


Equating  the  two  volumes  yields  the  relation 


o  «  j  (8  -  2) 


(72) 


so  that,  from  Equation  (69) 


Da  *  fR.l 8 

He  -  ^  (S  -  2)  K* 


(73) 


Reference  22  suggests  the  values  a  *  0.5  and  8  ■  3.9 
for  soil,  and  a  *  0.3  and  8  ■  3.1  for  rock.  However, 
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Equation  (72)  yields  the  value  6  ■  4.0  when  o  *  0.5  and 

A 

0  *  3.2  when  a  ■  0.3.  Such  a  variation  is  negligible  in 
view  of  the  disparity  in  test  results. 

Of  more  significance  are  experimental  values  for  the 
height  of  the  crater  lip.  From  Equation  (73)  the  analytical 
expression  is 


(74) 


According  to  results  tabulated  by  Vortman  (Reference 
16),  the  100**ton  Suffield  Experimental  Station  hemispherical 
shot  in  clay  yielded  a  crater  lip  height  which  was  27 
per  cent  of  the  crater  depth,  which  produces  a  value  of  3.1 
for  fL  On  the  other  hand,  for  the  500-tcn  Sailor  Hat  shot 
on  basalt  rock,  the  height  of  the  crater  lip  was  36  per  cent 
cf  the  crater  depth,  which  yields  a  value  of  3.4  for  B. 

These  results  imply  that  the  values  of  0  do  not  assume 
the  same  range  of  values  for  conventional  high  explosives 
as  for  nuclear  explosives.  Furthermore,  for  conventional 
explosives  the  variation  in  0  may  be  quite  small  for  changes 
in  earth  media. 

For  this  project,  a  reference  value  of  3.1  was  chosen 
for  0  for  soil. 

It  seems  plausible  to  assume  that  most  of  the  material 
will  move  radially  outward.  For  the  next  section  it  is 
necessary  to  know  the  positions  of  the  centroids  of  crater 
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and  ejecta  elements  that  aubtand  a  small  ngl.  ,e  (8#. 
figure  21) .  The  crater  and  eject,  volume  element,  are 


*  AV_ 

e 

(75) 

The  coordinates  to  the  centroids  of  these 

elements  are 

defined  as  follows: 

AVe  R„  *  A0 
t  c 

/*•  Is  *2  ax  dR 

(76a) 
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The  results  for  g  >  3.0  are 


nMMHi | 


CRATER  RIM 


TOP  VIEW 


SIDE  VIEW 


CENTROID  OF 
CRATER  ELEMENT 


CENTROID  OF 
EJECTA  ELEMENT 


CENTROID  OF 
.EJECTA  ELEMENT 


CENTROID  OF 
CRATER  ELEMENT 


Figure  21.  Centroids  of  Crater 
and  Ejecta  Elements 
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As  an  example,  the  centroid  for  the  ejecta  for  6  *  3.2 
is  given  by 


Re  - 
e  6-3.2 


(78a) 
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6-3.2 
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b.  Energy  Considerations 

Suppose  an  element  with  initial  velocity  vq  moves  a 
horizontal  distance  X  and  a  vertical  distance  H  as  shown 
in  Figure  22.  If  just  the  effect  of  gravity  is  considered 
in  the  equations  of  motion,  then  from  the  equation  of  the 
trajectory  it  can  be  shown  that 

2 

H  -  -  2 - - -  +  X  tan6  (79) 

v*  cos^6 
o 

For  a  conservative  estimate  on  the  energy  requirement, 

choose  6  with  H  and  X  considered  fixed  such  that  v  i r 

o 

a  minimum,  that  is,  set 
dv 

XT  ”  0  (80) 
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This  yields 


which  can  be  written  in  the  alternate  ton 


1 

- 3T 

cos  8 


«v  *  < 


(81) 


(82) 


By  substituting  Equations  (81)  and  (82)  into  (79)  we  get 
vo  -  gH  Cl  ±  /l  ♦  X2/H2  )  (83) 

Only  the  positive  sign  is  appropriate.  It  X/H  <  <  lf 
%#e  have  the  classical  vertical  Motion  relation 

Vq  *  2gH,  X/H  <  <  1  (84) 

whereas ,  it  X/H  >  >  1,  we  get 

vo  *  gX,  X/E  >  >  1  (85) 

Suppose  that  on  the  average,  an  element  Moves  from 
the  centroid  of  a  crater  eleMent  to  the  centroid  of  an 
ejecta  element.  Then 


(86a  » 


H  «  x  ♦  s 

c  c 


1  <8  -  2)21 

16  (8  -  1)  j 


(86b) 


For  the  range  3.1  <  8  <  3.6,  the  seal  lest  value  of  X/H 
occurs  at  §  *  3.6.  Bence 


5.5 


3.1  <  8  <  3.6 


(87) 


For  actual  shots,  it  is  generally  true  that  the  apparent  radius  is 
at  least  twice  the  apparent  depth  so  that 

|  >  11  (88) 

Because  of  the  diversity  of  experimental  data,  the  inequality 
associated  with  Equation  (85)  can  be  considered  satisfied 
and  hence,  with  the  use  of  Equation  (86a) 


(89) 


In  order  that  this  expression  always  be  positive,  we  must 
have  8  >  3  which  is  the  same  restriction  imposed  previously 
in  connection  with  centroids. 

If  the  mass  density  of  the  earth  media  is  denoted  by 
p,  the  total  initial  kinetic  energy  of  the  earth  elements 


is  approximated  by 


fVo 


2 

3 


(90) 
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It  should  be  emphasized  that  E„  is  not  the  same  as  the 

Vj 

energy  input  E  since  a  portion  of  the  latter  will  be 

B 

dissipated  into  the  ground  in  the  form  of  heat. 

c.  Dimensional  Considerations 

To  describe  the  shape  of  the  crater  it  is  necessary 
to  know  the  ratio  D  /R  .  Postulate  that  this  ratio 

cl  cl 

primarily  depends  (Reference  23)  on  the  following  par¬ 
ameters:  (1)  an  earth  media  viscosity  v',  (2)  mass 

density  of  the  earth  media  p,  (3)  the  kinetic  energy 
input  E  ,  and  (4)  the  apparent  crater  radius  R  .  Such  a 

8  A 

dependence  can  be  expressed  analytically  by 

Da/Ra  *  f'  (p'  Es '  v''  V  (91) 

I 

where  f  is  the  unknown  function.  If  [M] ,  [L]  and  [T] 
denote  the  fundamental  dimensions  of  mass,  length  and  time 
respectively,  then  the  dimensions  of  the  parameters  in 
Equation  (91)  are 


‘VEa 

]  *  1 

(92a) 

[p]  = 

[M]/ [L] 3 

(92b) 

IES]  - 

[M]  [L]2/[T]2 

(92c) 

[R  ]  - 

Cl 

[L] 

(92d) 

[v  ]  = 

[M]/[T]  [L] 

( 92e) 
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According  to  the  Buckingham  ir-Theorem,  the  latter  four 
variables  can  be  combined  into  one  non-dimensional  variable 
which  is  chosen  to  be 


<P 


v2»a 

E_ 


(93) 


where 


v  «  v'//p  (94) 

can  be  considered  a  generalized  viscosity.  Now  Equation 
(91)  can  be  given  as 


(95) 


where  f  is  a  non-dimensional  function  of  the  parameter 
4>  and  is  unknown.  Since  no  analysis  in  this  connection 
appears  to  be  available  as  a  guide  in  choosing  a  suitable 
form  for  f,  assume  a  simple  exponential  relation  of  the 
fo] lowing  type : 


vr«  ■  k*c 


(96) 


where  k  and  ;  are  constants.  Note  that  this  relation  will 
yield  the  maximum  slope  of  the  crater  wall  with  the  use 
of  Equation  (68) . 
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Suppose  that  for  one  earth  media,  values  of  parameters 
associated  with  a  reference  charge  are  denoted  by  a  super¬ 
script  zero.  Then,  according  to  Equations  (90)  and  (96) 

ES-  Jp«  *£<*$> 3  [M-1-]  w.) 

*  *  -  3  is' 


(97b) 


The  parameters  p,  §  and  v  depend  only  on  the  type  of 
earth  media  and  not  on  the  size  of  charge.  Hence,  these 
variables  do  not  have  the  superscript  sero.  On  the  other 
hand,  it  is  assumed  that  k  and  C  are  independent  of  both 
earth  media  and  charge  size. 

By  taking  appropriate  ratios  of  the  terms  in  Equations 
(90) ,  (96)  and  (97)  the  following  relations  are  obtained: 


(98a) 


(98b) 


These  equations  do  not  contain  the  terms  that  depend  on 
the  earth  media;  hence,  they  can  be  used  for  predicting 
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results  from  the  knowledge  of  one  surface  detonation  on 
a  given  soil  or  rock. 

The  following  analysis  is  developed  to  relate  the 
actual  kinetic  energy  of  the  earth  particles  E_  to  the 
kinetic  energy  delivered  to  the  surface  Eg. 

d.  Energy  Dissipation 

The  initial  kinetic  energy  EG  of  the  earth  particles 
will  differ  from  the  kinetic  energy  delivered  to  the  sur¬ 
face  Eb  by  the  amount  of  energy  ED  dissipated  in  the  form 
of  heat: 


E  -  E 

s  o 


(99) 


From  the  expression  in  the  incremental  theory  of  plasticity 
(Reference  24)  for  energy  dissipation  where  t  denotes 
stress  and  £  strain  rate 


t i  -  vV 


(100) 


it  seems  appropriate  to  assume  that  ED  is  linearly  pro- 

9 

portional  to  v  and  also  a  function  of  E  ,  p  and  R  , 

•  A 

that  is 


ED  -  V  *  (Eg,p,Ra) 


(101) 


Dimensional  homogeneity  for  the  equation  implies  that 


Ed  -  M, 


A.l1/2 


(102) 


where  b  ii  taken  to  be  a  non-dimensional  constant  and  v 


is  defined  by  Equation  (94).  For  a  reference  charge,  we 
get 


b  E 


E. 


so  that 


!• !» 
ES  *1 


1/2 


(103) 


(104) 


For  the  case  of  a  reference  charge,  it  is  convenient 
to  introduce  a  dissipation  ratio 


e: 


E' 


(105) 


which  is  simply  the  ratio  of  the  energy  dissipated  to  the 
energy  available  at  the  surface.  Hence  E^  must  assume  a 
value  between  0  and  1. 

From  Equation  (99) 


(106) 
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Substitute  Equation*  (104),  (105)  and  (61b)  into  (106) 
to  qet 


eg*e.  1 


E®  [  Ra)l/2 

■  1  ‘  —  K  -4 

L  E  (  *l\  J 


[-  N  '1 


(107) 


Also 


eg  "  <  »  -  # 


(108) 


so  that 


(i  -  E°  [— 

L 


(109) 


Hence,  with  the  use  of  Equation  (61b),  Equations  (98a) 
and  (98b)  can  be  written  in  the  alternate  form 


Ra) 


l  _  pD 

[  V 

1/5] 

s 

KR° 

J 

j 

I?-e3 


(110a) 
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(110b) 


With  R° ,  D° ,  c  and  E^  presumed  known,  these  two  equations 
give  the  apparent  radius  and  depth  as  a  function  of  K  and 
hence  as  a  function  of  the  yield  of  the  explosive.  Then, 
if  0  is  known.  Equation  (73)  can  be  used  to  predict  the 
ejecta  depth. 

For  a  given  earth  media  (alluvium  for  example),  p 

A 

and  0  are  assumed  known.  For  some  reference  energy  e”, 
the  apparent  radius  and  depth,  R®  and  D°  respectively  cm  . 
be  measured,  and  values  for  C  and  determined  experimentally. 
A  different  earth  media  will,  in  general,  yield  different 
values  for  these  parameters  which  will  be  designated  by  an 
asterisk.  For  a  new  reference  energy  E*°,  R*°  and  D*® 
can  be  measured  as  before.  The  new  density  p  can  also 
be  determined  experimentally  and  it  is  assumed  that  C  remains 
unchanged.  The  parameters  0*  and  E*^  can  be  determined 
according  to  the  following  analysis. 

From  Equations  (103)  and  (105) 

1/2 

v  (111a) 


.V. 


(111b) 


or,  after  eliminating  b 


(112) 


In  a  eimilar  manner,  Equation  (97b)  yields 


*o 


i+c 

E° 

c 

*1 

a 

s 

V 

k  1 

8 

.  v  . 

(113) 


By  eliminating  v  /v  between  Equations  (112)  and  (113), 
we  get 


S*D  -  eD 

-  8  "  ES 


1 

"5? 


(114) 


To  obtain  a  value  for  8*,  Equations  (97a)  and  (108) 
for  the  two  media  can  be  combined  to  get 
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Equation  (59)  can  be  ug.d  to  expreaa  E*°/E°  in  term,  of 
reference  charge  yield*  and  shape* t 


*  * 
CF0  Wo 


CF0  Wo 


\ 


(116) 


4. 


For  may  rafaranoa  in  concoction  with  tha  computer  proyraa, 
th«  pertinent  equations  will  he  suwaarised. 

For  ona  partueul.r  aarth  madia,  *o  danota.  a  rafaranoa 
char,,  for  which  tha  apparant  radio.  E°  and  app.r.nt  dapth  D 
ar.  known,  toy  othar  yi.ld  i.  .xpr.a.ad  i„  tan.,  of  tha  raf- 
ar.no.  chary,  by  Kan.  of  tha  non-dim.n.ion.1  paramatar 


T  «  2- 

wo 


(117) 


*h.  anaryy  E,  avails.  ,t  tha  aurfac  i.  .1.0 
dimensionally  by  swans  of  the  factor 


K-fi 


(118) 


where 


K  -  C'  T 

and 


(119) 


(120) 
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For  lack  of  a  pracisa  value,  WL  la  takan  to  ba  106  lbs. 
or  500  tons  of  TNT. 

With  the  dissipation  ratio  E®  for  the  reference  charge  Wq 
and  the  particular  earth  media  assumed  known  fro*  experimental 
sources,  the  apparent  radius  and  depth  for  charges  of 
various  values  are  determined  from 


(121a) 


(121b) 


The  parameter  c  is  assumed  to  be  the  same  for  all  earth  media 
and  is  chosen  so  that  the  theoretical  results  fit  the  exper¬ 
imental  data  as  closely  as  possible. 

With  6  known,  the  depth  of  ejecta  is  given  by 


(122) 


One  set  of  basic  reference  parameters  are  those  associated 


with  the  100  ton  shot  at  the  Suffield  Experimental  Station 


According  to  Reference  16,  the  coil,  a  silty  clay,  had  a  weight 

density  o£  94  lb/ft3  and  the  apparent  depth  and  radius  were 

21  and  70  feet  respectively.  As  mentioned  previously,  an 

.  appropriate  value  for  0  for  soil  is  3.1.  A  value  of  0.3  for 

both  C  and  E®  for  this  reference  charge  and  earth  media  appears 
to  give  reasonable  results. 

The  above  equations  are  also  used  for  a  different  earth 

■edia.  However,  the  new  apparent  depth  D*°  and  apparent 

radius  R  ®  must  be  determined  experimentally  for  the  new  reff- 
.  * 

erence  charge  WQ.  The  new  dissipation  ratio  is  given  by 


and  the  new  ejecta  parameter  0  can  be  found  from 


-  2 


-  3 
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(124) 


where 


(125a) 


(125b) 
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For  a  discontinuous  basalt  rock  of  weight  density  190 
lb/ft3,  Reference  16  states  that  a  charge  W*  -  500  tons  of 

TNT  yields  an  apparent  radius  R  ”  ■  79  ft.  and  an  apparent 

*o  *D 

depth  D  »  38  ft.  These  values  can  be  used  to  obtain  E 

«  m 

*  * 

and  6  for  this  type  of  rock.  Then  equations  (121a)  and 
(121b)  with  the  new  parameters  can  be  used  to  find  the  dim¬ 
ensions  of  a  crater  for  any  other  charge  below  500  tons  of 
TNT. 

The  above  outline  briefly  describes  a  theory  that  should 
adequately  predict  the  crater  and  ejecta  shapes  for  a  wide 
range  -»f  yield  for  conventional  explosives.  The  effect  of 
charge  shape  is  included  in  the  analysis. 

For  bomb  stacks,  the  theory  of  this  section  appears  to 
be  satisfactory  if  only  the  weight  of  the  high  explosive  is 
used.  For  stacks  that  are  of  the  same  order  of  magnitude  as 
the  barricade,  the  effect  of  the  barricade,  as  far  as  crater 
and  ejecta  shapes  are  concerned,  is  assumed  to  be  negligible. 
For  cases  where  the  bomb  stack  is  relatively  small,  the  sit¬ 
uation  is  quite  different  and  for  most  practical  purposes,  the 
crater  and  ejecta  shapes  are  not  too  significant. 

The  computer  program  associated  with  the  theory  of  this 
section  is  described  in  Appendix  III. 
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SECTION  V 


ILLUSTRATIVE  EXAMPLES 

1.  INTRODUCTION 

In  this  Motion/  a  representative  sat  of  curves  is  pre- 
santad  List  is  based  on  the  theory  of  the  previous  secticns 
and  the  associated  computer  programs.  Because  of  the  large 
number  of  parameters  that  are  present,  only  typical  values 
were  chosen  for  a  graphical  representation  of  the  computer 
output.  These  curves  are  intended  just  to  illustrate  the 
type  of  information  that  can  be  obtained.  In  many  instances 
the  computer  output  is  more  detailed  and  can  handle  several 
possible  situations  which  are  not  appropriate  for  a  graphical 
display. 

Some  of  the  results  of  available  experimental  data  are  also 
plotted  to  indicate  the  degree  of  correlation  between  predicted 
and  actual  values.  As  stated  previously,  a  certain  amount  of 
judgement  is  necessary  when  using  these  programs. 


I 

\ 

2.  RESULTS  FROM  BLAST  PRESSURE  PROGRAM 

As  a  typical  example,  parameters  associated  with  Phases  I 
and  II  of  the  "BIG  PAPA"  (Reference  1)  tests  were  used  as  input 
data  for  the  program.  The  bomb  stack  dimensions  were  30  ft. 
wide  by  50  ft.  deep  by  8.83  ft.  high  and  the  corresponding 
barricade  dimensions  were  100  ft.  by  70  ft.  by  11  ft.,  respectively. 

Associated  with  a  bomb  stack  are  several  conversion  factors 
which  are  listed  in  Reference  3.  The  first  replaces  the  weight 
of  the  explosive  in  a  bomb  by  an  equivalent  weight  of  TNT. 

According  to  Reference  1,  which  used  a  factor  of  1.23  in  convert¬ 
ing  tri tonal  to  TNT,  the  bomb  stacks  contained  an  equivalent 
weight  of  307,500  lbs.  of  TNT.  From  Reference  3,  a  factor  of 
0.6  was  considered  most  appropriate  for  replacing  the  bomb  stack 
by  a  bare  charge  (184,000  lbs.)  that  would  yield  the  same  blast 
pressure  characteristics.  A  different  factor  must  be  used  for 
impulse  but  since  the  procedure  is  quite  analogous,  impulse 
distributions  are  not  given. 

Pressure  versus  distance  is  plotted  in  Figure  23  for  the 
following  cases: 

a.  Bare  hemispherical  unbarricaded  charge, 

b.  Bare  hemispherical  barricaded  charge, 

c.  Rectangular  unbarricaded  charge,  and 

d.  Rectangular  barricaded  charge. 

Figure  23a  shown  rather  predictable  results  for  the  effect 
of  the  barricade  on  the  blast  pressure  from  a  hemispherical 
charge.  For  a  given  distance  from  the  center  of  the  charge, 
the  pressure  out  the  front  or  open  end  of  the  barricade  is 
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(a)  Hemispherical  Barricaded  and  Unbarricaded  (Bare)  Charges 

Figure  23.  Pressure  versus  Range  for  Various  Charge 
and  Barricade  Combinations 


t 
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(b)  Rectangular  and  Hemispherical  (Bare)  Urbarricaded  Chargee 

Figure  23.  (Continued) 
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higher  than  that  for  the  corresponding  unbarricaded  charge  while 
the  pressures  out  the  side  and  the  back  are  lower.  For  large 
ranges,  the  pressures  approach  those  of  the  unbarricaded  charge. 

The  pressure-distance  relations  for  a  bare  rectangular 
charge  with  the  same  weight  are  shown  in  Figure  23b.  When  com¬ 
pared  with  the  corresponding  result  for  a  hemispherical  charge, 
there  is  a  region  where  the  pressures  out  the  front  and  side 
of  the  charge  decrease  quite  rapidly.  An  explanation  for  this 
is  that  the  model  initially  assumes  that  the  pressure  waves 
propagate  in  directions  perpendicular  to  each  of  the  charge 
faces.  Thus,  rarefaction  waves  that  originate  at  the  corners 
of  the  charge  will  travel  parallel  to  the  faces  as  the  main 
pressure  wave  travels  out.  When  the  rarefaction  wave  reaches 
the  point  directly  out  from  the  center  of  a  charge  face,  a 
further  decrease  in  pressure  could  be  expected.  However,  this 
phenomenon  needs  more  study  and  as  more  data  become  available, 
the  appropriate  coefficients  in  the  computer  program  should  be 
changed. 

The  combined  charge  geometry  and  barricade  effects  are 
shown  in  Figure  23c  together  with  the  experimental  results  from 
"BIG  PAPA"  as  compiled  in  Figure  27  of  Reference  1.  The  correla¬ 
tion  between  the  "BIG  PAPA"  tests  and  the  predicted  pressure 

O 

distribution  out  the  back  of  the  barricade  (6  ■  180  )  is  not 
very  satisfactory.  However,  until  test  results  for  the  effects 
of  geometry  on  the  pressure  from  large  scale  charges  become 
available,  the  program  is  forced  to  use  the  data  from  the  limited 
range  of  charge  sizes  used  in  Reference  4. 
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Pressure  isobars  for  a  bars  hemispherical  barricaded  charge 
are  shown  in  Figure  24a.  The  isobars  for  the  unbarricaded  charge 
would  just  be  circles  and,  hence,  are  not  included. 

The  effect  of  charge  shape  geometry  is  dramatically  illus¬ 
trated  in  Figure  24b  which  gives  the  pressure  isobars  for  an 
unbarricaded  rectangular  charge  for  the  same  case  illustrated  in 

Figure  23b.  Similar  curves  for  the  rectangular  barricaded  charge 
are  shown  in  Figure  24c. 
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3. 


FRAGMENT  DISTRIBUTIONS 


To  illustrate  the  type  of  data  that  can  be  obtained  from 
this  portion  of  the  program,  the  bomb  parameters  given  in  Figure 
9. were  used  in  conjunction  with  the  stack  of  bombs  considered 
in  Part  2  of  this  Section.  In  addition,  the  following  infor¬ 
mation  was  assumed: 

a.  The  front  of  the  bomb  was  oriented  at  an  angle  of 

0 

270  counterclockwise  from  the  front  of  the  barricade 
(0B  «  270°) 

b.  Each  bomb  had  a  gross  weight  of  1500  lbs. 

c.  Each  bomb  contained  750  lbs.  of  TNT 

d.  There  were  333  bombs  in  the  stack,  and 

e.  The  number  of  effective  bombs  was  266. 

It  should  be  emphasized  that  the  bomb  characteristics  out¬ 
lined  above  are  purely  fictional  because  of  security  reasons  so 
that  it  will  not  be  possible  to  compare  results  with  experimental 
data-. 

Typical  fragment  trajectories  associated  with  these  bomb 
characteristics  are  shown  in  Figure  25.  This  particular  set 
of  curves  was  computed  using  a  10  degree  increment  for  beta. 

These  trajectories  indicate  that  fragments  ejected  by  the  bomb 

O  O 

for  approximately  5  <  3  <  50  will  impact  at  ranges  furthest 

from  the  bomb;  thus,  a  high  density  of  fragim  ts  is  to  be 
expected  at  these  impact  points.  However,  because  of  air 


friction,  the  impact  velocity  if  relatively  low. 

Impact  condition*  versus  range  at  selected  azimuth  angles 
ar^  shown  in  Figure  26.  These  curves  appear  as  piece-wise 
lines  of  constant  value  because  of  the  averaging  techniques 
used  in  the  program.  By  decreasing  the  increment  in  6,  the 
changes  would  not  be  nearly  as  abrupt  and  hence,  would  be  more 
realistic.  All  diagrams  indicate  a  high  fragment  density  at  large 
ranges. 

In  Figure  26a,  the  presence  of  the  large  impact  velocity 
is  due  to  the  absence  of  a  barricade  wall  at  that  azimuth  angle. 
Fragments  ejected  approximately  parallel  to  the  ground  strike 
the  ground  much  sooner  than  do  those  with  a  larger  departure 

l 

angle  0  and  hence,  these  large  impact  velocities  are  to  be 
expected. 

Figures  26b,  26c,  and  26d  show  the  variation  in  impact  con¬ 
ditions  around  the  stack  for  representative  values  of  the  szimuth 
angle.  For  similar  directions  of  propagation  from  the  barricade, 
the  impact  conditions  are  different  in  Figures  26e  and  26f  be¬ 
cause  of  the  different  fragmentation  properties  of  the  bomb  in 
these  directions. 


DISTANCE  FROM  BURST  POINT  (FT) 


(b)  AZIMUTH  ANOLC  ■  IOO* 

Figure  26.  Impact  Condition*  versus  Range  at  Selected  Azimuth  Angles 
for  Sample  Problem 
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AVERAGE  KUSS  OF  FRAGMENTS  (GRAMS)  AVERAGE  MASS  OF  FRAGMENTS  (GRAMS) 


4. 


APPARENT  CRATER  AND  EJECTA  DIMENSIONS 


The  apparent  radius  and  depth  of  craters  in  soil  according 
to  the  theory  of  Section  IV  and  the  computer  program  described 
in  Appendix  III  is  shown  in  Figure  27  for  a  range  of  10  to  500 
tons  of  TNT.  It  is  evident  from  the  figure  that  the  apparent 
crater  radius  and  depth  (70  ft.  and  21  ft. ,  respectively)  assoc¬ 
iated  with  a  100  ton  TNT  hemispherical  shot  at  the  Suffield  Ex¬ 
perimental  Station  have  been  used  as  the  reference  parameters  for 
this  program.  In  addition/  the  following  values  have  been  assumed: 

a.  The  dissipation  ratio  E®  *  0.3/ 

A 

b.  The  ejecta  parameter  6  •  3.1/  and 

c.  For  the  range  of  charge  sizes  10  to  500  tons  of  TNT, 

C  -  0.3  and  WL  «  500  tons. 

These  resultant  curves  of  Figure  27  show  significant  var¬ 
iations  from  scaling  laws  which  would  be  represented  by  straight- 
line,  relationships  on  this  plot  for  both  apparent  radius  and 
apparent  depth.  Because  of  the  attempt  to  include  fundamental 
quantities  in  the  theory ,  these  curves  should  predict  apparent 
crater  dimensions  more  accurately  than  the  curves  associated  with 
scaling  laws.  Of  more  importance/  the  parameters  and  C  can 
be  changed  if  more  data  warrant  such  an  adjustment.  Furthermore/ 
the  effect  of  charge  size  is  included  in  the  theory  and  changes 
in  shape  can  be  accounted  for  by  adjusting  W^. 

Strictly  speaking,  the  predicted  results  are  applicable 
only  for  the  earth  media  associated  with  the  reference  shot  where 
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Figure  27.  Apparent  Crater  Dimensions  in  Soil  versus  Charge  Height 


the  soil/  a  silty  clay,  had  a  weight  density  of  94  lb./ ft. 3 
However/  for  easy  reference/  results  from  shots  that  may  have 
been  in  a  slightly  different  soil  are  also  included  in  the  figure. 
It  is  assumed  that  bomb  stacks  can  be  replaced  by  a  bare  charge 
with  the  equivalent  amount  of  TNT  as  far  as  crater  dimensions  are 
concerned.  Thus/  some  of  the  “BIG  PAPA"  test  results  are  also 
included  in  Figure  27. 

The  depth  of  ejecta  can  be  obtained  from  Figure  28  where 
the  ejecta  depth  to  apparent  crater  depth  ratio  is  plotted  as 
a  function  of  the  ratio  of  distance  from  the  crater  center  (at 
the  surface)  to  the  apparent  radius  for  various  values  of  the 

A 

ejecta  parameter  8.  Hence ,  to  obtain  the  depth  for  a  particular 
location/  the  apparent  radius  and  depth  and  8  must  be  known. 

A 

When  the  earth  media  is  soil,  a  value  of  3.1  is  used  for  8. 

For  the  type  of  rock  encountered  in  the  Sailor  Hat  test, 
the  weight  density  was  190  lb. /ft. 3  A  hemispherical  charge 
of  500  tons  of  TNT  produced  a  crater  with  an  apparent  radius 
of  79  ft.  and  an  apparent  depth  of  38  ft.  Results  for  hemispher¬ 
ical  charges  of  other  sizes  are  shown  in  Figure  29  together  with 
the  predicted  values  of  8  and  the  dissipation  ratio. 
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Figure  2*.  Non-di,»n.ional  Ejecta  Shape,  a  a  a  Function  of  Earth  Media  Parameter 


Figure  29 »  Apparent  Crater  Dimensions  in  Basalt  Rock  versus  Charge  Weight 


5. 


SUMMARY 


Representative  sets  or  curves  have  been  given  based  on  the 
computer  programs  and  the  theory  of  the  previous  sections. 

For  given  charge  size  and  shape  and  barricade  dimensions,  the 
programs  yield  the  following  information: 

a.  Pressure  versus  range  and  impulse  versus  range  for 
various  azimuth  angles, 

b.  Pressure  isobars, 

c.  Impact  velocity,  number  of  fragments/ft.  ,  average 
fragment  mass  as  a  function  of  range  and  azimuth  angle, 

d.  Apparent  crater  radius  and  depth  as  a  function  cf 
charge  yield  and  shape,  and  earth  media,  and 

e.  Ejecta  depth  as  a  function  of  apparent  radius  and 
depth,  and  earth  media. 

Normally,  using  the  computer  programs  for  a  given  set 
of  input  data  will  be  the  most  convenient  method  for  design 
purposes.  However,  if  a  given  situation  occurs  repeatedly, 
then  it  would  be  more  convenient  to  construct  a  set  of  curves 
similar  to  those  illustrated  in  this  Section. 


SECTION  VI 


RECOMMENDED  INVESTIGATIONS 

1.  INTRODUCTION 

During  the  course  of  this  project  it  became  apparent  that 
a  large  amount  of  important  information  was  not  available.  On 
the  other  hand,  there  were  certain*  areas  such  as  cratering 
in  soil  for  which  a  great  deal  of  data  had  been  gathered.  This 
section  outlines  the  experimental  data  that  would  be  necessary 
to  corroborate  a  complete  theoretical  model  that  could  be  used 
with  some  degree  o.c  confidence.  The  requirements  for  data  have 
been  listed  in  the  same  order  as  the  topics  were  covsred  in  this 
report.  It  is  rather  obvious,  however,  that  more  than  one  type 
of  data  could  be  gathered  from  one  test. 

In  addition  tc  experimental  data,  the  overall  problesi  of 
safely  storing  bombs  suggests  a  corresponding  analytical  study 
in  linear  programming  where  such  factors  as  cost,  time  and  safety 
are  the  limiting  parameters.  It  is  believed  that  such  a  study 
would  be  extremely  useful  to  the  Air  Force  and  accordingly,  a 
brief  outline  of  the  approach  is  given. 


121 


2. 


PRESSURE  AND  IMPULSE  DATA 


a.  Single  Bombs 

Because  the  shape  of  a  bomb  is  significantly  different 
from  that  of  a  spherical  or  hemispherical  charge ,  peak  over¬ 
pressure  and  impulse  data  should  be  obtained  at  various 
angles  of  azimuth  and  distance  for  bombs  resting  on  the 
surface  of  the  earth.  Furthermore,  the  confining  effect 
of  the  bomb  casing  should  be  investigated  by  obtaining 
data  for  bombs  with  the  same  amount  of  charge  but  with 
different  case  thicknesses.  Intuitivelv,  one  might  expect 
the  initial  value  and  the  rate  of  decay  with  distance  of 
the  peak  overpressure  to  increase  as  the  casing  thickness 
is  increased.  However,  whether  or  not  this  is  true,  the  ex¬ 
tent  of  the  variation  should  be  investigated  both  analytical¬ 
ly  and  experimentally. 

b.  Bomb  Stacks 

Similar  pressure  and  impulse  measurements  should  be 
made  for  bomb  stacks.  Tests  should  be  conducted  with  the 
following  sequences:  (1)  Stacks  with  similar  shapes  but  with  different 
sizes,  and  (2)  Stacks  with  the  same  lumber  of  bombs  but  with  different 
shapes. 

Such  a  program  would  determine  whether  the  confining 
effect  of  several  bombs  tends  to  increase  the  initial  values 
of  the  peak  overpressure  above  that  expected  for  the  amount 
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of  explosive  present  or  whether  a  slight  time  difference 
in  detonation  of  the  individual  bombs  results  in  a  lowering 
of  the  expected  value.  Obtaining  values  for  the  unit  im¬ 
pulse  is  aiso  extremely  important. 

A  side  benefit  of  these  results  might  be  the  preference 
for  a  particular  stack  shape  and  size  based  on  pressure 
and  impulse  limitations  rather  than  stacking  convenience, 
c.  Barricades 

Very  little  pressure  and  impulse  data  appear  to  be 
available  for  the  region  immediately  outside  the  barricade. 
For  a  given  stack  size  and  geometry,  a  series  of  tests 
should  be  conducted  for  barricades  with  various  dimensions. 
Significant  differences  could  be  expected  close  to  the 
barricade  but  for  regions  farther  away,  the  results 
should  approach  those  of  the  unbarricaded  stack. 
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3. 


FRAGMENT  DATA 


There  appears  to  be  sufficient  data  on  the  fragnent  sizes, 
number  of  fragments  and  initial  velocities  of  fragments  for  indi¬ 
vidual  bombs.  One  aspect  that  could  be  handled  simultaneously  with 
the  pressure  data  of  the  previous  section  is  the  acquisition  of 
information  concerning  the  distribution  and  impact  velocities 
of  fragments  from  bomb  stacks  of  various  sizes  and  shapes.  These 
data  would  yield  essential  information  concerning  interaction 
effects  and  hence,  whether  or  not  the  simplified  theory  of 
Section  III  is  adequate.  If  the  theory  associated  with  the 
fragment  trajectories  is  fairly  accurate,  then  the  effect  of 
barricades  can  be  predicted  quite  confidently  and  hence  a  special 
experimental  program  considering  barricades  is  not  warranted  in 
this  connection. 
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CRATER  DATA 


A  vast  number  of  cratering  programs  have  been  conducted 
fpr  bare  hemispherical  and  spherical  charges  but  predicting 
crater  dimensions  with  the  use  of  scaling  laws  is  not  com¬ 
pletely  satisfactory.  Examples  of  unexpected  crater  sizes 
include  500  ton  shots  of  the  Suffield  Experimental  Station 
and  of  "Operation  Sailor  Hat".  It  is  believed  that  the  level 
of  the  underground  water  table  may  have  affected  the  results 
of  the  first  shot  mentioned  above.  The  variation  from  the 
normal  crater  shape  of  the  latter  shot  may  be  due  to  the  type 
of  rock  at  that  particular  location. 

In  light  of  the  large  amount  of  crater  data  available 
it  does  not  seem  advisable  at  this  point  to  conduct  more  tests 
of  the  same  type  until  there  exists  a  better  understanding  of 
the  effects  of  the  various  parameters  that  describe  the  earth 
media.  However,  in  connection  with  this  program,  a  number  of 
tests  do  seem  advisable. 

One  of  these  is  concerned  with  the  effect  of  charge  shape 
on  the  crater  dimensions.  For  example,  for  the  same  earth 
media  and  the  same  charge  weight,  a  series  of  shots  should 
be  conducted  in  which  the  charge  shape  is  varied.  Typical 
examples  would  be  rectangular,  cubical,  and  triangular  shapes. 

Another  series  could  involve  the  use  of  stacked  bombs, 
again  for  one  earth  media.  Such  tests  would  illustrate  the 
interaction  effects  of  bombs  and  the  effect  of  stack  size  and 
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shape.  Furthermore/  tests  of  this  kind  could  lead  to  comparisons 
with  the  results  of  bare  charges  of  the  type  mentioned  above. 
According  to  the  theory  developed  in  Section  XV  the  shape  with 
the  lowest  center  of  mass  will  produce  the  largest  crater  and 
the  basic  hypothesis  behind  this  theory  should  be  checked. 

Al*o  of  considerable  interest  is  the  effect  of  barricade 
size  on  crater  dimensions.  Experimental  data  in  this  area 
would  be  very  useful  for  both  theoretical  and  immediate  prac¬ 
tical  use. 
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5.  OPTIMIZATION  OF  STORAGE  AREAS 

In  deve loping  a  munitions  storage  area,  several  factors 
must  be  taken  into  consideration.  These  include  the  cost  and 
availability  of  land,  possible  methods  of  stacking  bombs  and 
building  barricades,  degree  of  safety  in  connection  with  frag¬ 
ments,  cratering,  peak  overpressure  and  impulse,  the  time 
available  to  stack  the  bombs,  and  so  on.  With  the  experimental 
and  theoretical  methods  available  it  appears  that  it  might  be 
possible  to  develop  a  computer  program  that  would  optimize  a 
given  parameter  under  a  given  set  of  circumstances. 

One  example  could  be  the  following!  Suppose  that  a 
particular  amount  of  munitions  had  to  be  stored  on  a  given  area 
and  the  safety  requirement  was  primarily  one  of  ensuring  that  no 
fragments  landed  outside  this  area.  The  problem  would  then  entail 
finding  the  appropriate  combination  of  munition  stack  shapes, 
dimensions  and  distributions  of  the  stacks  together  with  barricade 
shapes  and  dimensions  that  would  satisfy  this  requirement.  If 
more  than  one  combination  was  adequate,  then  additional  factors 
such  as  cost  and  time  could  be  included  in  the  program. 

Another  possibility  would  be  to  determine  the  "safest" 
possible  arrangement  that  could  be  developed  in  a  given  amount 
of  time.  To  determine  the  safety  aspect,  degrees  of  importance 
would  have  to  be  attached  to  each  of  the  hazardous  factors 
associated  with  a  munitions  dump.  Likewise,  the  parameters 
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describing  the  time  to  stack  the  bomba  would  be  included  and 
by  using  the  principles  of  linear  programming,  the  optimum 
configuration  would  be  predicted* 

Since  there  are  a  large  number  of  possibilities  that 
could  be  of  interest  to  the  Air  Force,  it  would  seem  that  the 
possibility  of  developing  a  program  that  could  predict  optimal 
arrangements  should  be  considered.  The  results  could  be 
significant  improvements  in  safety,  time  and  cost. 


SECTION  VII 


CONCLUSIONS 

■  Analytical  models  and  subsequent  computer  code*  have  been 
developed  for  typical  large  quantity  high -explosive  detonations 
of  those  types  of  conventional  Munitions  stored  by  the  Air  Force 
in  aboveground  barricaded  Modules.  The  paraaeters  that  are  pre¬ 
dicted  from  these  codes  include  peak  overpressure ,  unit  iapulse, 
distributions  of  fragnent  iapact  coordinates  and  velocities , 
crater  dimensions  and  depth  of  ejecta.  The  geometries  of  the 
boab  stack  and  the  barricade  are  taken  into  account  as  well  as 
the  type  of  earth  Media.  Results  of  large-scale  tests  had  in¬ 
dicated  that  the  burning  time  of  detonation  and  the  point  of 
initiation  of  a  bomb  stack  were  not  too  significant  and.  hence , 
these  parameters  are  not  considered. 

Because  of  the  lack  of  vide -scale/  definitive  experimental 
data,  it  was  considered  More  appropriate  to  use  a  combination 
of  analysis  and  empirical  curve  fitting  in  connection  with  the 
models.  Accordingly,  some  engineering  judgement  must  be  used 
with  the  computer  codes  since  there  nay  be  some  disparity  in  the 
values  of  parameters  for  an  actual  problem  and  the  values  that 
are  used  in  the  model. 

Evo.ry  effort  was  made  to  use  the  latest  experimental  and 
theoretical  results.  However,  there  are  several  areas  in  which 
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vital  gaps  in  fundanental  knowledge  exist* 
of  testing  is  necessary  to  obtain  this  information  and  the  types 
of  tests  that  should  be  run  are  listed  in  this  report. 

.  Of  more  immediate  benefit  to  the  Air  Foroe  is  the  development 
of  computer  codes  which  would  optimize  a  base  layout  under  a  given 
set  of  conditions.  Such  codes  could  produce  the  safest  possible 
combination  of  munition  modules  for  a  given  area  or,  alternatively 
the  cost  and  amount  of  time  required  to  construct  a  safe 
munitions  storage  area  could  be  minimized. 

The  codes  that  have  been  developed  should  be  extremely  help¬ 
ful  in  predicting  the  danger  zone  for  a  munitions  storage  area. 

In  selecting  models  that  yield  governing  parameters  with  respect 
to  the  safety  of  a  base,  a  conservative  approach  has  been  adopted. 
Thus,  as  better  technical  information  becomes  available,  these 
codes  can  be  adjusted  and  better  use  of  available  ground  apace 
for  storing  munitions  can  be  made. 
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APPENDIX  I 
BLAST  EFFECT  PROGRAM 

1.  •  FORTRAN  PROGRAM  DESCRIPTION 

The  computer  code  described  in  this  Appendix  follows  the 
development  of  Section  II  on  Blest  Effects-  The  computations 
are  performed  in  the  same  order  as  outlined  in  the  development 
of  the  analytical  motel  whenever  possible-  A  flow  chart  for 
the  computer  code  is  included  in  Appendix  1-2-  A  printout  of 
the  computer  code  is  included  in  Appendix  1-4- 

The  first  section  of  the  -computer  code  establishes  (a) 
the  coefficients  for  the  various  polynomials  needed  during  the 
computations ,  (b)  overpressures  to  be  solved  in  isobar  option, 
and  (c)  degrees  of  different  polynomials  require.!  in  computa¬ 
tions,  and  (d)  lists  format  statements - 

Each  problem  to  be  run  requires  a  control  card  to  be  read- 
If  wove'  data  are  to  he  read,  the  control  card  indicates  the  samnt 
to  be  read  and  the  variables  that  are  to  be  reassigned-  The 
control  card  also  sets  the  parameters  describing  bomb  stack  and 
barricade  and  indicates  azimuth  angles  for  which  the  pressure 
and  impulse  are  desired-  After  the  control  card  is  read,  the 
heading  on  the  output  is  printed  to  record  the  parameters  of 
the  problem  for  future  reference - 

The  solution  is  obtained  by  computing  the  pressure  and 
positive  impulse  for  a  bare  hemispherical  charge  of  TNT  at  the 
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scaled  distances  Z  -  4.0,  4.1,  4.2,  ...,  480,  490,  500.  The 
incremental  steps  of  Z  are  as  follows » 

AZ  ■  0.10  for  4.0<Z<5.0 

AZ  »  0.20  for  5.0<Z<10.0 

AZ  *  1.0  for  10.0<»<50.0 

AZ  -  2.0  for  50. 0£Z< 100.0 

AZ  *  10.0  for  100.0<l<500.0 

Pressure  and  impulse  are  evaluated  from  the  polynomials  pre¬ 
viously  fitted  to  the  In  P-ln  z  and  In  I-ln  z  representations 
ot  Figure  8.  If  a  scaled  distance  Z  falls  within  the  confine 
of  the  barricade  the  pressure  and  impulse  are  not  computed. 

This  evaluation  of  pressure  and  impulse  for  a  bare  hemispherical 
charge  is  completed  previous  to  statement  number  335  in  the  code. 

If  the  charge  of  TNT  is  rectangular  rather  than  hmd- 
spherical,  the  effects  of  the  change  in  geometry  are 

evaluated  according  to  the  theory  of  Section  II-4.  The  evalua¬ 
tion  of  effect  of  charge  geometry  on  pressure  and  impulse  dis¬ 
tributions  is  performed  between  statements  335  and  460  in  the 
computer  code.  The  procedure  used  is  to  evaluate  the  pressure 
ratio  and  positive  impulse  ratio  perpendicular  to  the  faces  of 
the  rectangle  and  along  a  line  through  the  center  of  mass  of 
the  charge  and  the  corner  where  the  faces  of  the  rectangle  meet. 
The  values  of  pressure  ratio  and  positive  impulse  ratio  for 
angles  other  than  those  perpendicular  to  the  faces  and  through 
the  corners  are  ob  ained  by  linear  interpolation  between  the 
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two  known  valuta  that  span  tha  angle.  La  tar  in  tha  program  tha 
prasaura  ratio  and  poaitiva  impulse  ratio  will  ba  uaad  to  multi¬ 
ply  tha  corraaponding  valuta  of  praaaura  and  poaitiva  impulaa  at 
a  diat&nce  2  for  a  bara  hamiapherical  charga  to  obtain  praaaura 
and  poaitiva  impulaa  for  tha  ractangular  charga  baing  mode  lad. 

If  a  barricade  ia  praaant  around  a  atack  of  TNT,  4 to  ahapa 
ia  aaaumed  to  ba  ractangular  with  tha  cantar  of  masa  of  tha 
atack  of  TNT  vary  naarly  corraaponding  to  tha  gaomatric  cantar 
of  tha  barricada  in  tha  plan  viaw.  Tha  walla  of  tha  barrioada 
art  attuned  to  ba  parallel  to  tha  aidaa  of  tha  atack  of  TNT. 

Tha  affact  of  tha  praaanca  of  a  barricada  ia  evaluated 
according  to  tha  theory  of  Section  II-5  between  atataaant 
numbers  465  and  555  in  tha  computer  code.  Tha  praaaura  and 
poaitivo  impulaa  ration  (BP  and  BI) ,  that  will  ba  multiplied  by 
tha  prep  jure  and  poaitiva  impulaa  for  a  bara  hamiapherical  charge 
at  the  corraaponding  diatancea,  can  be  evaluated  directly  for 
directiona  perpendicular  to  tha  walla  of  tha  barricada  and  out 
tha  rear  comer  of  tha  barricade  by  interpolation  from  tha 
input  data.  Tha  ration  for  anglea  not  equal  to  thoae  juat  men¬ 
tioned  are  obtained  by  a  linear  interpolation  acheme  from  tha 
valuea  obtained  for  direction  perpendicular  to  tha  walla  and  out 
tha  rear  corner  of  the  barricade.  It  ia  only  necaaaary  to  deter¬ 
mine  which  two  known  directions  the  angle  theta  liea  between  to  choose  tha 
proper  data  for  interpolation. 

This  concludes  the  evaluation  of  the  pressure  and  positive 
impulse  ratios  that  are  needed  to  estimate  the  effective  change 
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in  pressure  and  impulaa  diatribution  produced  by  tha  change  in 
oharga  gaomatry  from  a  hamiapharioal  to  a  ractangular  ahapa 
and  tha  praaanca  of  a  barricada  around  tha  oharga. 

To  obtain  tha  aatimatad  praaaura  and  poaitiva  impulaa 
distributions ,  it  ia  only  nacaaaary  to  multiply  tha  praaaura 
and  poaitiva  impulaa  diatributiona  for  tha  bars  hamiapharioal 
oharga  by  tha  appropriata  praaaura  and  impulaa  ratioa  aa 
darivad  abova.  Thaaa  computationa  ara  parformad  and  printad 
out  in  tha  cocnputar  coda  between  atatamanta  560  and  565. 

If  iaobar  plots  ara  daairad  to  indicata  diatanca  to  linaa 
of  oonatant  praaaura  tha  atatamanta  630  through  680  ara  axacutad 
by  aatting  I  BO  -  0,  otherwise  I  BO  -  1.  Tha  diatancaa  to  pointa 
of  aqual  praaaura  along  linaa  aaparatad  by  5°  incramants  ara 
avaluatad  by  starting  at  tha  axtrama  distance  (Z  -  500)  and 
comparing  expected  pressure  to  isobar  pressure  for  decreasing 
valuss  of  S  until  expected  pressure  is  greater  than  tha  isobar 
pressure.  When  tha  points  on  either  side  of  the  distance  where 
isobar  praa«afw  quels  expected  pressure  ara  found,  tha  distance 
to  tha  isobar  pressure  is  found  by  linear  interpolation.  If 
tha  expected  pressure  along  a  direction  of  propagation  does 
not  attain  the  isobar  pressure,  tha  symbol  PNR  is  printad  in¬ 
stead  of  the  distance.  This  indicates  that  the  isobar  pressure 
is  not  reached  along  that  direction.  This  method  of  evaluation 
continues  until  all  the  distances  to  the  isobar  pressures  have 
been  evaluated  along  the  directions  indicated. 
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Mshrio  *  m 

READ  19  • 

— ~~~wWTr~% 

H8  lF_'(NtMET 

ms  nnn  *  i8i 
_ jn  *  t_  _* 

CO  TO  ITS 
270  WWW  «  NTH 

JN  *  NTMETA  ♦  1 

C  " _ 

C  _  MAJOR  LOON-  INP 
C  NHEN  NTMETA  IS 
C  TMETA  IS  TO M 


nmwnTfmil 


lilf  NT|f:iA_«A0  IN  XJ  -1 
ONLY  FOR  ANCLE  NTMETA 
[Oe«ECS  SET  KEEN  t  A  NO  ltt 


8tB%Sl 

Vofi%2l' 


089%?% 

908%26 


980%33 

008436 


880%%1 

989%%2 

008%%% 

000%%7 

000%88 

000%59 

900%97 

000%60 


000%6% 

000%70 

008%7t 

000%71 

990%7? 

900%72 

000%73 

009%73 


TO  NEOUCE  THE  ANOUNT  Of  PRINTOUT  WHEN  SOL  WING  PON  PRESSURE  ISOOARS 
SUPPRESS  PRINT  OJLfRESSVRES.  MQ  IMPULSES  VS.fiXST.POR  ANCLES 
COMPUTED  BY  NARIR51IM1S*  21 1x119*2 St  *963*  MHO  C21  INTO 


NTMETA  *  NN  -  1 


SET  FIRST  INCREMENT  OF  01STANCE 


t  *3.9 

288  PRINT  18t 


NTMETA 


PRINT  OUTPUT  COLUMN  HEAOINCS 
289  PRINT  182 

298  PRINT  US  “  "I _ ~ 


EfflIE 


M  *  11  _ 

S  *  8.18  _ 

9PMIl*ATAN?f RL9 80)  •  RA0_ 
SPMI2»ATAN2<RLt-90)»  RAO 
NUN  *  8" 

RR  *  NTMETA 

THETA  *  PR  t  RAO 

NOEG  *  NTMETA  7  9_4  $ _ 1 

SECONO  MAJOR  LOOP  INOE 

380  00  569  «*1,M 


nn 


NUM  *  MUM  ♦  1 
9P  *  8.0 
FG  *  0.9 
G!  *  0.0 
91  *  0.0 
TOP (MUM)  * 
9S0PI (NUN) 


9.0 
*  9.0 


•••47% 

•01%77 

miii 

m««4 


ITftfit 

lltfiS 

rmir 

ntf» 

wmrr- 

•••939 

•••939 


POPfNMCtNON)  »  •  •• 
7307! U026,NUH>  •  1.1 

2*  7  ♦  S 

2*fN6C6|NU'n  ■  i 
2ft  *  i  •  cos  tfHEfit 
Td~*  WSIfC! 


“  h  swfwmr - z  T 

trcimi-  NtflEY  j>F~Sl9~,1iS,3i5 
3U  t7  lii.iti.m 


ifrir  nc 


njito  — 


e 

c 


3*1 1  •  Ktdl(2f 
ft*  1.1 


MfT  7*ESSt*F  ¥T 


••1937 

1! —  - 

00  329  J  *  ltNOOP 

1119%9 

N  *  J  -T 

ST09%6 

*29  C  •  C  ♦~WIJ1  •  W  *»Ni 

110993 

OOPfNUXr  *  ExPtftf* 

101996 

0  *T.l 

•11997 

•••999 

IMHI 


ft* 

ft 


m«*  it#uHTW~r 


00  331  J  *  IVMOSOPf 
N  *  J  -  1 


111977 

010692 

011617 

11161% 

•01621 

•00622 

000623 

00067% 

00062% 

000629 

000629 

000626 


ft 

C 

C 

ft 

ft 


331  0  *  0  ♦  SOPtfJI  **  Fx  76  Nt 

IfOTTTWWr**  tkPIS)  »Wfl 


OECtWNlNC  07  ftlONtfiV  OPTION  -  17  0CSIN70»THIS  SECTION  SOLVES  TON 
TN0  74CT0NS  WTNO  «JWIChTnE  TO  9C  NOLTIPLIEO  IT  SANE  PtESSIPE 
•no  nine  Impulse  Include  nectancjlon  chance  ceonetny  eftect 


'  ITTNCfl  -  II  331,61mm - 

339  171  "2  -  91.11  3%0,%FOY%6a - 

3%1  Ml,  *  1.1/C  «0/3»D  TF1  ♦  1  ¥'(SB/Sll  1 
•NO  •  1.1/1  (SL/ftff  *  f.O  ♦  2.0MSL/SOH 
•1  «  1.0'  “ 

*2  •1.1 - * - 

•3  *  l.fl 


r 


91 

92 

93 


0.0 

0.0 

0.0 


ft 

ft 

c 

c« 

ft 


PNE  •  1.0 
SI»P  *  0.0 


evaluate  pocv?w*TOSHrry  'at,  n,  it  2 


000627 

0006%n 

ft 

no  L00  L  * 
J  *  L  -  1 

l.NOGPO 

0006%l 

9906%6 

•1  *  A(l.L) 
•2  «  A (2, L) 

•17  •*  J>*  A1 
*12  77  JYTTS  ~ 

000693 

%00 

r 

•9  «  A<3,1» 

•tZ  •*  J>*  »3 

000662 

I712.LE.26. 

0>  GO  TO  %09 
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ITffjmtW  03  MOT 
A  Of  THE  CHARGE*  CHOOSE 


TtAFOLATE  FOII  THE  PRESS ^  AMO 


SiAtn  ♦  Sial 


^ffirWflOM  SbtVES  FOE 
txnito  OT  SAME  PFESSUOE 


irr  octkEtw  ggri  aggnm  the  appropriate  values  of 
AIARtcAPE  WtfSS.HfcPfiCT  (»I)  to  EVALUATE  THE  POLYNOMIAL 


♦_fFE«i,J) 


iNTEMFoL  nr  aFTVftiriw euRVKFo*  t«e  appropriate  values  of 

OAWtTTWr  fflRPlSE  EFFECt  CStEI  TO  EVALUATE  The  FOLYMOMI AL 
FOA  TflWWTtTABr  THFUWE  WtO  (HI). 


00  SOS  J  ■  UMflE 


011136 

111HT 

011190 

001193 


El  *  CTR>*  (OtEtE.JT  -  ♦  9IE(i*JI 

919  91  •  «I  ♦  El  •  (2  *•  l) 

- TO  INSURE  AEASOWIBUE  RESULTS  GUARANTEE  THAt  ?6Rl  6AARICA0E0 

CHARCT  THE  RE  SUITS  APPROACH  THOSE  FOA  A*  UHOAAATCAOEO 


IF I  7.LE.4.5)  910*915 
910  T9I  ■  FT 
TIP  ■  FF 


011171 

0113"? 

011213 

091323 


mm 


BO  Ilf  2 

mwi 


60  TO  9*0 


r 

C- 

0‘ 


FOLLOWING  STATEMENTS  TO  595  APPLY  F OR  THETA  CHEATER  THAN  PI/2 


920  8THETA  «  ATANIBD/BL)  *  RAO  ♦  90. 

- TB  «  OOTkgTl  S - 5 - 2SJ- 


e 

c 

c 

TT 


CHECK  L PC AT IPX  OF  THETA  WITH  RESPECT  TO  AH6LE  THROUGH  CORNER 

_ OP  BARRICADE  TO  QCTCRHINE  APPROPRIATE  SET  OF  EQUATIONS 

__ _ TO  INTERPOLATE  FOR  EFFECT  OP  RARRICAPE  ON  PRESS.  ANO 

_ IHPtX.SC  RATIOS  (  BP  ANO  II  » _ _ 


0012*2 

0912*5 

101163 

00126* 

Hilt* 

001271 

001277 

001313 

00131* 

nnrr 

101329 

00i33t 

001333 


_IP  JTB  8THETA)  529.925.9*0 


525  OO  530  J 
l«  J  -  1 


l.NOBPE 


TN  •  NTHETA  -  98 _ 

-Pg  ■  < <TR)  'JBTjglA^AO  >M9PgJ3i JL-yE (2, J) ) »9PE (2,U) 


c 


530  8P_*_BP  ♦  PE  »  C2  ♦»  L> 


oo  s35_  J?  l.NOBIE 
L«  J  -  1 


■ .  II  9  C<TN1/(BTHETA  -  51.)  )*(BIE(3.  J)»BIE(2.J) )  ♦  BIEI2.JI 
535  BI  «  81  ♦  El  •  12  ••  L» 


001*31 


001*33 


mro- 

001352 

001353 

001356 


KHW 
0 01*01 
001*02 
001*05 

001*13 

nuff 

001*21 

001*22 

001*23 


ipiSlclalo?  bp  «  i!a 

60  TO  560 


5*0  OO  5*5  J«  1,  NQBpT 
L*  J  -  1 


PE  *  (JTB  -  BTHET M/JJ 8 0  .-OTMET  AJ)»  IBPE^.Jt  -  BPE(3,J>>  ♦  B»E( 

13,  J)  _ _ _ 

5*5  IF  ■  BP^  ♦  PE  »  f  2  »»  LI _ _ 


"BS  5«~  J~«  l.NftBfi 


L  *  J  -  1 


El  -  ((TB  -  8THETA)/(lt0.  -  5TH£T«})  VJBTFU,  J) -BIE (3„j)  )40IE(3aJ) 
550  BI  *_BI  J  El  •  <2  U _  _ 


tP<8I.6E,l»0J  ?I  »  Iti. 

IF (BP. 6E. 1.0)  *F  «  1,0. 


60  TO  560 

555  BI  *  1.0 
BP  «  1.0 


C _ 

C - HULTfPLT  OVER  PRESSURE  FOR  A  _8ARC_  HEMISPHERICAL  CHARGE  BY  EFFECTS 

C  OF  rcONETRT  (FG)  ANO  BARRICADE  (BP>  TO  OBTAIN  THE 

C  PR  CTED  VALUE  FOR  PRESSURE  AT  OISTANCE  2  AND  ANCLE 

C  THl *  A 

C 

560  F0P(N0F6»HUH)  »  BOP  (HUH)  *  BP  »  FG _ _ _ _ _ 

_  _ 

C - HULTIPLT  INPULSE  FOR  A  BARE  HENISPHERICAL  CHARGE  BT  EFFECTS 

C  OF  GEOMETRY  (Cl)  ANO  BARRICADE  (BI)  TO  OBTAIN  THE 

C  PREOICTEO  VALUE  FOR  IMPULSE  AT  OISTANCE  Z  ANO  ANGLE 

C  THETA 

r 


C 

r- 

c 


FSOPI  (N0f6,'NUH)  ■  OSOPt(NUH)  •  BI  *  61 
■  CONVERT  2  TO  FEET 
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HP-«* 


011436  RD(NOEG,NUM>  •  2R(NDCG, MUNI  •  C*»Q 

001440  963  PRINT  101,  ZR(N0EG,NUN)  ,flOP(NUN)  ,E0P(N0EG,  NUN)  ,9S0PI(NUN)  ,FSOPI(NO 

1EG.NUH)  ,RO(MOEC,NUm 
001473  969  CONTINUE 

00T676  IF  (7  -  4.0)  970,979,979 

101900  970  N  *  40 

wmri - 8*  S“I7m - 

011903  GO  TO  Jo 0 _ '  _ _ 

011903  979  IF  CZ  -  9.9»  900,919,905 

001906  910  H ~m  09 

001907  ‘  S  £  67f9 - 

OOWir  GOToHO" 

001911 - nrir  <2  -  44.41  CTnnnrwv . .  ~ 

001914  990  H  «  40  '  ' 

001919" -  -  S  *  1.0 - 

001917  '  GO  TO  300 

001917  600  IF  72  -  99.n“6f9^lO,Ti0  - 

001927  609  N  «  25  -  - 

001923 - 9  STTTT - 

001929  GO  TO  300 

001929  610  IFlZ  -  094.01  6l9,7?i;6?0~  " 

001930"  "619  H  *~40 - - 

001931  S  «  ll.3~  -  - 

001933  ~  GO  TO  300  “  - 

191933 - 121  POINT  123 - - - 

001936  "629 " CONTINUE”  “ 

C  t  SOB  AO  COORDINATE  OPTION  fHtS  SECTION  SOLVES  EON  THE  7  AT  THE 
c  "ANGCerO  -  100  FoffHr~WS§U#fS  INPUT  AS  P(N1  IN  THE  DATA 

c  '  " 


001961 

001943 

001949 

001947 

001594 

C 

630 

TE  TIBfl  -  TT  6TT,H3,ni9 

00  60f  N  *  1,72 

PP  *  P(Nl 

PRINT  199  ,  PP 

P1TNT  199 

— 

u01*>57 

PRINT  156 

.  _ _ 

.  . 

001 56? 

00  600  J  »  1,37  " 

. 

-  -  - 

001964 

NN  .*  (J  -  11  *  5 

001966 

LL  «  142  —  -  ~ 

001967 “ 

DO  670  L  *  1  ,  iVi 

0ol67i 

XL  *  It  -  l 

001972 

IE  (  F0P(J,LLl  -  PP  ) 

670 

,  660  ,  645 

001607 

645 

LI.L  *  LL  ♦  t 

001610 

72  »  ZR ( J,LLLI 

001613 

71  «  ZR(J»LL) 

001615 

P2  »  FOP ( J»LLL) 

0016l6 

pi  *  tup  f  jvm 

001620 

71 ( J,N)  *  (122-711  * 

(Pi  - 

PP  l»  /  (Pl-I 

001676 

GO  TO  675 

0' 1 627 

660 

21  ( J,N)  *  7RCJ,U 

6.. .636 

GO  TO  675 

001637 

670 

CONTINUE 

001641 

PRINT  138,  W  - 

—  . 

*  "  1  " 

001644 

GO  TO  600 

001647 

675 

PRINT  157  ,  HN  ,  ZI(J 

♦  Nl 

001661 

640 

CONTINUE 

c 


C  TO  ®UN  A  SINGLE  PROOLEH  TNE  FOLLOW  INC  THO  CARDS  SHOULD  OF 

C  645  CALL  EXIT  - - -  ‘  '** 

«*  END 

r  TO  PUN  A  SERIES  OE  PROBLEHS  TH*  FOLLOWING  THREE  CAROS  SHOULD  BE 
r,  545  GO  TO  199 
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5. 


INPOT  FORMAT 


The  input  to  the  Blest  Program  normally  is  one  card  (See 
Figure  30)  consisting  of  the  following  parameters  and  control 
options  in  their  respective  order  on  the  input  card: 

a.  '  SL  -  F6.2  Format — stack  length  in  feet,  measured 
parallel  to  the  open  end  of  the  barricade,  or  perpendicular  to 
the  line  for  theta  (d)  equal  to  zero  degrees  (See  Figure  31). 

b.  SD  -  F6.2  Format — stack  depth  in  feet,  measured  per¬ 
pendicular  to  the  open  end  of  the  barricade  or  parallel  to  the 
line  for  theta  (6)  equal  to  zero  degrees  (See  Figure  31). 

c.  SH  -  Ffil2  Format — stack  height  in  feet. 

d.  BL  -  F6.2  Format — barricade  length  in  feet,  measured 
parallel  to  the  open  end  of  the  barricade  (See  Figure  31). 

e.  BD  -  F6.2  Format — barricade  depth  in  feet,  measured 
perpendicular  to  the  barricade  open  end  (Sea  Figure  31). 

f.  QUAN  -  F8.0  Format-- quantity  of  explosives  in  explosive 
stack  in  equivalent  pounds  of  TNT.  The  value  of  QUAN  is  to  be 
determined  by  converting  the  weight  of  explosives  in  the  stack 

to  its  equivalent  weight  of  TNT  and  multiplying  by  the  appropriate 
bomb  factor. 

g.  N THETA  -  14  Format — for  0  <_  NTHETA  <  180  the  computer 
code  will  solve  for  peak  overpressures  and  positive  impulses  at 
selected  values  of  scaled  distanoe  Z  (4.0  <  Z  <  500.0)  only  along 
the  angle  specified  (See  Figure  31) .  For  NTHETA  ■  -1  the  com¬ 
puter  code  will  solve  for  peak  overpressure  and  positive  impulse 


i 


•t  -Lctad  valua.  of  diatanoa  I  (4.0  <  ,  <  ,oo.„>  .ion, 

Una.  from  that,  a  aqual  to  aaro  dagraa.  through  thata  a  aqual 
to  180  degrees  in  increaents  of  5  dcfnu 

h.  MOPD  -  14  format— n«w  gantry  pra.aur.  d.t.  optlon. 
for  hgpd  -  o  tha  computar  coda  uaaa  intarnal  polynomial 
coaffioianta  to  account  for  ataok  gaomatry  affacta  on  paak  ovar- 
praaaura  along  lina.  parpandicular  to  th.  vartlcal  atach  f.o.a 
.nd  through  th.  .tack  can tar.  for  NGPD  >  0  th.  computar  cod.  will 

r*.d  “W  P°l3'K~i*1  to  aooount  for  gaomatric  a,f.ct.. 

MGPD  muat  to  m  to  a  valua  aqual  to  th.  higha.t  dagra.  of  th. 
naw  polynomiala  whoa,  co.ffici.nt.  «.  to  b.  ra.,.ig».d.  th... 

co.ffici.nt.  will  b.  for  Polynesia  oorraapondin,  to  c,  and 

D  in  Equation  (7)  or  *1,  *2,  and  A3  in  tha  ooaputar  coda.  Tha 

raathod  for  obtaining  tha  appropri.t.  valua.  for  th.M  ooaffiolant. 
i*  -ascribed  in  Section  II. 

To  input  now  pra.aur.  data  for  th.  appropriat.  polynomial 
coafficiant.  a  total  of  NGPD  ♦  1  data  card.  ar.  rm,uirad  for 
polynonialc  of  maximum  ordar  aqual  to  NGPD.  |.0h  data  card  will 
contain  on.  coafficiant  for  aach  polynomial  c,  and  D  in  a  3*14.7 
Format.  Th.  fir.t  card  will  contain  tha  aaroth  ordar  oo.ffioi.nt. 
*or  B,  c,  and  D  raapactivaly.  Th.  a.oond  oard  will  contain  th. 
fitat  ordar  coafficiant.  for  B,  C.  and  D  raapactivaly,  ate.  Th. 
final  card  in  thi.  data  act  will  contain  th.  NOPD**  ordar  co- 


efficients  for  B,  C,  and  D  respectively 
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‘  1  3! 


t  r  .  • 


' 


i.  N6ID  -  14  Format — new  geometry  impulse  data  option. 

* 

For  N6ID  ■  0  the  computer  coda  uses  intarnal  polynomial 

> 

'  i  t  t 

coaffioiants  to  account  for  atack  gaoraatry  affact  on  impulse 
along  a  lina  parpandicular  to  the  vertical  stack  faces  and  through 
the  stack  canter.  For  NGID  >  0  the  computer  code  will  read  new 
polynomial  coefficients  to  account  for  geometric  effects.  The 
number  of  coefficients,  number  of  data  cards,  and  the  values  for 
the  coefficients  can  be  determined  as  in  h  above.  The  variations 
from  part  h  is  that  the  coefficients  will  be  for  polynomials  F, 

G,  and  H  in  Equation  (7)  or  Bl,  B2,  B3  in  the  computer  code. 

j.  NBPD  -  14  Format-new  barricade  pressure  data  option. 

For  NBPD  ■  0  the  computer  code  uses  internal  polynomial  co¬ 
efficients  to  account  for  barricade  geometry  effects  on  peak 
overpressure  along  lines  perpendicular  to  the  barricade  walls 
and  through  the  barricade  center,  and  along  lines  through  the 
comer  where  the  barricade  walls  meet  and  the  barricade  center. 

For  NBPD  >  0  the  computer  code  will  read  new  polynomial  coefficients 
as  in  h  above.  These  new  coefficients  will  be  for  polynomials  to 
define  pressure  ratios  in  directions  A,  B,  C,  and  D  in  Figure  8. 

The  number  of  coefficients,  number  of  data  cards,  and  values  of 
coefficients  can  be  determined  in  a  manner  similar  to  part  h 
above.  Each  data  card  will  contain  one  coefficient  for  each 
polynomial.  The  data  card  Format  will  be  4E14.7. 

k.  NBID  -  14  Fortnat— New  barricade  impulse  data  option. 

For  NBID  ■  0  the  computer  code  uses  internal  polynomial  co¬ 
efficients  to  account  for  barricade  geometry  effects  on  impulse 


ratios  in  directions  A,  B,  C,  and  D  in  Figure  8.  The  number  of 
coefficients,  number  of  data  cards,  and  values  of  coefficients 
can  be  determined  in  a  manner  similar  to  part  h  above.  Each  data 
card  will  contain  one  coefficient  for  each  polynomial.  The  data 
card  Format  will  be  4E14.7. 

l.  NGO  -  14  Format- -no  geometry  option.  For  NGO  ■  0, 
the  program  solves  the  problem  according  to  the  model  develop¬ 
ment  and  includes  explosive  stack  geometry  effects.  For  NGO  »  1, 
the  program  solution  excludes  explosive  stack  geometry  effects. 

m.  NBO  -  14  Format — no  barricade  option.''  For  NBO  *  0, 
the  program  solves  the  problem  according  to  the  model  develop¬ 
ment  and  includes  barricade  effects.  For  NBO  *  1 ,  the  program 
solution  excludes  barricade  effects. 

n.  IBO  -  14  Format  — isobar  option.  For  I  BO  ■  1,  the 
program  computations  exclude  the  isobar  option.  For  IBO  -  0, 
the  computations  include  the  isobar  option  which  also  requires 
that  NTHETA  *  -1  (See  part  g  above) . 

o.  ,  NEFPRD  -  14  Format — new  edge  to  face  pressure  ratio 
data  option.  For  NEFPRD  -  0,  the  computer  code  uses  internal 
polynomial  coefficients  to  account  for  stack  geometry  effect 
on  the  ratio  of  pressure  along  a  line  through  the  comer  where 
the  vertical  faces  of  the  bomb  stack  meet  and  the  center  of 
the  bomb  stack  to  pressure  along  a  line  perpendicular  to  the 
vertical  faces  and  through  the  stack  center.  For  NEFPRD  >  0,  the 
computer  code  will  read  new  polynomial  coefficients  as  in  h  above. 
The  number  of  coefficients,  number  of  data  cards,  and  values  of 


coefficients  can  be  determined  in  a  manner  similar  to  part  h 
above.  Each  data  card  will  contain  one  coefficient  in  an 
E14*7  Format. 

p.  NEFIRD  -  14  Format — new  edge  to  face  impulse  ratio 
data  option.  For  NEFIRD  ■  0,  the  computer  code  uses  internal 
polynomial  coefficients  to  account  for  stack  geometry  effect  on 
the  ratio  of  impulse  along  a  line  through  the  corner  where  the 
vertical  faces  of  the  bomb  stack  meet  and  the  center  of  the  bomb 
stack  to  impulse  along  a  line  perpendicular  to  the  vertical  faces 
and  through  the  stack  center.  For  NEFIRD  >  0,  the  computer  code 
will  read  new  polynomial  coefficients  as  in  h  above.  The  number 
o;  coefficients,  number  of  data  cards,  and  values  of  coefficients 
can  be  determined  in  a  manner  similar  to  part  h  above.  Each  data 
card  will  contain  one  coefficient  in  an  E14.7  Format. 


6. 


OUTPUT  FORMAT 


The  output  from  the  Blast  Program  depends  on  the  options 
that  are  chosen  on  the  input  control  card.  The  output  will  be 
as  follows: 

a.  Geometry  effect 

(1)  If  the  geometric  effects  of  ^he  stack  SL  X  8D  are 
to  be  considered,  the  program  will  print  out 

GEOMETRY  OPTION  USED 

STACK  LENGTH  -  SL  STACK  DEPTH  -  SD 

(2)  If  the  geometric  effect  of  the  stack  is  not  con¬ 
sidered,  the  program  will  solve  the  remainder  of  the  prob¬ 
lem  for  a  bare  hemispherical  charge  and  print  out 

GEOMETRY  OPTION  NOT  USED 

b.  Barricade  effect 

(1)  If  the  effect  of  the  presence  of  a  barricade  BL  X  BD 
around  the  stack  is  to  be  considered,  the  program  will 

. print  out 

BARRICADE  OPTION  USED 

BARRICADE  LENGTH  -  BL  BARRICADE  DEPTH  -  BD 

(2)  If  the  problem  is  to  be  solved  without  a  barricade, 
the  program  will  print  out 

BARRICADE  OPTION  NOT  USED 


c.  Charge  Weight 

Th«  equivalent  amount  of  tot  (i.e.,  „  ib,.)  that  j, 
repreiented  in  th.  .tack  of  ia  printad  Qut  aj 

OUAMTITT  CP  EXPLOSIVES  IN  STACK  -  m  LBS  OF  TNT 
a.  Preaaura  and  Inpulaa  Diatributiona 

(1)  The  azimuth  angla  9  ia  firat  printad 
AT  THETA  -  0  DEGREES 

(2)  The  following  headings  are  then  listed: 

2  OVERPRESSURE  SwmEUH  SUSlSE 

(3)  undar  the.,  title.,  the  following  counted  v.luaa  are 

printad.  (a,  .clad  di.tance,  (b)  praaaura  from  a  bar. 
hamiapherioal  char,.,  (c)  aatimatad  pra.aur.  from  th. 
actual  charg.  (ahapad  and/or  barricaded) ,  (d)  poaitiva  im- 

pule.  from  a  bar.  hami.pharioal  charga,  (.)  a.ti«.t.d  im- 
pulae  from  th.  actual  charga,  and  (f)  diatanc.  in  faat. 

«•  Pressure  Isobars 

U)  If  tha  iaobar  option  ha.  ba.n  uaad  th.  program  computa. 
the  acalad  di.tance  to  th.  iaobar  pra.aur..  i„  incr.mant. 

Of  5  degrae.  batwaan  the  azi«.th  .ngle,  of  0  and  no  degree,, 
scaled  diatanc.  will  b.  aolvad  for  iaobar  pra.aura.  of  1000, 
800,  600,  400,  200,  100,  80,  60,  40,  20,  10,  9,  8,  7,  6,  5, 

3,  2,  1,  0.5,  and  0.1  psi. 
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The  first  line  to  be  printed  under  this  option  is 
THE  ISOBAR  COORDINATES  FOR  AN  OVERPRESSURE  OF  1000  ARE 

(2)  The  next  line 

(PNR  -  PRESSURE  NOT  REACHED) 

defines  the  symbol  which  is  used  instead  of  the  radius  when 
the  pressure  given  in  Part  (1)  is  not  attained. 

(3)  The  following  headings  are  then  listed t 

THETA  Z 

(4)  under  these  titles,  computed  values  of  the  asimuth 
angle  and  scaled  distance  for  the  prescribed  pressure  are 
printed. .  If  the  prescribed  pressure  is  higher  than  any  of 
the  computed  values  for  the  particular  line  of  propagation, 
then  the  symbol  PNR  is  printed  under  X. 

The  sequence  of  output  given  by  Parts  (1)  to  (4)  is 
repeated  for  each  of  the  pressures  listed  in  Part  (1). 


t 


APPENDIX  II 


FRAGMENTATION  PROGRAM 

)  •* 

1.  FORTRAN  PROGRAM  DESCRIPTION 

This  computer  cods  corresponds  to  the  analytical  model 
described  in  Section  III. 

The  first  step  is  to  define  the  problem.  A  control  card 
is  read  to  instruct  the  computer  to  retain  previous  sets  of 
data  for  barricade,  fragment,  or  fragment  dispersion  on  multiple 
runs  or  to  read  in  new  data  for  barricade,  fragment,  or  fragment 

dispersion.  The  control  card  also  instructs  the  computer  to 

.  ( 

(  (a)  solve  the  problem  with  or  without  a  barricade,  (b)  sets 

> 

.  the  increments  of  theta  (asimuth  angle)  and  beta  (angle  of 
departure) ,  (c)  solve  for  a  certain  angular  orientation  of 

bosbs  in  the  barricade,  and  (d)  sets  ths  height  of  the  canter  of  the  bomb 
stack. 

When  the  barricade  data  are  read  into  the  code,  the  input  is 
converted  to  radius  to  barricade,  R3(J),  and  height  of  barri¬ 
cade,  ZW(J) ,  for  each  angle  of  6,  6(J).  J  is  an  integer  where 
1  <  J  <  (360/ITH)+l.  ITH  is  the  increment  of  theta. 

Initial  conditions  for  the  fragments  are  read  into  the 
program  in  terms  of  polar  angle  measured  from  the  nose  of  the 
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bomb.  The  input  increments  of  polar  angle  must  correspond  to 
the  increments  of  azimuth  angle  for  the  barricade  (cylindrical) 
coordinates  (i.e.,  if  A®outpUt  ■  20  ,  values  of  fragment  mass, 

initial  velocity,  and  number  of  fragments  must  be  read  for 

•  •  •  •  • 
values  between  0  and  180  at  the  angles  0  ,  20  ,  40  . . . , 

180*). 

To  include  the  effect  of  a  multiple  number  of  bombs  existing 
in  the  stack,  a  correction  factor  is  included.  This  factor  has 
been  named  "EFNB*  for  "Effective  Number  of  Single  Bombs  in  Stack." 
EFNB  must  be  determined  experimentally  such  that  tho  fragment 
dispersion  pattern  for  the  "Effective  Number  of  Bombs"  is  the 
same  as  for  the  detonated  stack. 

If  the  current  problem  being  solved  has  experimental 
results  to  compare  with  the  computer  output,  a  permanent  record 
can  be  printed  with  the  output  by  using  Table  4  — Experimental 
Fragment  Dispersion  Data.  To  use  this  option,  NCD4  is  given 
a  number  on  the  Control  Card  corresponding  to  the  number  of  data 
cards  that  are  to  be  read  into  Table  4.  The  input  format  will 
be  discussed  in  Appendix  II-5.  No  computations  are  performed 
using  these  data.  The  information  is  read  in  and  immediately 
printed  out  to  establish  a  permanent  record. 

So  far  in  this  section  the  eaphasis  has  been  to  organise 
the  available  information  in  a  usable  form  for  the  subroutine 
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TRAJ3.  The  subroutine  in  this  program  computes  the  trajectory 


for  the  fragment  "launched"  at  tbs  azimuth  syle,  a, 
and  departure  angle,  6;  considers  properties  of  the  fragment 
and  forces  acting  on  it;  estimates  the  range,  velocity,  number 
of  fragments,  and  fragment  mass  at  the  impact  point.  Table  5; 
and  estimates  the  distribution  along  the  directions  of  propaga¬ 
tion,  Table  6. 

To  produce  the  distribution  pattern,  the  differences  in 
angle  for  the  directions  of  propagation  are  established  by  the 
value  of  ITH  read  from  the  control  card.  The  program  is  de¬ 
signed  to  solve  for  the  distribution  360  around  the  boob  stack 
in  increments  of  ITH.  TTA  is  the  name  given  to  the  angle  mea¬ 
sured  from  the  center  of  the  open  side  of  a  three-sided  barri¬ 
cade  to  the  radial  direction  of  the  trajectory.  TTA  will  take 

•  • 

on  the  values  0  £  TTA  <  360  and  will  increase  in  increments 
of  ITH.  NOTE;  This  convention  for  TTA  is  established  as  a 
standard.  Actually,  barricades  of  any  shape  that  can  be  approxi¬ 
mated  by  one  or  more  broken  lines  of  finitely  swny  linear  seg¬ 
ments  can  be  handled  by  reading  in  the  angles  and  distances  to 
the  ends  and  the  height  of  the  segments.  The  format  for  reading 
in  the  barricade  data  will  be  described  in  Appendix  1 1-5.  In 
general,  the  origin  for  measuring  theta  can  be  placed  in  any 
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position,  but  the  standard  has  been  chosen  to  produce  computer 
results  that  can  be  easily  compared. 

The  computer  code  will  solve  for  the  increments  of  TTA 
just  described  when  the  CO  3tatement  on  the  major  loop  of 
TRAJ3  is  as  follows: 

DO  1800  J  -  1,  JTMM1 

JTMMl  is  the  integer  value  of  360  divided  by  increment 
Theta,  ITH,  where  any  resulting  decimals  are  dropped. 

JTM"1  -  HI 

Example  1.  If  the  radial  directions  of  the  trajectory  are 

t 

described  in  20  increments,  ITH  ■  20,  the  program  must  go 
through  the  DO  loop  360/20  or  18  times. 

Example  2.  For  ITH  -  14  the  program  must  go  through  the  DO 
loop  360/14  *  25.7  or  25  times. 


If  the  distribution  pattern  is  desired  in  a  specific 
direction  6  it  is  necessary  to  compute  the  number  of  increments 
of  ITH  in  6,  rounding  the  result  to  the  nearest  integer  value. 


Por  future  discussion,  call  this  intagar  L.  Tha  Do  atataiiant 
in  tha  major  loop  of  TRAJ3  would  ba 

DO  1800  J  -  L,  L 

Exarapla  3.  If  tha  radial  distribution  ia  daairad  for  tha 
( .rection 

0  *  120*  and  ITH  -  20* 

120*/20*  -  6 
OR 

L  -  6 

Example  4.  If  0  -  131*  and  ITH  -  20* 

131*/20*  -  6.55 
OR 

L  -  7 

This  means  that  tha  computer  would  solve  tha  problem  for 
9  «  140*  instead  of  tha  angle  131*  desired.  To  avoid  this  dif¬ 
ference,  it  would  be  necessary  to  change  tha  increment  of  theta 
to  a  more  suitable  value  to  produce  an  angle  exactly  equal  or 

very  near  to  the  desired  angle.  Note  that  whan  tha  distribu- 

tion  in  a  direction  e  is  computed,  tha  area  effected  by  tha 

«*» 

fragments  in  that  trajectory  is  between  the  angles  0  -  ITH/2. 
and  0  +  is H/?. 
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For  ouch  .aimuth  angl.  m  to  b.  ^  m  progr„  „lu 
•olv.  for  .  .aria.  0f  d.partur.  angl..  fran  b,rrlo,d.( 

BETA.  The  values  of  BETA  will  v-*,  . 

beta  Will  be  between  the  angles  -89«  and 

+”‘  *inee  +90*  “d  -90‘  ««  no  intaraat  and  o.»a.  probl.m. 
with  aom.  of  th.  trigomatrie  funotiona.  Th.  di.tinet  valua, 

of  BBT*  will  b.  datarminad  by  th.  incramant  of  bat.,  IB,  on 
th.  control  card.  *  larg.  valua  for  18  will  produo.  ooara. 

approximations  p  but  reauire  fa  wav  nn.n  *.  .  , 

require  fewer  computations  than  a  small 

value  for  IB. 

“»pla  5.  For  IB  -  20*  tha  program  will  computa  JBMX  traj.c- 

toria.  for  a.ch  angl.  tt*.  jbmx  io  computad  in  th.  program  .. 
follows : 

1)  Th.  numbar  of  inoramanta  of  ib  in  89*  dropping  any 
r**ulting  decimal  values  is  NB. 

NB  -  89*/20*  -  4 

2)  Th.  maximum  numbar  of  inoramanta  of  BBT*  that  will  b. 
aolvad  batwaan  -89*  and  +89*  ia  JBMX. 

JBMX  -  2  x  NB  +  1 


For  ..oh  valua  of  tt*  the  tr.J.otory  will  b.  computad  for 
BBT*  -  “80* ,  -80.,  -40*,  .20*,  0*,  20*,  40*,  80*  and  80*" 
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Example  6.  For  ZB  ■  10* 


1) 

NB  - 

89*/10  -8 

2) 

JBMX 

-2x8  +  1 

-  17 

FOr  etch  value  of  TTA  the  trajectory  will  be  confuted  for 
BETA  *  -85*,  -80*,  -75*,  ...,  +80*,  and  +85*. 

To  change  the  program  to  compute  trajectories  for  BETA 
from  0*  to  89*  would  require  the  following  card  changes 

Statement  No. 

1640-12  JBMX  -  NB  +  1 

1640-1  BETA  «  0. 

Increments  of  beta,  IB,  would  remain  as  previously  discussed. 

To  compute  a  trajectory  the  governing  parameters  for  the 
fragment  must  be  set.  Since  the  fragment  parameters  are  in 
terms  of  polar  bomb  coordinates,  the  angle  between  the  initial 
direction  of  the  fragment  and  the  bomb  centerline  out  the  nose 
of  the  bomb  must  be  computed.  This  angle  has  been  named  GAMMA 
and  will  be  a  function  of  TTA,  BETA,  and  the  angle  to  the  center- 
line  out  tne  nose  of  the  bomb,  ANGD.  This  computation  appears 
in  the  program  as  statement  1640  +  17. 


1640  +  17  GAMMA  -  (ASIN((SBS  +  CBS  *  SNS)  **.5))  *57.3 
where 

SBS  -  SIN2  (BETA) 

CBS  -  COS2  (BETA) 

SNS  -  SIN2  (THETA-ANGD) 

Since  GAMMA  can  take  on  any  value  between  0*  and  180*  and  the 
fragment  parameters  are  read  in  at  discrete  angles *  the  values 
at  the  discrete  angle  nearest  to  GAMMA  are  applied  as  the 
conditions  at  GAMMA. 

Each  trajectory  is  computed  as  a  broken  line  of  finitely 
many  linear  segments  as  indicated  by  Equations  (26a)  to  (26c) . 

The  initial  conditions  at  the  beginning  of  each  linear  segment 
are  the  initial  conditions  at  the  bomb  stack  or  are  the  end 
conditions  from  the  previous  segment. 

The  difference  in  angle  a  between  adjacent  segments  is 
DEL.  The  value  of  DEL  is  computed  in  the  program  as  a  function 
of  position*  velocity*  and  direction  of  fragment  velocity. 

The  initial  value  of  DEL  is  computed  so  that  the  horizontal 
distance  that  the  fragment  has  traveled  after  the  first  step  is 
slightly  beyond  the  barricade  wall.  This  appears  in  the  program 
at  Statement  No.  1640  +  21  as 

1640  +  21  DEL  -  G*l. 01  *DISTB/ (VZBR*VZER) 
where 

G  -  gravity  constant 

DISTB  ■  distance  to  barricade 

VZER  ■  intial  velocity  at  bomb  stack 
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Succeeding  values  of  DEL  ore  computed  in  such  a  way  that  either 
the  increase  in  range  or  change  in  height  is  50  ft.  depending 
on  whether  or  not  BETA  is  less  than  or  greater  than  45*  respec¬ 
tively.  In  the  computer  program  the  statements  appear  as 
follows:  To  increment  the  range  in  50  feet  increments  state¬ 
ment  1752  is  used  _  . 

1752  DEL  •  G  *  50./VSQ 

Check  if  absolute  value  of  BETA  is  greater  than  45* f  if  it  is, 
DEL  is  changed  by  the  two  statements  following  1725: 

1725  +  1  VAL1  -  ABS  (ALPH1) 

1725  +  2  (VAL1.GE. 0.7853)  DEL  -  ABS (DEL*COSALl) /SINAI) 

where 

ALP  HI  ■  the  angle  the  segment  of  the  trajectory  being 
considered  makes  with  the  horizontal. 

COSALl  -  COS  (ALPH1) 

jXSAl  -  SIN(ALPHl) 

Two  exceptions  exist  that  cause  DEL  to  take  on  different 
values  than  just  described.  The  first  exception  occurs  when 
the  value  of  DEL  that  is  computed  as  above  is  used  in  Statement 
No.  1646+4  and  the  resulting  velocity,  V2,  is  negative.  V2  is 
the  predicted  velocity  at  the  end  of  the  linear  segment  being 
considered.  In  this  case  DEL  is  repeatedly  cut  in  half  and  the 
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velocity  recalculated  until  V2  remains  positive.  This  check  and 
modification  is  done  with  Statement  Numbers  1649-3  to  1649-1. 

1649-3  IF (V2.GT.0.)  60  TO  1649 

1649-2  DEL  -  DEL/2. 

1649-1  GO  TO  1646 

The  second  exception  occurs  near  the  top  of  the  trajectory 
when  the  speed  of  the  fragment  is  less  than  100  feet/second. 

The  value  of  DEL  should  be  much  less  than  one.  Since  DEL  is  a 
function  of  one  divided  by  speed  squared,  small  speeds  can  pro¬ 
duce  values  for  del  on  the  order  of  one  hundred.  To  avoid  this 
possibility,  set  DEL  equal  to  0.1  for  speeds  less  than  100  fps 
to  establish  the  proper  order  of  magnitude  and  let  the  routine 
in  the  preceding  paragraph  make  any  further  modifications  that 
may  be  required.  This  adjustment  to  the  program  variable  DEL 
is  checked  and  modified  with  statements  1725-3  to  1725-1  as 
follows : 

1725-3  IF (VI. GT.  100.)  GO  TO  1725 

1725-2  DEL  -  .1 

1725-1  GO  TO  1645 

The  last  bit  of  information  required  before  computing  the 
trajectory  is  the  Ballistic  Coefficient  C  which  is  a  function 
of  the  fragment  mass  and  hence  it  must  be  calculated  for  each 
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trajectory  since  fragment  mass  varies  with  the  direction  of 
"launch".  C  is  computed  in  statement  1645-1. 

1645-1  C  ■  RHP*CD/(rM(K)  **.33333) 

where 

RHP  ■  a  combined  set  of  constants  for  a  fragment 

CD  ■  drag  coefficient 

FM(K)  -  fragment  masa  in  the  direction  of  launch 

To  indicate  if  the  top  of  the  barricade  haa  been  cleared 
or  not  cleared,  a  dummy  variable ,  NBCL,  is  used.  Initially  the 
value  of  NBCL  is  -1  indicating  that  the  barricade  has  not  been 
cleared.  If  there  ic  no  barricade  in  the  radial  direction  of 
the  trajectory  or  if  the  fragment  clears  the  top  of  the  barri¬ 
cade,  NBCL  is  set  equal  to  +1  and  the  trajectory  of  the  frag¬ 
ment  is  computed. 

The  trajectory  of  the  fragment  is  approximated  by  a  broken 
line  of  finitely  many  linear  segments  where  the  angle  between 
the  segment  and  a  horizontal  line  is  continuously  decreasing. 

The  initial  conditions  for  each  segment  are  known  from  the  term! 
nal  conditions  of  the  previous  segment  or  the  initial  conditions 
at  the  bomb.  The  initial  conditions  for  a  segment  include  range 
XI,  altitude,  Yl,  velocity,  VI,  and  angle  between  the  tra¬ 
jectory  segment  and  a  horizontal  line,  ALPH1.  The  terminal  con¬ 
ditions  include  range,'  X2,  altitude,  Y2,  velocity,  V2  and 
direction  of  terminal  velocity,  ALPH2.  The  terminal  conditions 


for  each  segment  are  computed  as  the  trajectory  is  mapped  until 
the  fragment  strikes  the  ground  inside  the  barricade ,  strikes 
the  barricade  wall,  or  clears  the  barricade  and  impacts  the 
ground  at  some  distance.  The  printed  output  from  the  program, 
discussed  in  detail  in  Appendix  II-6,  will  include  information 
to  indicate  that  the  fragments  do  not  clear  the  barricade  or 
will  record  the  impact  range  and  velocity  for  each  combination 
of  launch  angles  TTA  and  BETA. 

The  number  of  fragments  landing  at  each  impact  point  is 
assumed  to  be  the  same  as  the  number  initally  launched  in  the 
trajectories  by  a  single  bomb  multiplied  by  the  number  of 
effective  bombs,  EFNB.  The  number  of  fragments,  NUM(x) ,  in  a 
trajectory  will  be  dependent  on  the  angles  TTA  and  BETA  and  the. 
increments  of  each  of  these  angles  ITH  and  IB.  NUM(K)  can  be 
approximated  using  equation  III-36.  In  the  program  this  appears 
as  statement  1645-10. 

i 

1645-10  NUM  (K)  «FPS  (I GAMA)  *RITH*  (SIN  (RBTA2) )  -  SIN  (RBTAl) )  *EFNB 
where 

K  -  indicates  the  angle  of  departure,  BETA 

IGAMA  -  indicates  the  angle  GAMMA 

FPS  (  )  -  fragments  per  steradian 

PITH  -  increment  of  theta*  TTA,  in  radians 

RBTA2  -  radians  to  angle  BETA  +  IB/2 

RBTAl  -  radians  to  angle  BETA  -  ZB/2 

EFNB  -  effective  number  of  bombs  in  stack 


With  the  impact  conditions  known  along  ths  radial  direc- 
tions  from  tha  barricade,  tha  avaraga  distribution/aq.  ft.  can 
ba  estimated  by  tha  theory  of  Section  lll.4f.  Tha  data  along 

a  radial  direction  are  ordered  according  to  tha  increasing 
magnitude  of  impact  ranges  using  statements  1740  +  1  through 
1790.  Tha  assumed  impact  area,  tha  number  of  f ragments/sq .  ft. 
and  weight  of  f ragments/sq.  ft  are  computed  and  than  printed 
as  Table  6. 


2.0  FRAGMENTATION  PROGRAM  FLOW  DIAGRAM 


PRINT  HEADING  ON  TABLE  2 
BARRICADE  DATA 


KEEP 3  -  1 


PRINT  HEADING  ON  TABLE 
FRAGMENT  DATA 


DO  FOR  EACH  CARD  IN  TABLE  3 


READ  FRAGMENT  INITIAL  \ 
CONDITIONS 

TH1,  AFM1,  VOX,  FPS1 

INDEX  FRAGMENT  INITIAL 
CONDITIONS  IN  TERMS  OF 
OUTPUT  ANGLES 

PRINT  INDEXED  FRAGMENT \ 

INITIAL  CONDITIONS 

TH2,  AFM(I),  VO (I),  FPS(I) 


0 


8 


PRINT  HEADING  ON  TABLE  4 
EXPERIMENTAL  FRAGMENT 
DISPERSION  DATA 


l 


DO  FOR  EACH  CARD  IN  TABLE  4 


READ  EXPERIMENTAL  FRAGMENT^ 
DISPERSION  DATA  XC,  YC, 

XH,  YV,  NF1,  Wl,  UW1,  NF2 , 
W2 ,  UW2,  NF3,  W3,  UW3 


PRINT  EXPERIMENTAL  FRAGMENT  ^ 
DISPERSION  DATA  XC,  YC, 

XH,  YV,  NF1,  Wl,  OW1,  NF2, 

W2 ,  UW2,  NF3,  W3,  UW3 


CALL  SUBROUTINE 
TRAJ3 


READ  CD/  NEFB 


3. 


LIST  of  variables  to  computer  program 

MODEL 


for  analytical  fracmentat ion 


APM(J) 

ALPH1 

ALPH2 

ANG 

ANGD 

AREA 

BETA 

BETA1 ,  etc. 

C 

CBETA 

CBS 

CD 

COSAL1,  etc. 

D 

DB 

DELDG 

DELTA 

DENS 

DISTB 


Average  fragment  mats  emitted  by  bomb  at  an 
angle  GAMMA 

Direction  of  fragment  velocity  at  start  of  a 
segment  of  trajectory. 

Direction  of  fragment  velocity  at  end  of  a 
segment  cf  trajectory 

Number  of  radians  in  TH1. 

Equal  to  ©B  in  Section  in 

Representative  area  of  impact  tone  for  a  trajectory. 

Angle  between  horizontal  and  initial  velocity  of 
fragment. 

Angles  defining  region  containing  trajectory, 
used  to  calculate  number  of  steradians. 

Ballistic  coefficient 

Cos (BETA) 

Cos2 (BETA) 

Drag  coefficient 

COS (ALPH1) ,  etc. 

Distance  along  barricade 

Increment  BETA  divided  by  2. 

DEL/G 

Angle  between  R3  and  barricade  wall. 

Number  of  fragment. /sq.ft.  in  impact  zone 
Distance  to  barricade  from  center  of  bomb  stack. 
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DTH 

DTHB2 

Dl 

EFNB 

ETA 

PITH 

FLl,  etc. 
FM 

FPS(J) 

G 

GAMMA 

HB 

HCB 

HEITB 

IB 

IGAMA 

ITH 

JBMX 

JTHMX 
KANG 
KEEP  2 
KEEP  3 


Angle  between  ends  of  barricade  segment. 

Increment  of  THETA  divided  by  2. 

Intermediate  value  to  solve  for  SII  and  ALPHA. 

Number  of  single  bombs  to  produce  same  distri¬ 
bution  as  bomb  stack. 

Difference  between  barricade  and  bomb  coordinates. 
Floating  point  increment  of  THETA. 

Floating-point  number  of  TH1,  etc. 

Fragment  mass. 

Average  fragment/s teradian  in  direction  GAMMA, 
bomb  coordinates. 

Gravitational  constant. 

Angle  between  VO  and  bomb  centerline. 

Height  of  barricade  for  segment  defined  by 
TH1-TH2. 

Height  of  center  of  bomb  stack  about  ground. 

Height  of  barricade. 

Increment  of  departure  angle  BETA. 

Integer  number  of  increments  of  THETA  in  GAMMA. 

Increment  of  azimuth  angle  THETA. 

Number  of  increments  of  BETA  (IB)  between 
-89*  and  89f. 

Number  of  integer  values  of  ITH  in  360*. 

Integer  value  for  ANGD. 

IF  1/  Hold  Table  2  Data  from  previous  run. 

IF  1,  Hold  Table  3  Data  from  previous  run. 
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KEEP  4 


KFLAG 

KNG1 

KTH1 

KTH2 

NB 

NBCL 

NCD2 

NCD3 

NCD4 

NF1,  etc. 
NTH1 

NTH  2 

NUM(K) 

PHI 

RA1 

RA2 

RBETA 

RBTAl,  etc. 

RETA 

RHP 


IF  1,  holJi  Table  4  Data  from  previous  run. 

Indicator  to  initialize  R3  and  ZW  to  zero. 

Integer  number  of  increments  ITH  between 
barricade  and  bomb  coordinates. 

Integer  degrees  to  start  of  barricade  segment. 

Integer  degrees  to  end  of  barricade  segment. 

Max.  number  of  increments  of  BETA  in  89°. 

Index  to  indicate  when  barricade  has  been 
cleared. 

Number  of  cards  to  be  read  for  Table  2. 

Number  of  cards  to  be  read  for  Table  3. 

Number  of  cards  to  be  read  for  Table  4. 

Number  of  fragments  in  impact  area  1,  etc. 

Nearest  integer  number  of  increments  of  ITH  in 
TH1. 

Nearest  integer  number  of  increments  of  ITH  in 
TH2. 

Number  of  fragments  in  a  given  trajectory. 

Angle  between  barricade  and  R2. 

Radius  to  inner  bound  of  representative  impact 
area. 

Radius  to  outer  bound  of  representative  impact 
area. 

Number  of  radians  in  BETA. 

Number  of  radians  in  BETAi,  etc. 

Number  of  radians  in  ETA. 

Constants  for  computing  drag  force. 
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PIB 


RITH 

RKTHl 

RTH1 

RTH2 

R1 

R2 

P.3(J) 

SBETA 

SBS 

SINAI 

SII 

SN 

THl 

TH2 

TS 

TTA 

TTA1, 


Number  of  radians  in  increment  BETA. 

Number  of  radians  .in  ITH. 

,  etc.  Number  of  radians  in  THl,  etc. 

Radians  in  T*il. 

Radians  in  TH2. 

Radius  to  beginning  of  barricade  segment  being 
defined. 

Radius  to  end  of  barricade  segment  being  defined. 
Distance  to  barricade  at  angle  corresponding  to 
J  increments  of  azimuth  angle  THETA. 

Sin (BETA) 

Sin2 (BETA) 

Sin (ALPH1) 

Angle  between  barricade  and  Rl. 

Sin  (ETA) 

Angle  to  beginning  of  barricade  segment  being 
defined. 

Angle  to  end  of  barricade  segment  being  defined. 
Tan  (SI) 

Azimuth  angle,  THETA,  using  cylindrical  barri¬ 
cade  coordinates. 

etc.  Angles  defining  region  containing  trajectory, 
used  to  calculate  number  of  steradian*. 
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UF1,  etc. 

Maximum  fragment  weight,  grams,  etc. 

VAL1 

Absolute  value  of  ALPH1. 

VO(J) 

Average  fragment  velocity  in  direction  GAMMA. 

VSQ 

Velocity  squared. 

vsqd 

Velocity  squared  divided  by  DEL. 

W(I) 

Velocity  vector. 

VZER 

Initial  velocity. 

VI 

Velocity  at  start  of  a  segment  of  trajectory. 

V2 

Velocity  at  end  of  a  segment  of  trajectory. 

WTPA 

Weight  per  sq.  ft.  at  impact  areas. 

Wl ,  etc. 

Minimum  fragment  weight,  grams,  etc. 

xc 

X  coordinates  of  impact  area. 

XH 

WIDTH  of  impact  area. 

XV(I) 

Range  vector. 

XI,  X2 

Initial  and  end  conditions  for  a  segment 
trajectory. 

of 

YC 

Y  coordinates  of  impact  area. 

YV 

Length  of  impact  area. 

Yl,  Y2 

Initial  and  end  conditions  for  a  segment 
trajectory. 

of 

ZW(J) 

Height  of  barricade  at  distance  R3(J). 
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HHit 

mm 


11*1  ot»t 
101041 

iimi 

moki 

r  *  .<19 


>4  or  m  ’ 
0-  *  *  .  1 


Hlflit 

01001.1 


iimt 

mm 

1 1904 1 


nmi 

inm 

01101.1 

11101.1 


00101.1 

Him 


111041 

11111? 


n  ■y<nf.'» 

Hill*. 

111115 

mm 


«p4flP 


®006 oan  ftaom  <  input,  output  i 

1IHPNSJON  ®'iMC21l  / 

<'0NH0N  ?Mnf.1».  AfH<3*1»,  IflUlli  PPSt’Gl),  ®TC361>,  J.1M  . 
t  NNAX  .  HC0  »  ANGO  .  JTNNX  .  JINX  •  10  .  ITH  ,  •!!«  .  *IfH, 

*TT»t  .  TT \9  .THdlO)  ,  *Vlt*11  ,  mill)  .  HUH(lll)  ,  TNttH)  » 

»  TKiioi  .  txihii  , 


•rOHOT  *t*T''HFNTS 


511  F04HA1  (  1111  1 

105  P01HAT  (  IHIf  ?5X,  tOM 

111  fOOHAT  C  IT®,  F0.8,  fO.l 

51?  C03HAT  1  //TIN 


9 

y 

4 

0 

4 

T 

1 

4 

t 


35X,  hh 
45X .  ’OH 
’OH 
/10H 
>/11H 
’OH 
TOM 
4?H 


.?IX,  P5.1  1 


,  .  >• 

DISTANCE,  PT  ®A* 

/,  15X,  F5,t,  15X,  Pl.t. 


//  1 
1 

TA4|.e  1.  CONTOOt  OATA  ,  /» 

THCt  NUH4E0  t  /. 

?  T  l  ,  /. 

FflOO-OATA  OPTIONS  UPHOLD  .  11X t  3  (4X «T  1) 
KUO  GOODS  INPUT  TMI5  PJIOIILEH  ,11X,3C3X,T?> 
PiOftllOP  HT*  1  1«NC  I  ,?’X,  T3*  /t 

TNCPPHENTS  OP  TMfTA,  OEGOFES  ,?*X,  13,  /, 
INC°€H€HT0  Of  OFTA,  OFOOEP.S  ,?TX,  13,  /, 
ANGLE  1ST KEEN  COOPOINATfS  Of  1010  ANO  PANP 
OHTl'Hf.  TX,  F5.1,  /, 

TON  HEIGHT  Of  STACK  CE NTEO ’ 

OlO  fOOHAT  1  Of 10.1  1 

01*  fOOHAT  I  //TOH ’  TAPLF  ?.  OAOTACAOf  OITA 
t  S?H  TPETA,  CfOPEES 

?  POHNICIOE  Mftf.NT,  fT  ,  /  1 

51?  FO®MAT  (  15X,  P5.1,  15X,  ff.t,  1?N,  fl.l, 
t  l?x,  *5.1  ,  *  1 

4?«  FO»HAT  {  if H.1  I 

0?*  rOPHAT  (  //11H  T»r>LP  p,  FPAGNENT  OAT* 

t  5?H  TMEfA,  OEG 

P  -  ,t  -ItH  *®AGS  ,  f, 

3  5?H  1NPT  SOHO  COOTOl 

4  11ME°  STPOIOfAN  ,/  I 

5?7  fOOHAT  i  *X,  *10.1,  OX,  ffd.O,  IX,  flO.1,  3X,  flO.O  1 
5?A  POOHAT  C  /  1 

030  POPHAT  <  4F* .1,  3<  T5,  ff.i,  ff.l  )  1 

530  P01HAT  1  //«?H  TA"LP  i.  EXPPPTNENTAL  FtAGNfNT  OISPCOSXON  DATA 
t  /5?M  COORDINATES,  fT  NUN 

?  POHNETIHT  PANGE  ,/, 

1  5?H  XP  TC  XH  TV 

4  10M  (G«AH3),  /  1 

03*  POPHAT  <  IX,  41  »X,  f*.l  1,  3X,  13,  *X,  f*.t,  IX,  IN.,  IX,  f*. 

11,  /  »  ?!43x  ,  I?  ,  bt  ,  fl.l  ,  3H  -  ,  F5.1  ,  /  1  1 

534  POPMAT  (  1H1  1 


AVC  fPAG 


PASS,  G0AO3 


INITIAL 
V*L,  PP5 


p.....p»OG01P  ANO  PPOOLTO  lOENTTf TCATfON 

r 

5®1  HEAD  510,  »  OUNf II ,  T  »  1,  10  1 

P®INT  515,  f  OUN1I1,  X  «  1,  10  1 

r 

r, - TNPUT  T44LE  t  -  rONTOOL  OMA 

t* 

»p»1  511,  KPPP?,  KFEP3,  KEEP4,  NCI?,  Nil’,  NC04,  KfL AG, 

1  TTH,  10,  ANG1,  HC1 

KANG  *  ANSI 

If  <  t’M.EO.i  »  GO  ’0  044 

°®INT  51?,  KEEP?,  KEEPS »  XErP4,  NCO?,  NED?,  NC14,  KPL AG, 
1  ITH,  IS,  ANG1,  HCO 

join  ■  i 

JTHHX  •  (  3*1  )  ITH  1  *  1 
JINX  ■  1  100  /  TO  1  P  1 


011150 

111101 

1111*5 


191 


!1  <1*1*0 
'Ills'* 


wi69 
II'ICM 
111171 
1 V174 
101149 


111177 

"11747 

911794 

1117?! 


9H74S 

111749 

11*747 

111?5> 


11179** 

117799 


117797 

919769 

11179? 

11*7** 


1117S4 

911766 


111?€7 

11^771 

177?74 

7i"7*»<; 

7i«ms 

117?7f 
117*1 l 
717*19 


7TTH  a  TTM 
>I*M  ■  fTTH  /  97.11 

r, 

r . I*  <FHC  «  n  INITTAl.IZr  11  ANO  I«  TO  ?*#0 

r 


IF  (  <Fl»-,  .  £0  .  I  I  60  TO  997 
G 0  TO  994 

99?  99  941  T  *  1  ,  JTMN* 

>*111  *  1.1 

99*  7M(T»  •  1.1 

99/  IF  (  KF*P?  .  F.f)  .  t  )  GO  TO  6*? 
r 

r . TN"UT  TABLE  E  -  ,M*7tH0r  8ATA 

r 

r  — -  THETA  19  TMC  ftNGL.*’  NEA9UBFH  FMOH  THE  CFNTEMLINE  OF  THE  OPEN  SIOF 

<• - OF  t\  TM»*E  9IOE1  BAMMICAOC 

r 

C— —  4NGO  TS  TH«-  ANGLE  FMON  THETA  *  7**0  TO  TMF  CENTERLINE  Or  TH  BOMP. 
*  NM9UMF  ING1  IN  9AM*  OI»FCTTON  19  UUS  THETA.  ( 9.LE.«MC0.LE.*360  > 


r 


r 

C 

r 


r 

r. 

c 


4 

r 


"STMT  916 

90  671  1  *  1.  NCI)? 

9E"0  919,  TNI,  41*  THE,  ME*  MM 

®®IMT  917,  TM1*  »1,  MB,  THE,  •?,  MM 

(  i  .  *  t  '  r  •**  *  •  •  • 

CONVEPT  IMTFGER  VALUE9  IM  OFCREES  TO  NEAPEST  MULTIPLE  OF  ITM  ClN- 
9P5MFMT9  OF  POLAR  ANGLE  > 

\  .  ’  ,  •  -  H  -  ‘  I=r 

HTML  *  |  TM1  /  FXTM  )  7  6.9 

NTHE  *  f  THE  /  FTTH  I  ♦  0.9 

9T«1  *  MTMt  •  TTM 
KTM?  •  NTH*  •  TTM 

,  i  < 

90NVEMT  TN»<JT  ANGLES  TM1  AMO  TM?  INTO  RADIANS 

I  <  w  >  »*’ 

TTMt  *  TM1  /  97.10 
4THP  •  TM Z  f  fT.10 


CONVENT  ANGLES  MTMt ,  RTM7,  T()  t«0IANr 

v  r 


ELI  •  KTM1 
FL?  ■  NTH? 

”  PMTMt  ■  *Ll  / 
KTME  •  FL?  / 

90NVE»T  IMTrGF9  VALUES 
IM9FKING  PURPOSES 


»aui •-  in,  put-f 

*T.50 

97.57 

,  ,  ‘  *  •  * 

OF  ICTMI  AN9  MTM7  TO  INTEGER  INCREMENTS  FOR 


II  *  MTM1  7  l 
J?  *  NTM7  7  t 


COMPUTE  ANGLES  OF  T*I ANGLE  MO*  MT  Ml,  •?,  AMO  MAMMICAOF 


OTM  *  NTH?  -  MTM1 

SI  «  0.5  •  <  1.14  -  OTM  ) 

T9  *  TAM  (  SI  ) 

01  •  ATAN  C<<  *1  -P?  )  /  I  ME  ♦  Ml  J  )  •  TS  ) 
"HI  «  91701 
911  *  SI  -  01 

0  «  Ml  •  (  SINIOTMI  F  SIMCPMI)  > 

IMG  «  MXTM1 


9ET  HEIGHT  o*  SAPPTCAO*  EOM  EACH  INC*  OF  TMfTA  F»OM  JTMl-JTM? 
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11072(1  <50  615  J  *  II,  J2 

f.  —  .OOlPUTF  AN5 L€  MADE  XV  R»  ANO  9ARRI1AOr  MAIL 
/» 

0 1  ^ 2 2  OrLTA  «  1.14  *  <  PHt  ♦  f  RTH2  -  AN5  »  I 

r> 


S€T  MFICHT  Of"  14PRICA0F  FCR  PtRTANCP  p*IJ» 


111171 

HATH, 


100x7  6 
ni4ii 
o*»i4t4 
oi04i7 
<1114?*' 
100422 


»«<  11  =  *2  *  <  flNIPMT  I 

?m(  ji  ■  nr 


/  STN  I  DELTA  »  ) 


0- 

<* 


■CHPCK  F0»  1U°LIC4TE  VALUES  OF  J— TF  THr PF  ARE  TMO  VALUES  OF  J  f 
IMO  VALUES  Or  TH1  0»  TH’I  USE  VALUE  THAT  YIELDS  SMALLEST  R3CJ) 


1113X4 

IF  C ( J- IML01.5T.0.07. 

111345 

IF  <  PUJi.LT.HLOPI  1 

111Y«ft 

®XfJ»  «  MLORX 

1 11X5 I 

*M ( J|  «  HL17M 

11135* 

50  TO  614 

119?«i 

51* 

MLORX  s  R1(J» 

11"355 

ML02M  *  7U<Ji 

111357 

1NL1  »  J 

110351 

51 L 

ANG  «  AN5  ♦  PI*N 

11036* 

615 

'•ONTTN'I*' 

111165 

f* 

671 

CONTINUE 

.f AiOiif  rim  (T  i  _  raifyruf  i 

f* 

111367 

672 

tc  C  XEPP3.E0.1  )  GO 

1H341 

* 

p«TN*  97$ 

111374 

00  631  N  «  1,  N103 

SEAT  525,  GAMMA,  AFN1,  vCf,  FPS1 

N2  «  (  GANNA  /  FITH  )  ♦  1.5 

AFM(N2»  *  AFH1 

V1<N7|  «  V01 

PPSIM’I  *  FP51 

GAMMA  *  <  N2  -  1  I  •  TTM 


0  1*»4?5 
111441 


niH^i 

111/^0 

111451 
in  45? 

IIPSPJ 

o 11 9;  ; 

010547 
11  J55" 
1H551 
1H551 
"11554 


*»RTNT  5*7  ,  GAHMA  ,  AFM  C N 2 1  ,  VOCN*»  ,  FPSCN?) 

611  D0MTTMI1F 

c 

'•——INPUT  TABLE  4  -  EXPERIMENTAL  FRAGMENT  9T?PCRSI0N  DATA 

m 

651  TP  (  <FFP4,rl.  1  I  f,0  TO  661 
PRINT  9*9 

r 

DO  655  5*1,  NC04 

RE AO  531*  XC,  VI,  YM,  VV,  NF1,  Nl,  U«1 ,  NF2*  M2,  UM2,  NF1, 
1  MX,  UMX 

P»INT  5*6  »  XO,  V",  XM,  VV,  11*1,  Mt,  UMt ,  NF2,  M2,  UM2,  NF3, 

1  Ml,  UM3 

S«9  CONTTN'If 

r 

661  CALL  TPAJ3 

50  TO  1*1 
RR«*  EONTIN'IP 

"PINT  519 
TN9 


IP 
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f 


1)7*S1 


117*00 

117*5* 


117*6* 

1175S* 

117R65 

1175*5 

117555 


117*55 

*  17M6 
017*60 


007*61 


107572 
1 0  *57* 
007600 


hmgi? 


HHsm 

007607 

017610 

01761? 
H76i6 
117621 
0 1 7*  ? 1 

o 176?7 
0 ?76 16 


S'J0®0UYTNE  TRAJ1 

—  S'PROUTTNF  TO  '‘ALCOLATF  TRA IECTQPV  07  FRAGMENT* 

COMMON  ZM<T61),  AFM(161) «  ¥0(161).  Fps(lM),  71(160).  JMftX  • 

1  NHAX  ,  H»1  ,  »NC0  ,  JTMMX  »  JRHX  .  I*  ,  ITH  ,  RTTM  ,  FITM, 

7TTAt  .  TTA*  »7M( t 00)  ,  XV(tM)  ,  VV(IOI)  .  MUN(IOO)  ,  TM(IRO)  , 

7  TM(IOO)  ,  TX(tOI)  ,  TV (100) 


t«1/  70®m*T  (  7»0.2  ) 
1*11  FORMAT  (  //60M 
1  *?M 


*ARLr  *.  RANGE  AMO  VeL00tTV  OATA 


*?M  THETA  BETA 

20HUM  AVG  FRAG  MA**,G  t  / 


RANGE. FT  TNR  VEL.7P* 


t*?i  FORMAT  (  OK,  2(  *X,  7*.l  ),  SK,  FR.t,  6X,  77. t.  2X,  16.  2X.  77.1  ) 
1*10  70®mAT(0*.?(2X, 7*. 1I,16M  FRAGMENT  OOF*  NOT  CLEAR  9AR°ICA1F  ) 

t*61  FORMAT  (  ///////  .  RX,/w>FM(  ,  T1  ,1M)  »  ,  712.1  .  OMIG^AMS)  ) 

t*5fl  FORMAT  (  *K  ,  INC  ■  ,  €10.7  » 

1*61  FORMAT  (  *K  *  6MREf|  «,  F6.1  .  1M(0EGREFSI  ) 

RHP  I*  0  COMNINEO  *ET  OF  CONSTANT*  THAT  INCLUDES  «IP  OENSITV 

- IMP  *  .5  •  ®H0  •  ».10 

RMP  FOR  A  TUMSLTNG  CUBE  IS  ARRROXTHATELV  0.0000*5 

OEL  IS  INCRFHENT  FOR  ALPHA  (•AQTANS)  TN  OrTERHINING  TPAJECTORV 
0»  fragmfnt 


■SET  INITIAL  VALUE* 

G  «  12.2 

RMP  •  .0000*5 

0TMO2  «  ITH  /  2. 

CO  *  DRAG  COEFFICIENT  70*  7PAGMENT 

TO  GET  OPOIMAPV  FRICTION  FRFf  TRAJECTORY  SET  CO  *  0.0 

■EFNP  *  EFFPCTIVF  NUMSER  07  MONRS  IN  BOND  STACK 
•E7N*  «  t  MILL  SOLVE  FOR  TME  EFFECT  OF  A  SINGLE  ROM* 

READ  1506  .  CO  .  EFNB 


•VA"V  POLAR  ANGLc»  THETA 

NR  •  I*  /  TO 
J5MX  *  2  •  H*  ♦  1 
JHKM1  *  JTHMX  -  1 


...................  MA J0°  LOC*R - - - 

■  INCREMENT  A2IMUTH  ANGLE,  THETA,  FROM  0  TO  16*  OEGRFES 

00  1010  J  «  1  ,  JNXN1 


c  INITIAL  *£T A 
r 

1CTA  «  -NR  •  TR 
1660  1ISTR  «  »1(J) 

HFITO  *  7H(  1) 

r - SET  CONSTANTS 

RIR  «  If)  /  S7.3 
IT 0  «  (J-l)  •  ITH 
TTA1  *  T|A  -  1THR? 
TTA?  «  TTA  ♦  0THR2 
ETA  »  TJA  -  ANGO 
RFTA  «  ETA  t  17.7 


v 


n*6?s 

"i*f*i 

n»6** 


M7f*«: 
117*** 
"1*64" 
"1*64? 
il*64i 
"1 *646 


11 7*4? 
"1**57 


11**5* 

11465? 

11*664 


11*666 


11*66* 

11*67? 

11151/ 

11*6*5 

11*6** 

mini 


117*1? 

mil* 

11**14 

11**16 

11H15 


ii**?n 

titt?i 

10**?1 

11**11 

n**3i 

114136 

"177',? 

11*744 


55  «  SIN^'.M) 
555  *  «N*SN 
P1I5*  ill" 


C— ....... ............  HTMQ9  LOOP  —————  —  — - 

C— —  I5C*PMPNT  np»1irg®€  ANGLE ,  0«*TA,  F10M  -*1  TO  ♦  0r59EF« 

*• 

10  1740  5*1,  JIM* 

®1FT4  *  6£TA  /  54.? 

C*E*A  #  nimFMi 
SOFT*  »  SlNdlET*) 

COS  «  P4FTI  •  C3?T* 

S95  •  S"ET*  •  5151* 

<• 

C— 5*MH*  T5  *M?  ANGLE  5rTM«-fW  THE  INTTHL  VFLOCTTV  OF  THE 
0  FP46HFNT  AH')  THE  CEN*F*LTNE  0*  THE  4CH? 

GAON*  «(AS!0((  SIS  ♦  CIS  •  5NS  )  ••  .5  I  >•  57.1 
C  PPINT  15/1  ,  GANNA 

1641  F05MIT  (  tiv  ,  P11.4  » 
r 

f— — COMPUTE  N«f*9F5T  NUH1EP  “F  INCOpMFNTS  OF  THETA (ITH)  IN  G*HHA 
C 

IG«HA  *  GAOOA  /  FITH  +1.5 
C  ®*TNT  16*5  ,  IGAM  A 

1S<0  FORMAT  (  11Y  ,  l<i  i 
FNOO  «  ApO<IGAKA> 

C 

C— — SET  TNTTIAL  VELOCTTy  FOC  THE  FRAGMENT 

r. 

V*PP  *  VO (TGAOA) 

0 - -SET  OEL 

C 

OEL  *  f,  •  1.1t*OTST*»  ✓  CV*P®  •  V7”  ) 

01  *  TP  /  ?. 

OrTAl  *  flpTA-  0* 

0*TA?  »  «ET"f  10 
OOTAt  *  PFT«1  /  57.3 
POT*?  *  »ET*2  f  54.3 

r 

1—  COMPUTE  THE  NUM1EP  OF  FtAFNENTS  FJfCTFO  TN  OIOECTTON  OF  TPAJECTOPY 
— INCLUOE  THP  EFFECTIVE  NUM»FP  0*  MHOS 

c 

NUHCKI  »  ”5(1010*)  •  PIT)*  •  (5IN(*0TA?)  -  SIP  (10TA1H  *  fFNO 
*IPM1  «  POF*» 

VI  =  VZPP 
VI  *  1. 

VI  *  HOP 

#» 

0— POOL  «  -1  FRAGMENT  0*5  HOT  <*LEA®*0  6»»PTC»0r 

r - NOCL  »  *1  FPA5HFNT  H*S  PLEAPEO  1*PPIC*0E 

* 

POOL  *  -1 

IF  (  HriT*».LT.1.t  *  MlCt  «  ♦  ! 

r« 

COS*Lt  «  C?5(*LPH1  ) 

V50  *  VI  •  VI 
r 

r.  *  pop  *  cn  /  f  rp(K)  •*.31311  ) 

c  POINT  1551  ,  - 

1S/S  SINAI  «  STN(ALPHl) 

164®  1EL0G  *  9^1  *  5 

VST)  *  V51  •  1FLCC 
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i17Tf; 

AL°M7  *  *L3W1  -  1*1 

1 1’747 

01SAL7  *  '‘OMALPMEl 

i ITT*! 

V2  r  (  vi  *  OOSALt  -  4  •  VI  •  VSQO  1  /  CCSAL2 

oi?T5s 

IF  (  V7.GT.".  I  60  TO  1640 

11?760 

Oil  «  OEL  /  2. 

oo*?6i 

60  TO  14Z5 

10*7A7 

1640 

X*  =  X 1  ♦  V^OO 

1177*4 

r\ 

V?  =  VI  ♦  VSOO  •  SINAI  /  rOSALl 

r\ 

r. 

CHECK  TO  SEE  IF  PARTICLE  MAS  HIT  THE  GROUND  INSIOE 

a 

* 

HAS  STRUCK  THE  SARRTCAOE  MALL. 

117/47 

/* 

IF  CNOCL  »  1660*1700.1744 

117771 

1650 

IF  CV2)  1641*1670*1470 

417771 

1661 

®RINT  1074  ,  TTA  ,  SETA 

114441 

XVCK1  *  XI 

114044 

VV(K»  «  VI 

111404 

60  TO  1771 

114047 

16?1 

IF(9TSTO  -  Y?|  1640*1720*1771 

114017 

1611 

IF(MFITO-T?»  16O0(i  1664*  1560 

otiits 

1600 

NSCL  «  1 

111014 

50  TO  1720 

11111? 

1700 

TFCY21  1710,1720,1771 

111021 

1710 

"RtNT  1520  ,  TTA  ,  SETA  ,  VI  ,  VI  *  NUNOO  *  FNCKI 

111041 

XV  (Kl  *  XI 

110042 

vvm  «  vi 

111444 

60  TO  1 7*0 

110070 

1470 

•L®M1  *  AL°H? 

114444 

COSAL1  «  COSAL2 

111451 

VI  *  V2 

114001 

XI  *  XX 

111153 

u  «  v r 

111454 

VSO  *  VI  •  VI 

1110*0 

TF  (  V1.GT.100.  »  GO  TO  14*0 

11ifl41 

OEL  «  0.1 

114042 

60  TO  1645 

111052 

1725 

OEL  *  G  •  50.  /  VSO 

114064 

VAL1  *  4BSC4LRMU 

111064 

I*  (  VAL1.5r.7.7«51  »OEL*ASS(OEL*COSALl/SINAl) 

Til 474 

50  TO  1645 

111475 

1774 

OFT®  *  qpT4  ♦  IS 

010104 

1741 

CONTINUE 

HMf? 

CALL  0R0FO 

1111"1 

1410 

CONTINUE 

011106 

RETURN 

111146 

ENO 
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*S? 

,vV 


i  irl  ^  ’it-"w'-  - 


SU®»OUTTN"  0®,0P® 


C 

4®?4f'GF  <*4NSrS  TM  INOrr45JNS  OME®  Off  M1SMITUOF 

<• 

•i 

Ul?’9  COMMON  7W(7G1I.  AI-'HMGOI,  VOdSII,  FPSdfl).  93(3611,  JMiX  , 

1  NMAX  ,  HCM  ,  1MCO  •  JtMMX  ,  WMX  «  T6  ,  ITM  f  *>ITM  ,  FITM, 

7TTA1  ,  TT 42  ,FMU1fl)  *  *VC1*0I  #  VVCtOH  ,  MUKHOI  ♦  TNC110I  * 

3  TMCtlll  ,  TXU91I  ,  TVI131) 


110315 


M"334 

"11375 

111737 

ni7/i 

n«3<»? 

111365 

0113/6 

111155 

111757 

011751 

11H62 

111363 

111361 

01*371 

11137? 

11136/ 

*113*6 

11*615 

*l*/n< 

*11607 

11*611 

111/1’ 

111671 

*116*7 

11*631 

11*67? 

111673 

01*636 


111661 
*  1*667 
11166" 
11*666 
*1*651 
11*6*57 


11*657 


11*66* 

11166? 
11*67  7 
U15H 
011503 


*575  FO®M*T|/7*M  TAMLP  6.  0IST9I9UTT0N  04T4  ,/,5X,9MTHPT4 (II «  , 
1*5.’  ,  ?1M1CG®FF5  T«rTA(?|«  .  F5.7  ,  6MOCG9EES  ,  /  ♦ 

75X  «  7M9«M*.E  1  «  11M  ®INGE  7  «  51M  TM®  VEUF®*  HUM/SO  rT 
7  G9IMS/SQ  FT  I 

7576  F07M4T  (  2<  ,  2F11.1  ,  *13.1  .  ?(M  ,  E9.  1  II 
00  ?7<*0  HU  -  2  ,  |9MX 
NNM1  *  MN  -  1 
00  7350  H  ■*  1  •  NNMt 

IF  (  XVCMMI  -  VV(  Ml  >  7751  .  7711  ,  ®79l 
2750  MPl  *  N  ♦  1 

OO  2693  II  *  MPl  ,  MN 
ITM1  *  It  -  1 
TVIIII  «  WdlMil 
TXdll  *  KVdfMl) 

TMCIIl  «  MUM  dt  Ml  I 
7757  TMCIIl  «  PMdtMlI 
VVCNt  *  yVCMMI 
XVCMI  *  XVCMMI 
M'JMCMI  «  MUM(MM1 

pmcni  *  pmcmm) 
no  7756  I  S  N°1  «  NM 
Will  *  TV<T> 

XVCII  *  TX(TI 
M-IMCII  *  TN ( II 
7756  PMCII  «  TNCTI 
7691  CONTINUE 

P’TNT  7575  f  TTU  ,  T T|9 
941  «  0.0 

no  779 C  T  *  l  ,  I5MX 
IP1  «  I  ♦  1 

7«2  *  (  xy(TI  ♦  XVdPlI  )  t  2. 

IF  C  Ml. PI. 0.  I  GO  TO  7705 

C - CXLC'ILXTF  »’F1  of  L*N9INF  4®E»  4M0  MUHtTF  OF  F9AGMENTS 

C  ®E9  SO.  FT.  (  rt»A5  /  SO.  FT.  I 

49EA  *  XVdl  •  9ITM  •  1  PI?  -  911  I 
IFd.EO.JIMXI  GO  TO  779’ 

GO  TO  27«6 

979’  91?  *  7.  *  xym-9it 

•  9e1  *  XVdt  *  RI’M  *  {  o»7  -  911  » 

2796  IF  C  4PF4.LT.1.7  I  G0  TC  77«5 

C  COMPUTE  F91GMFMT5  /  10  *T  (  F94G/S0  FT  I 

OEMS  «  MUM (II  f  99EI 


C~ - COMPUTF  MEIGMT/Sn.  FT.  C  G9IMS/SO  FT  | 

C 

MT®  4  *  nFT?  •  FMCI1 

<> 

®9TNT  2576  ,  941  ,  ®*2  «  VVdl  ,  OEMS  •  MTPA 
7795  941  r  94? 

®796  COMTINU® 

9’TU9M 
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5. 


INPUT  FORMAT 


The  input  cards  to  a  problem  will  depend  on  the  problem 
being  run  and  the  number  of  variations  to  be  considered  in  a 
run.  The  input  to  a  single  run  or  the  first  data  set  of  a 
multiple  run  will  include  a  card  to  describe  the  problem,  a 
control  card  to  instruct  the  memory  to  read  data  for  the 
problem  and  set  parameters,  a  series  of  cards  containing  data 
to  define  the  barricade  geometry  surrounding  a  charge,  a  series 
of  cards  containing  fragment  initial  condition  in  terms  of 
polar  bomb  coordinates,  and  a  card  containing  the  typical  drag 
coefficient  for  a  bomb  fragment  and  the  number  of  effective 
bombs  in  the  stack. 

Subsequent  data  sets  on  multiple  runs  will  include  a  card 
to  describe  the  problem,  a  control  card  to  instruct  the  memory 
storage  to  retain  information  from  the  previous  run  or  to  read 
new  data  for  some  or  all  of  the  data.  If  the  barricade  geometry 
is  to  be  changed,  a  series  of  cards  will  be  included  to  define 
the  new  geometry;  otherwise,  these  cards  are  omitted,  if  the 
geometry  of  the  bomb  stack  is  to  be  changed,  a  series  of  cards 
containing  the  fragment  initial  conditions  in  terms  of  polar 
bomb  coordinates  will  be  included;  otherwise,  these  cards  are 
omitted.  A  card  containing  the  drag  coefficient  and  the  umber  of 
effective  bombs  will  be  included  in  all  data  sets. 
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a.  Description  Card 

Provisions  are  included  to  use  the  first  card  of  each 
data  set  to  describe  in  an  alphanumeric  format  the  problem 
being  solved.  The  field  for  this  format  will  be  80  spaces 
on  one  card. 

b.  Control  Card 

The  second  card  in  a  data  set  will  contain  information 
as  follows: 

(1)  KEEP2,  15  Format 

In  Columns  1-5  enter  a  zero  {0)  in  Column  5  if  barricade 
data  is  to  be  read  during  execution  of  the  current  problem. 
Enter  a  one  (1)  in  Column  5  if  barricade  data  is  to  be  re¬ 
tained  from  the  previous  problem  or  if  the  problem  is  to  be 
solved  without  a  barricade. 

(2)  KEEP3,  15  Format 

In  Columns  6-10  enter  a  zero  (0)  in  Column  10  if  the 
fragment  initial  conditions  sure  to  be  read  from  data  cards 
during  execution  of  the  current  problem.  Enter  a  one  (1) 
in  Column  10  if  fragment  initial  conditions  are  to  be  re¬ 
tained  from  the  previous  problem.  A  one  (1)  in  Column  10 
only  applies  on  multiple  runs  and  a  zero  (0)  will  be  necessary 
on  the  initial  problem  of  multiple  runs  and  all  single  problem 
runs. 

(3)  KEEP4,  15  Format 

In  Columns  11-15  enter  a  one  (1)  in  Column  15  if 
experimental  fragment  dispersion  data  is  to  be  retained 
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from  the  previous  problem.  Otherwise #  enter  a  zero  (0)  in 
Column  15.  A  one  (1)  in  Column  15  only  applies  on  multiple 
runs. 

(4)  102,  NCD3,  NCD4,  15  Format 

In  Columns  16-20/  21-25#  and  26-30#  enter  the  number 
of  data  cards  that  are  to  be  read  containing  barricade  data# 
fragment  initial  conditions#  and  experimental  dispersion 
data#  respectively. 

(5)  KFLAG,  15  Format 

In  Column  35  enter  a  zero  (0)  when  reading  new  barri¬ 
cade  data  or  to  solve  the  problem  without  a  barricade.  If# 
on  a  multiple  run#  barricade  data  from  the  previous  run  is 
to  be  retained#  enter  a  one  (1)  in  Column  35. 

(6)  ITH,  15  Format 

In  Column  36-40  enter  the  increment  size  for  the  azimuth 
angle#  TTA.  A  small  value  for  ITH  will  require  larger  num¬ 
bers  of  fragment  data  cards  to  be  read  and  more  computations 
than  large  values#  which  are  less  accurate  but  faster.  In 
general#  a  value  of  20*  for  ITH  will  produce  useful  results 
and  require  reasonable  computer  time. 

(7)  IB,  15  Format 

In  Columns  41-45  enter  the  increment  si  ze  for  the 
departure  angle#  BETA.  Small  values  for  IB  will  yield 
more  trajectories  and  impact  points  and  produce  a  more 
accurate  distribution#  but  more  computer  time  will  be  required. 
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In  general/  a  value  of  10°  for  IB  will  produce  useful  results 
and  require  reasonable  computer  time. 

(8)  ANGD,  F5.0  Format 

In  Columns  46-50/  enter  the  angle  between  TTA  *  zero 
(0)  and  the  centerline  out  the  nose  of  the  bombs  in  the  bomb 
stack.  ANGD  will  be  measured  CCW  and  lie  between  0°  and  +360°. 

c.  Barricade  Data,  5F10.1  Format 

Barricade  data  cards  need  to  be  included  on  single  runs, 
initial  problem  of  multiple  runs,  and  on  occasional  problems 
of  multiple  runs  when  it  is  desired  to  change  the  barricade 
configuration.  Each  card  will  contain  five  (5)  values  to  be 
read  in  a  5F5-1  Format.  The  data  on  the  card  will  be  the 
information  necessary  to  define  the  finite  number  of  linear 
segments  representing  the  barricade.  Each  segment  will  require 
one  (1)  card.  The  card  for  each  segment  will  contain  the 
following  information  in  the  order  presented. 

(1)  1HL,  F5.1  Format 

THl  is  the  angle,  degrees,  to  the  start  of  the  segment 
being  defined. 

(2)  Rl,  F5.1  Format 

Rl  is  the  distance,  ft.,  from  the  center  of  the  bomb 
stack  to  the  start  of  the  segment  being  defined. 

(3)  TH2,  F5.1  Format 

TH2  is  the  angle,  degrees,  to  the  end  of  the  segment 
being  defined. 

(4)  R2,  F5.1  Format 

R2  is  the  distance,  ft.,  from  the  center  of  the  bomb 
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stack  to  the  end  of  the  segment  being  defined. 

(5)  KB,  F5.1  Format 

HB  is  the  height,  ft.,  of  the  segment  being  defined. 

d.  Fragment  Initial  Conditions 

The  cards  described  in  this  section  are  require  ~  only  on 
single  runs,  initial  problems  of  multiple  runs  end  on  inter¬ 
mediate  problems  of  multiple  runs  when  the  fragment  data  is  to 
be  altered. 

One  data  card  will  be  required  for  each  polar  angle  where 
the  fragment  initial  conditions  are  to  be  read  in.  The  input 
angles  will  begin  with  0*  and  increase  in  increments  equal  to 
the  increment  of  TTA,  ITH,  previously  chosen.  The  number  of 
card  containing  fragment  initial  conditions  is  equal  to  HCD2. 
Each  card  will  contain  the  following  information! 

(1)  GAMMA,  F10.0  Fermat 

.Enter  in  Columns  1-10  the  polar  angle,  degrees,  measured 
in  spherical  bomb  coordinates  from  the  nose  of  the  bomb. 

(2)  AFM1,  F10.0  Format 

In  Columns  11-20  enter  the  average  value  over  the  region 
GAMMA  +  1/2  ITH,  for  the  fragment  mass,  grams,  in  the  direction 
GAMMA. 

(3)  VOl,  F10.0  Format 

In  Columns  21-30  enter  the  average  value,  over  the  region 
GAMMA  +  1/2  ITH,  for  the  fragment  velocity  ft. /sec.,  in  the 
direction  GAMMA. 


(4)  *PS1,  P10.0  Format 

In  Columns  31-40  enter  the  average  value,  over  the  region 
GAMMA  +  1/2  ITH,  for  the  number  of  fragments  per  steradian 
ejected  by  a  single  bomb  in  the  direction  GAMMA. 

e.  Experimental  Fragment  Dispersion  Data 

When  experimental  results  are  available  for  a  problem 
being  solved  by  the  program,  the  experimental  measurements  may 
be  read  into  the  program  to  produce  a  printed  record  for  com¬ 
parison  with  the  results  predicted  by  the  code.  There  will 
be  one  data  card  for  each  fragment  recovery  area.  The  infor¬ 
mation  on  each  data  card  will  be  as  follows: 

(1)  XC,  F6.1  Format  1 

Enter  in  Columns  1-7  the  distance,  ft.,  from  tht  center 
of  the  bomb  stack  to  the  center  of  the  recovery  area  measured 
parallel  to  the  direction  TTA  -  zero. 

(2)  YC,  F6.1  Format 

Enter  in  Columns  7-12  the  distance,  ft.,  from  the  center 
of  the  bomb  stack  to  the  center  of  the  recovery  area  measured 
parallel  to  the  direction  of  TTA  ■  45° . 

(3)  XH,  F6.1  Format 

Enter  in  Columns  13-18  the  dimension,  ft.,  of  the  recovery 
area  in  the  direction  of  XC. 
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(4)  YV,  F6.1  Forest 


Enter  in  Columns  19*24  the  dimension,  ft.,  of  the 
recovery  area  in  the  direction  of  YC. 

(5)  NF1,  NF2,  NF3,  IS  Fon»t 

Enter  in  Columns  25*29,40*44,  and  55*59  the  number  of 
fragments  recovered  in  the  recovery  area  in  the  heaviest 
weight  class,  second  heaviest  weight  class,  and  third 
heaviest  weight  class,  respectively. 

(6)  Wl,  W2,  W3  FS.l  Forest 

Enter  in  Columns  30*34,  45*49,  and  60*64  the 
fragment  weight,  grams,  in  heaviest  weight  class,  second 
heaviest  weight  class,  and  third  heaviest  weight  class, 
respectively. 

(7)  UW1,  UW2,  UN3  FS.l  Fonsst 

Enter  in  Columns  35-39,  50-54,  and  65*69  the  maximum 
fragment  weight,  grajus,  xn  the  heaviest  weight  class,  second 
heaviest  weight  class,  and  third  heaviest  weight  class, 
respectively.  More  than  three  weight  ranges  can  be  included 
by  changing  Formats  530  and  536  and  by  increasing  the  number 
of  variables  in  Statements  650+3  and  650+4. 

f.  Drag  Coefficient  and  Effective  Number  of  Bombs  in  Stack 
(CD  and  EFNB) 

(1)  CD,  F8.2  Foxmat 

In  Columns  1-8  enter  the  value  of  the  drag  coefficient 
to  be  read  in  a  F8.2  Format.  This  card  is  required  for  each 
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problem  of  a  multiple  run  and  obviously  for  a  single  run. 
Setting  drag  coefficient  equal  to  zero  will  solve  for  fric¬ 
tion  free  trajectory  In  general,  CD  will  have  a  value 
between  0.3  and  2.0. 

(2)  EFNB,  F8.2  Format 

In  Columns  9-16  enter  the  number  of  single  bombs  that 
must  be  exploded  one  at  a  time  to  produce  the  same  fragment 
dispersion  pattern  as  the  bomb  stack  being  modeled.  In 
general,  EFNB  will  be  less  than  the  actual  number  of  bombs 
in  the  stack. 

g.  Blank  Cards 

After  the  last  data  set,  include  two  blank  data  cards  to 
terminate  the  run. 

Notes  on  data  sets 

Normally  a  single  run  or  the  first  problem  of  a  multiple 
run  will  require  data  cards  a,  b,  c,  d,  and  f  as  a  minimum. 

The  inclusion  of  e  is  optional  on  any  problem.  These  data 
cards  will  describe  the  problem  and  solve  the  dispersion 
pattern  for  a  number  of  bombs  in  a  specific  barricade  configu¬ 
ration. 

If  on  the  second  problem  of  a  multiple  run  the  barricade 
geometry  is  to  be  changed  and  the  bomb  stack  remains  the  same 
as  the  first  problem,  data  cards  a,  b,  c,  and  f  would  be 
included  in  the  data  set. 
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If#  instead  of  changing  the  barricade  on  the  second 
problem#  the  type  of  bomb  is  to  be  changed#  data  cards  a# 
b,  c#  and  f  would  be  Included  in  the  data  set. 

h.  typical  Set  of  Input  Data 

The  data  shown  in  Figure  32  ere  for  a  single  run.  All 
data  cards  discussed  in  this  section  are  included  with  the 
exception  of  e. 


COLUMN  NUMBER 
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Figure  32.  Typical  Input  Date  for  Prograa 


6. 


OUTPUT  FORMAT 


a.  Description  of  Problem 

Across  the  top  of  the  first  page  of  printout  for  each 
problem  one  line  of  alphanumeric  print  is  used  to  describe 
the  problem  being  solved. 

b.  Table  1.  Control  Data. 

This  table  prints  a  record  of  the  variables  and  program 
options  used  for  each  problem. 

(1)  Table  1.  Heading 

Lines  lf  2  and  3  print  the  heading  for  the  table 
columns.  These  headings  appear  as  follows: 

TABLE  1.  CONTROL  DATA 

TABLE  NUMBER 

2  3  4 

(2)  Table  1.  Data 

Under  these  headings  are  printed  the  number  0  to 
indicate  that  new  data  are  read  into  Tables  2,  3,  or  4  an 
the  current  problem  or  a  1  to  indicate  that  data  from  the 
preceding  run  of  a  multiple  run  is  retained.  If  data  is 
read  into  Tables  2,  3,  and  4  on  the  current  problem,  the 
program  prints 

PRIOR  DATA  OPTIONS  (1  -  HOLD)  000 
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(3)  Table  1.  Hunter  of  Cards 

••hen  data  are  read  into  Table  2.  3,  or  4,  the  «*er 
of  crd.  read  into  aach  table  i,  printed.  „  3#  ^ 

7  data  card.  ar.  read  into  Table.  2,  3,  and  4,  re.pec- 
tively ,  the  program  will  print 

NOM  CAROS  INPUT  THIS  PROBLEM  , 

3  10  7 


(4)  Table  1.  Barricade  Bata 

If  Barricade  Data  are  not  read  into  the  current 

problem  or  to  either  omit  the  previon.  barricade  and  read 

n~  barricade  or  eolv.  for  the  dletribotion  without  . 

barricade  a  1  or  0,  reep.ctiv.ly,  i.  printed  with  the 
following  statement. 

BARRICADE  DATA  (l  «  NO) 

0  or  1 

(5)  Table  1.  Increment  Theta 

The  number  of  degree,  uawd  in  the  current  problem, 
for  the  increment,  on  azimuth  angle  theta  i.  printed  a. 

INCREMENTS  OP  THETA,  DEGREES 


(6)  Table  1.  Increment  Beta 

The  number  of  degree,  uaed  i„  the  current  problem 
for  the  increment,  on  departure  angle  bet*  i.  printed  a. 


INCREMENTS  OP  BETA,  DEGREES 
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(7)  Table  1.  Angle 

The  angle  measured  from  Theta  ■  0*  to  the  center* 
line  out  the  nose  of  the  bombs  in  the  stack  is  printed  as 

ANGLE  BETWEEN  COORDINATES  OF  BOMB  AND  BARRICADE  — 

(8)  Table  1.  Height 

The  height  of  the  geometric  center  of  the  bomb 
stack  is  printed  as 

HEIGHT  OF  STACK  CENTER 

c.  Table  2.  Barricade  Data 

This  table  is  printed  only  on  problems  where  barri¬ 
cade  data  is  read  from  data  cards .  If  data  are  retained 
from  the  previous  problem,  this  table  is  not  printed 

(1)  Table  2.  Heading 

The  heading  for  the  table  to  contain  data  defining 
the  geometry  of  the  barricade  is  printed  as 

TABLE  2.  BARRICADE  DATA 

THETA  DEGREES  DISTANCE,  FT.  BARRICADE  HEIGHT,  FT. 

(2)  Table  2.  Data 

The  data  used  to  define  a  segment  of  the  barricade 
geometry  are  printed  in  pairs  for  each  segment.  Data  for 
segment  1,  2,  ...,  N,  ...,  are  contained  in  lines  1  and  2 
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3  and  4,  “*,  2N-1  and  2N,  ***,  respectively.  The  first 
line  of  information  for  a  particular  segment  will  print 
the  angle  to  the  start  of  the  segment ,  the  radial  dis- 
"  tance  from  the  bomb  stack  geometric  center  to  the  start 
of  the  segment#  and  the  height  of  the  segment.  The 
second  line  of  information  for  a  particular  segment  will 
print  the  angle  to  the  end  of , the  segment/  the  radial  dis¬ 
tance  from  the  bomb  stack  geometric  center  to  the  end  of 
the  segment/  and  the  height  of  the  segment.  The  height 
of  any  segment  is  to  be  constant  between  its  end  points. 
When  segment  K  meets  segments  K+l,  the  printed  output  for 
the  first  line  of  segment  K+l  will  be  the  s me  as  the  last  line  of 
segment  K  with  the  possible  exception  of  the  height. 

d.  Table  3.  Fragment  Data 

This  table  is  printed  only  on  problems  where  the 
fragment  data  is  to  be  read  from  data  cards.  If  data 
is  retained  from  the  previous  problem,  this  table  is  not 
printed. 

(1)  Table  3.  Heading 

This  heading  for  the  table  to  contain  data  defining 
the  initial  condition  for  the  fragments  ejected  by  a  single 
bomb  as  a  function  of  polar  angle  measured  from  the  nose  of 
the  bomb,  is  printed  as: 

TABLE  3.  FRAGMENT  DATA 

THETA  ,DF.G  AVG  FRAG  INITIAL  FRAGS 

(WRT  BOMB  COORD)  MASS, GRAMS  VEL,FPS  PER  STERADIAN 


(2)  Table  2.  Data 


For  each  angle  where  date  are  input,  starting  with  0* 
and  increasing  in  steps  equal  to  the  increment  of  theta  to 
nJaxisnura  of  180*,  the  program  will  print  out  the  angles 
where  the  data  are  input.  The  fragment  mass  in  grams, 
initial  velocity  in  ft. /see.,  and  the  number  of  fragments 
per  steradian  are  printed  beside  the  appropriate  angle  to 
indicate  the  average  initial  condition  used  in  the  current 
problem  at  that  angle. 


e.  Table  4.  Experimental  Fragment  Dispersion  Data 

This  table  is  printed  on  problems  where  there  is 
experimental  fragment  dispersion  data  available  that  have 
been  read  into  the  program.  If  data  are  retained  from  the 
previous  problem,  this  table  is  not  printed. 

(1)  Table  4.  Heartily 

The  heading  for  Table  4  containing  experimental 
fragment  dispersion  data  is  printed  as: 

TABLE  4.  EXPERIMENTAL  FRAGMENT  DISPERSION  DATA 

COORDINATES,  FT.  HUM  WEIGHT  RANGE 

XC  YC  XH  YV  GRAMS 

(2)  Table  4.  Data 

For  each  recovery  area  of  dimensions  XH  by  TV  at 
coordinates  XC  and  YC  the  number  of  fragments  per  weight 
range  is  printed  for  the  three  heaviest  fragment  ranges. 
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f.  Table  5.  Range  and  Velocity  Data 

This  table  if*  printed  for  all  problems  run. 

(1)  Table  5.  Reading 

The  heading  for  Table  5,  Range  and  Velocity  Data, 
is  printed  as: 

TABLE  5.  RANGE  AND  VELOCITY  DATA 

THETA  BETA  RANGE, FT.  IMP  VEL,FPS  NUM  AVG  FRAG  MASS,G 

(2)  Table  5.  Data 

For  each  combination  of  azimuth  angle,  theta,  and 
departure  angle ,  beta,  the  range  to  the  impact  point,  ft., 
the  impact  velocity,  ft. /sec.,  number  of  fragments  expected 
to  impact  at  that  range,  and  the  average  fragment  mass, 
grams,  are  printed  in  tabular  form.  There  will  be  one 
Table  5  for  each  azimuth  angle,  theta,  where  trajectories 
are  confuted  and  one  entry  in  Table  5  for  each  value  of 
departure  angles,  beta. 

The  entries  for  any  single  line  of  Table  5  will  be 
in  one  of  two  forms. 

(a)  If  the  fragment  does  not  clear  the  top  of  the  barr¬ 
icade  vail  tor  the  launch  angles  €  and  B,  the  program 
will  print  the  launch  angles  and  indicate  that  the  frag¬ 
ment  does  not  clear  the  barricade.  This  appears  as: 

0  $  FRAGMENT  DOES  NOT  CLEAR  BARRICADE 

i 
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(b)  If  the  fragment  clears  the  barricade  or  if  there 
is  no  barricade  in  the  direction  6  for  the  launch  angles 
6  and  6 ,  the  program  will  compute  the  range  to  the  impact 
point,  XI,  impact  velocity,  VI,  number  of  fragments  ex¬ 
pected  to  impact,  NUM,  and  average  fragment  mass,  FM. 

The  printed  output  appears  as: 

6  0  XI  VI  NUM  FM 

g.  Table  6 

This  table  is  printed  for  all  problems  run. 

(1)  Table  6.  Heading 

For  each  azimuth  angle,  theta,  used  in  the  current 
problem  solution  a  Table  6  will  be  printed.  Each  Table  6 
will  give  the  distribution  in  a  direction  theta  for  an 
area  bounded  by  lines  symmet-ic  about  theta  and  separated 
by  the  angular  increment  of  theta,  ITH.  Let  6  be  the 
azimuth  angle  and  61,  62,  be  the  symmetric  lines  about  6. 
Then  the  Table  6  heading  will  be  printed  as: 

TABLE  6.  DISTRIBUTION  DATA 

THETA (1)  »  61  DEGREES  THETA (2)  -  62  DEGREES 

RANGE  1  RANGE  2  IMP  VEL,FPS  NUM/SQ  FT  GRAMS/SQ  FT 

(2)  Table  6.  Data 

Along  an  azimuth  angle  the  impact  ranges  are  computed 
at  dis^^ete  points.  From  these  discrete  impact  points 
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a  representative  impact  area  is  defined  by  the  angles  61, 
and  62  mentioned  in  g. (1)  and  average  ranges  between  impact 
points.  For  impact  point  number  1  call  these  ranges  RA1 
and  RA2.  Only  impact  areas  where  RA1  is  greater  than 
zero  (0)  are  of  interest.  Therefore,  for  the  first  im¬ 
pact  area  0  <  RA1  <  RA2.  In  the  printed  Table  6  this  would 
appear  as: 

RA1  RA2  VI  NPA  GPA 
where  PA1,  RA2 ,  NPA,  and  GPA  are  numbers  corresponding 
to  ranges  defining  the  representative  impact  area,  average 
impact  velocity,  number  of  fragments  per  sq.  ft.,  and 
number  of  gram,  of  fragment,  per  sq.  ft.  in  the  impact  area. 
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3  TEST  CASE  FOR  FRAGM  CHECKOUT 


table  i*  control  data 


PRIOR-DATA  OPTIONS  (UhOLD) 
nun  Caros  input  this  problem 

RARRICAOE  OATA  (  l.N0  | 

increments  of  theta*  oeorees 
increments  of  beta*  oeorees 

Hrl«MT8JJWX?,.CO0*0lNmS  0F  bomb 
heioht  or  stack  center 


TABLE  number 

2  3  4 

o  0  1 

3  in  o 

o 

20 

10 

ANO  BARRICADE  270 

♦•A 


TABLE  ?•  SaRIACAOC  data 

theta,  degrees  distance,  ft 


61.3 

132.0 


55.9 

67.3 


132.0 

223.0 


67.3 

67.3 


226.0 

296.6 


67,3 

55,9 


barricade  heibmt.  Ft 

U.o 

U.o 

U.o 

U.o 

.U.O 

u.o 


table  3.  fragment  data 


theta*  deo 

(HAT  BOMB  C00R0) 
0 

20 

40 

60 

BO 

100 

120 

140 

160 

ISO 


AVO  FRAG 
NASS*  SRAMS 


Initial 
vel*  fps 


fraos 

PER  steradian 


20 

17 

16 

10 

6 

9 

16 

23 

31 

AO 


AOOO 

ABBS 

5776 

666A 

7552 

7216 

5662 

AlOB 

255A 

1000 


5000 

AOOO 

3000 

3000 

•000 

liooo 

3000 

5000 

7000 

10000 


table  5.  Range  ano  velocity  data 

theta  beta  range, ft  imp  vel*fps 


0,  -80,0 
o.  -70.0 
-60.0 
0.  -50.0 

o.  -40.0 
0.  -30,0 

o.  -20.0 

U.  -10.0 
0.  0. 
o.  10,0 


0, 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

lloo.o 

1917.6 


7216.0 

7552.0 

7552.0 

7552.0 

7552.0 

7552.0 

7552.0 

7552.0 

♦87.8 

101.3 


HUM  AVO  FRAG  MASStG 


30929 

44289 

64739 

83322 

99177 

112119 

121655 

127*95 

129A61 

127A9S 


*.0 

6.0 

6.0 

6.0 

6.0 

6.0 

6.0 

6.0 

6.0 

6.0 
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3. 

20.0 

1970.6 

113.4  121655 

6.0 

«. 

30.0 

1«79. 4 

116.3  112119 

6.0 

0. 

40.0 

1699.7 

117.0  99177 

6*3 

«. 

50.0 

1462.4 

117.2  83222 

4.1 

0. 

60.0 

1177.5-  - 

UT.l _ 64739 

_ 6tQ _ 

0. 

70.0 

031.1 

117.3  44209 

0. 

80.0 

♦87.0 

125.5  30929 

9.0 

TABLE  6*  I 

distribution  data 

thftaiii* 

-lO.OOOOEOReES  THETA(2)»  10.80IDE6RCES 

RANSE  I 

RANGE  2 

IMP  VeL«PPS 

NUM/Sft  FT  _ 

NlNSaiPT 

244 

659 

125.5 

4.379E-01 

3.941(400 

659 

966 

117.3 

4.9I1E-01 

2.959(410 

966 

1139 

487.8 

1.944(450 

1.167(401 

1139 

1320 

117.2 

6.6I9C-01 

5.214(401 

1320 

1561 

117.2 

4.250E-O1 

3.750(400 

isai 

1790 

117.0 .... 

_ M1II41  .  -  _ 

4.9111*46  - 

1790 

1699 

116.3 

1.900(400 

9.489(400 

1699 

1944 

101.3 

4.177C«60 

2.586(401 

1944 

1997 

113.6 

3.339E400 

2.883(401 

TABLE  5.  RANGE  ANO  VELOCITY  DATA 


theta 

BETA 

RANGE. FT 

IMP  VCL.PP1 

HUM  AVO  PRA0  MASS.  6 

20.0 

•00.0 

0. 

7552.0 

22403 

6.6 

20.0 

-70.0 

0. 

7552.0 

44289 

6.0 

20.0 

-60.0 

0. 

7152.0 

44739 

6.6 

20.0 

•50.0 

0. 

7512.6  — 

- M - 

20.0 

•40.0 

0. 

7552.0 

9»l77 

6.0 

20.0 

•30.0 

0. 

7552.0 

112119 

6.6 

20.0 

-2U.0 

0. 

7552.0 

121655 

6.0 

20.0 

-10.0 

0. 

7552.0 

127495 

6.0 

20.0 

0. 

1150.0 

-  609.6 

48048 

10.0  - 

20.0 

10.0 

19 I 7*6  - 

—  -10  U3  — 

027495- 

- 2*0 - 

20.0 

20.0 

1970.6 

113.4 

121655 

6.0 

20.0 

30.0 

1079.4 

114.3 

112119 

6.0 

20.0 

40.0 

1699.7 

117.0 

90177 

6.6 

20.0 

50.0 

1462.4 

117.2 

03222 

6.0 

20.0 

60.0 

1177,9 

117.2 

44739 

-  6*6 

20.0 

70.0 

031*1 

H7.3 

A_* 

20.0 

00.0 

429.3 

117.3 

22493 

6.0 

table  6*  distribution  data 

theta (D*  IO.OOOOEQREES  TMETA<2) ■  co.ooodeqrees 


RAmGE  1 

RANGE  2 

IMP  VfL.PPS 

NUM/SO  FT 

0RANS/SQ  FT 

215 

630 

117.3 

3.613E-0I 

-2*162(400 —  - 

630 

991 

117.3 

4.237E-01 

2.542(400 

991 

1164 

609.6 

4.9T4E-01 

0.974(400 

1164 

1320 

117.2 

1.008E«00 

0.046(400 

1320 

1581 

117.2 

4.2501-01 

3.756(400 

1501 

1790 

117.0 

0.8195.01 

4,6111400 

1790 

1099 

112*3 _ 

— ~  _ _ 

•tifOf.ll 

1099 

1944 

101.3 

~  4.177(400 

2 .566(401 

1944 

1997 

113.4 

3.339C.00 

2.003(401 

table  5.  Range  ano 

VELOCITY  DATA 

theta 

SETA 

RANGE. FT  IMP-  VCUEP5  -  -MU&  -A2A. FBA6  MASS. 6 _ 

40.0 

-83.0 

0. 

7552.0  22493 

6.0 

40.0 

-70.0 

0. 

7552.0  44209 

0.0 
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,  r  $  T'A'*  * 

$  V  $  ^  -  y  „  ■»>  f.  T 


~*^V  ' 


4«.0 

-60.0 

40,0 

-50.0 

40.0 

-40.0 

40.0 

-30.0 

40.0 

-20.0 

40.0 

-10.0 

40.0 

0. 

40,0 

10.0 

40.0 

20.0 

40.0 

30.0 

40.0 

40.0 

40.0 

50.0 

40.0 

60.0 

40.0 

70.0 

40.0 

80.0 

TA%E  *. 
THFUCi), 
PAnGE  i 
215 
630 
691 
1164 
1433 
162.3 
2057 
2177 
2226 


0. 

7552.0 

64739 

0. 

6664 , 0 

31208 

0. 

6664.0 

37191 

o. 

6664.0 

42044 

0. 

6664.0 

45620 

®  • 

6664.0 

47810 

1150. 0 

689.3 

48548 

2199.1 

111.5 

47610 

2253,0 

122.6 

45620 

2155.0 

126.2 

42044 

1956.3 

127.3 

37191 

1649.0 

127.6 

3120B 

1177,9 

117.2 

64739 

831.1 

117.3 

44209 

*29.3 

117.3 

22*93 

distribution  data 
>  30«0000£UrcES 

RANGE  2 


630 

991 

1164 

1433 

1023 

2057 

2177 

2226 

2200 


XHP  VEL.rps 
117.3 

117.3 

609.3 

117.2 
127.6 

127.3 
126.2 

111.5 

122.6 


THETA  121 «  SO.OOODEGREES 
NUM/SQ  FT 
3.613E-01 
4.237E-01 
6.974E-01 
S.8S4E-01 
1.357E-01 
2.326E-01 
4.641E-01 
1.2ilE«00 
1.076E«00 


6.0 

10.0 

10.0 

10.0 

10.0 

10.0 

16.0 

10.0 

10.0 

10.0 

10.0 

10.0 

6.0 

6.0 

6.0 


GRAMS/SQ  FT 
2.166E*00 
2.542E400 
9.764Et00 
3.5l2EtOO 
1.357E*0C 
2.326E«00 
4.641E*00 
1.281E*01 
1.076E*01 


TAOLC  5.  RAM3E  ANO  VELOCITY  DATA 

THETA  BETA  RANGE. FT  IMP  VELtFPS 


60.0 
60.0 
60.0 
60.0 
60.0 
60.0 
60.0 
60.0 
60.0 
6  0.0 
60.0 
60.0 
60.0 
60.0 
60.0 
60.0 
60.0 


-00.0 

•70.0 

-60.0 

•50.0 

-40.0 

-30.0 

-20.0 

-10.0 

0. 

10.0 

20.0 

30.0 

40.0 

50.0 

60.0 

70.0 

60.0 


fragment 

fragment 

fragment 

fragment 

fragment 

fragment 

fragment 

fragment 

fragment 

2354.5 
2427,0 

2324.6 
2119.9 

1606.7 
1355.6 

"20.1 

379,4 


GOES 

DOES 

DOCS 

does 

docs 

DOES 

DOCS 

DOES 

DOES 


not  clear 
not  clear 
not  clear 
not  clear 
not  clear 
not  clear 
not  CLEAR 

not  clear 
not  clear 

117.2 

129.2 
133.0 

134.5 

127.6 
127.6 

117.3 
117.3 


TABLE  6. 
THETA  111 
range  1 
190 
600 
1000 
1521 
1903 
2222 
2340 
2391 


distribution  data 

so.ooooeurccs 


NJH  AVG  FRAG  HASSfG 

BARRICADE 

bARRlCAOE 

BARRIFAOE 

BARRlCAOE 

BARRICADE 

barricade 

BARRlCAOE 
BARRICADE 
BARRlCAOE 
47010  |4.o 

45620  14,0 

42044  14.0 

37191  14.0 

31208  10.0 

24277  10.0 

44289  6.0 

22493  6.0 


RANGE  2 
600 
1006 
1521 
1903 
2222 
2340 
2391 
2463 


imp  70.0000EBREES 

IMP  VELtFPS  NlJH/SQ  FT 
JJT* J  4.142E-01 

*17.3  3.170E-01 

Jfl’f  1.1B4E-01 

}*:••  1.3B7E-01 

J^4.5  1.576E-01 

}”.0  4.417E-01 

W136E*00 
I29»2  7.429E-01 


GRAMS/SQ  ft 
2.485E+00 
1 .9Q2E*00 
1 • 1S4£*00 
1.387E«00 
2.206E+00 
6.1S3E«00 
1.590E»01 
1.040E*01 
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APPENDIX  III 


CRATER  PROGRAM 

1.  FORTRAN  PROGRAM  DESCRIPTION 

The  program  on  cratering  follows  the  theory  sumearized 
in  Section  IV-4.  The  basic  values  for  the  parameters  that  form 
part  of  the  program  are  either  assumed  (tfL  -  500  tons,  ■  0.3, 

B  -  0.3)  or  result  from  the  detonation  of  a  100 -ton  hemispherical 
bare  charge  on  a  silty  clay  at  the  Suffield  Experimental  Station 
(Dm  «  21  ft.,  R  -  70  ft.,  soil  density  -  94  lb/ft3). 

New  reference  values  for  the  charge  weight,  crater  dim¬ 
ensions  and  soil  density  are  read  into  the  program  by  statement 
45.  The  following  statements  compute  new  values  for  the 
dissipation  ratio  and  the  ejecta  parameter  according  to  the 
formulas  given  by  Equations  (123),  (124),  (125a)  and  (125b). 

For  charge  weights  that  are  part  of  the  input  data,  the 
program  computes  apparent  depths  and  radii  according  to  Equations 
(117),  (119),  (120),  (121a)  and  (121b).  The  only  difficulty  in 
this  procedure  is  solving  the  transcendental  equation  (121a)  for 
the  ratio  R  /R°.  The  solution  is  obtained  by  means  of  the 
subroutine  TRANS.  Since  a  traditional  iterative  procedure  did 
not  work,  it  was  necessary  to  write  the  equation  in  an  alternate 

A 

form.  Suppose  we  define  a  function  FR(R)  by 

FR(R)  -  ^§pr  [x  -  V*’  (r4+?>]  *  R  u11-1’ 

8 


where 


(III-2) 


I 

Then  after  cone  algebraic  manipulations,  it  can  be  ehown  that 
Equation  (121a)  ie  eatisfied  if  pr(r)  is  zero. 

Eeeentially,  all  that  the  subroutine  TRANS  does  is  to 
find  that  value  of  fc  for  which  PR  is  zero  for  given  values  of 
K,  E®  and  c. 
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*4,v 


AND  APPARENT  DEPTH 

FOR  EACH  CHARGE  HEIGHT 

- 1 - 


LIST  OF  VARIABLES  IN  CRATER  PROGRAM 


a .  Main  Program 

Program  Variable 

BETH 

BETHST 

BETREF 

Bi 

BIST 

B2ST 

CDEN 

CDEN1 


Description 

a 

6 

$* 

Reference  value  for  $ 

6-2  _  8 
6  -  3  15 

A  ♦ 

6-2  8 

r* - 

6-3  15 


8  -  3 


DEN  Weight  density  of  earth  media,  pg 

(lb. /ft.3) 

DENREF  Reference  value  for  density  of  earth 

media  (lb. /ft.3) 

DENST  New  weight  density  of  earth  media, 

p *g  (lb./ft.3) 

DREF  Reference  value  for  apparent  depth  (ft.) 

DT  D  /D° 

A  A 
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gm 


Program  Variable 


EDS REF 

EDSST 

K 

NDATA 


Description 


DV(I) 

Apparent  dianeter 

DSA 

D°  (ft.) 

DZAST 

D*°  (ft.) 

Dl 

8 

“T? 

D2 

1 

~r 

D5 

CFo/CFo 

D6 

(1  -  E*°, 

/  (1  -  E 

D7 

l) 

1*5# 

*  a7 

EDS 

e° 

Reference  value  for  dissipation  ratio 
*D 


IF  NDATA  is  positive ,  new  reference 
values  are  to  be  read  in 
Nuaber  of  charge  weights  that  are 
to  be  read  in 


RATDST 

»‘X 

RATRST 

*>‘a 

RREF 

Reference  value 

for  apparent  radius  (ft.; 

RT 

V*. 

RV(I) 

Apparent  radius 

vector 

*  <y 

RZA 

R°  (ft.) 

k 

Program  Variable 


Description 


RZAST 

TZETI 

WL 

WREF 

WV(I) 

WZ 

WZST 

ZETA 

ZET1 

ZET4 

ZET4I 

b.  Subroutine  TRANS 
Subroutine  Variable 

A 

ALPHH 
ALPI1 
ALP4II 
ALP  4  4 
B 

DEL 

EDTT 

FNEW,  FN 
FP 


1/2 

WL  (tona) 

Reference  value  for  W 

c 

Charge  weight  vector 
Wo  (tona) 

WQ  (tona) 

C 

C  +  1 
c  +  4 
l/u  ♦  4) 

Deecrlption 

K/(E°)2 

C 

C  ♦  1 
1/U  +  4) 

C  +  4 

(1  -  eJ)/KC+1 
Absolute  value  of  FR(R) 

*D. 

New  value  of  FR(R) 
Previoua  value  of  FR(R) 
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Subroutine  Variable 


KK 

R 

RP 


New  value  of  R 


KNEW,  RH 


.jO<>n*i 

••coo**! 
i:  0  r  o  « 1 

•;0"n4i 


«0>)0M 


uOOrtM 


»0'*041 
II  Oi;  1)4  3 
00<)(l«4 
il0O04ri 
1)00147 

00005*) 

0000b4 

0000*1*) 
i»001C4 
OOOlrS 
OOUl'7 
0001 1 3 
'»00l-?l 

00U125 

000|13 

00014  i 
000144 
OOOUft 

.101115* 
000 151 
O0O15S 
uOulbfj 
O001<-.0 
00)161 
OOd  16.1 

(>0017.) 

>•001 1 1 
0001 1? 
(1001  1 4 

J0«l7S 
hhi  *r 

OOOPii" 

000*0*' 


c 

c 

r 

f 

r 


c 

r. 


»’R* -ram  Ctii\r’:  I  HPtjT  •  jilfPUT  ) 

.-1‘-  IS10M  V  (,>,), Rv<25>*  >v  <*5J 

ft«*L  K 

5  »  ?«  AT  (2110) 

1»  »  >'AT  IhH  «2» 

15  ♦  J-  .ftl  »lMl.l0AO4HHE»*f*t:>iCe  PA«4*n.Tl«S  for  THIS  run  // 

1  4«,  l7HHEf  EPf.  )CC  VlELl)  ■  fF7,*,5rt  TONS  *1SX» 

4*WL  «  •f-7*2*VI  Tf)*)S  //SXtiOHA^PARENl  DEPTH  a  , 

'  f4*2*4H  M  i  I  iA«17HA»»EA|<CNT  RA,)|MS  ■  »p6.2*3h  ►  T  // 

4  6A,l»HnlSNlMATI0N  RATIO  •  #F5»*t8X(6HZETA  a  f 
*■  f  ««2*4X*6'<'jtl(ji  *  iffctZ/// 

«•  ii*HAppA»ft it  Radius  a*o  depth  versu$  viclu  // 

/  »»*1HI#H**/ (•lT0MS)»ifJ*»  bMOEPTH  (9X*4HRAI)IUS1 

20  rQ'-AT  If  ID. V 14,2) 

RtTE-cfICE  Da7i  f  jn  SOIL  A  it)  RROcRa  '  CONSENTS. 

Da  1 4  H /bl).  ./,  yz/ioo./,  UZA/21 »/}  RZA/70./,  ZCTA/0.3/, 

A  t’S/^.Vt  oMM/3,1/*  oEN/q4./.  o  1  /0 .533333/,  02/3.J33333/ 


/KNa/fr**). 

/EI»*ZETA*4. 

/E>‘*Ial  ./<El4 
rZ‘ II«1«/I2*#Z» I  A) 

01  ia2#-<  (MZ/«LJ##D2> 

40  **t4  • 

If  <  lOAtA)  u»3  »45 

45  lf.w/ST,./i  ST*PZASr»OENST 

^Al  jSTaoZ<,ST/;>/„ 

''A^STaR?*'/ 

f  OS>T«F ns*  I  (*A  I .il*T#RATRS r ) ••T2ETI  > 

>(».•( 

OOa,l..PllSSr,/(l#.E0S) 

•)7*  ClFN/DtM'iT)  •  (RATRST*4j«) /RATi)S r 
..•l*MHtTH-2.)/(^cTH-3.1)  *ol 
i131bI)54D6*;74mj«r2ST/*</  ’ 
H2StaMi«;T4ni 

"L I  ’STa(3.«*.2Sl-2.)/(«2<.r-l,) 

*Rt*a«l2ST 
.f*-FF«ntNST 
.'•RE^.O/AST 
"'8fc>  ARZaST 
fcl)Wf «F0SS1 
>E  I -EFaHETMiT 
CUE  I*  >EF/.U**P2> 

I)  5«. 

SO  •Rt*aWZ 
LE  '><EF*oCl* 

I'Ht*  sD  it 
•  NHraR  Ik 
L  'ta 'FF»FnS 
E 1 *EF awr  Th 


55  COt  lafr.f..  #  *t> 
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mmm 


J002H- 

000220 

00‘J2?4 

000233 

000241 

000245 

ooo?si 
000272 
000311 
000 312 
000313 


RC'«  •  10, (4V<!) •lal«N«) 

i»0  70  I  a  1 ,  N  « 

('»«?•-(  tWV(l»/f<L|44n2>)«4V(I)/COENl 
<T, (HANS (4 F,ZET*»Zfcil ,ZET4»ZET4J) 
jT« (rT»*ZETi )/ (4<h»ZETA) 

KV(l)«RT*RRrF 

i'V(lK«0T*UREF 

**01\T  i5«WRe.F»«L.0RrF*RREF»E0SRCF*ZCTA*aeTREF 
PRINT  20,(I»0V(]i,fjv<ll»RV(l)*|«l  ,Ntf) 

00  JO  40  _  . 

CALL  EXIT 
END 


function  tra’J3(*k#edtt»alphh,alpii*alP44*alR4II) 


000677 

«E*U  KK 

000677 

FRH)»A#(  (l.-U4(P*nLP44l  )**2»0)-R 

000713 

4 

fo«s4t  aHO»io«*i6HiT£p4iiRN  failed  i__ 

000713 

A,\K/(tnTT«E0TT, 

000714 

4*<i ,-E0TT)/(^K44ALPll) 

00072C 

r* 

II» 

000721 

W 

HP  «  0, 

000721 

FP  «  A  _ 

000723 

HH.(l./M»**0LP4lI 

000726 

FN  ■  -RN 

000727 

5 

II*U*1 

000731 

IF  (20.11)  6,6,7 

000733 

6 

PR  I  T  4 

000742 

TK*NS*0, 

0007*3 

RETURN 

000743 

7 

RNE,  a  RP*FR»(HS.RP)/|FP.FNi 

00075(1 

F!)E-»aFR(RNE") 

0007*2 

OE^ABSIFNE*) 

000754 

IF  (UEL.0.0^2)  10,10,15 

000756 

10 

TRAxSaRfjErf 

0007,, 0 

HE  URN 

000760 

15 

IF  (FNE„)  14,14,20 

000762 

16 

HN«hNE* 

0007„3 

FN«FNE4 

000745 

GO  TO  5 

000766 

20 

MPaWNEW 

0007,,7 

F-ps**K|£  W 

000771 

GO  .0  5 

('00772 

C 

RClJRN 

n 00 774 

I  J4IIUIIUIPI 


5. 


INPUT  FORMAT 


Each  set  of  Input  Data  consists  of  either  two  or  three 
9roups  of  information.  The  integers  NDATA  and  NW  are  read  in 
on  the  first  card  according  to  the  format  2110.  If  NDATA  is 
a  positive  integer,  then  the  following  new  reference  values  are 
read  in  on  the  next  card: 

a.  WZST  -  the  new  reference  charge  W*  (tons) , 

b.  DZAST  -  the  new  reference  apparent  depth  D*°  (ft.), 

c.  RZAST  -  the  new  reference  apparent  radius  R*°  (ft.), 

d.  DENST  -  the  new  weight  density  of  the  earth  media, 

p*g  (lb. /ft.3). 

The  format  for  these  variables  is  4F10.2.  If  NDATA  is  negative 
or  zero,  this  card  is  omitted. 

The  last  group  of  input  data  consists  of  NW  charge  weights 
for  which  predicted  values  of  apparent  radii  and  depths  are  de¬ 
sired.  The  format  for  the  charge  weights  is  (8F10.2). 
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f 

6 .  OUTPUT  FORMAT 

The  output  format  is  illustrated  by  the  sample  runs  shown 
on  the  following  pages.  In  addition  to  the  reference  values  that 
formed  part  of  the  input,  values  for  the  dissipation  ratio 
and  the  ejecta  parameter  6  (beta)  are  given.  The  apparent 
depth  and  radius  for  each  of  the  charge  weights  are  shown  in 
the  table  following  the  reference  parameters. 

To  obtain  the  ejecta  depth  at  a  given  position,  0,  and 
the  apparent  depth  and  radius  are  used  in  connection  with  Figure 
28  of  Section  V. 


i 
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PEEEhl  mC-_  PAPA.fTCHs  FOR  THIS  WJ* 
»•£►»  WC  lCE  tCFui  •  ion..io  j  v|$ 


*i-  »  500.03  TONS 


*'»hapf»it  otPTi  *  21.00  rr 
•  >tSSIPAT!0.«  Nufio  S  ,j, 


APPARENT  HAD INS  ■  70.00  FT 
ft TA  «  ,30  HETA  a  l.lo 


APHAPf.  •  <?  »«.»0IUS  AN  »  >»TH  VEHSns  ^ItLD 


r 

I 

? 

3 

A 

5 

6 
7 
H 
9 

In 

U 


*  <T0>  S) 
10.  ..1 
PO.OO 
30.00 
ao.  on 
S0.no 
7j.no 

mn.ro 
«?0o.nn 
•oo. no 
>•00.00 
300.00 


depth 

14.7^ 

16.7#, 

I?. a* 

Id.  6b 
19.25 
20-11 
21*00 
22. b  6 

?3.5  / 
?A.lb 
2A.62 


H*f>I  JS 
32.7b 
A2.U 
*8.27 
52.3b 

56.  d  A 
63.11 
TO.  00 
8A.d3 
9A.  Jti 

100. rb 

105.bT 
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(Dletrlbutlon  Limitation  Statement  No.  2) 

^Analytical  models  and  aubaequant  computer  codea  have  been  developed  for  predicting 
peak  overpressure,  positive  unit  Impulse,  the  distribution  and  impact  velocity  of 
bomb  fragments,  efrater  dimensions  and  ejecta  thickness  from  the  detonations  of  typi¬ 
cal  bomb  stacks  used  by  the  Air  Force,  These  models  consider  aboveground  barricaded 
stacks  with  an  equivalent  net  weigh*;  high-explosive  range  of  10  to  500  tons  of  TNT. 
The  peak  overpressure  and  impulse  from  a  detonation  are  obtained  by  modifying  the 
known  results  of  a  bare  hemispherical  charge  to  take  into  account  the  stack  and 
barricade  geometries  and  the  interaction  effect  of  bombs.  Fragment  dispersion 
patterns  are  predicted  by  combining  experimental  results  for  single  bombs  and  using 
the  trajectory  equations  for  the  motion  of  a  steel  fragment  in  air.  By  using  basic 
principles  and  experimental  data,  crater  and  ejecta  dimensions  are  predicted.  Based 
on  output  from  tha  computer  codes,  illustrative  examples  are  given  together  with 
recommendations  for  future  tests  to  obtain  needed  data.  Programs  for  optimising 
munition  storage  areas  are  also  suggested. ( 
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